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Absorption curves have been obtained for cosmic rays at low altitude, with absorbers of 
carbon, iron, and lead. The separate intensities of fast mesotrons, slow mesotrons, and electrons 
have been deduced from the differences between the absorption curves. The number of slow 
mesotrons found is much greater than the number which may have descended from the meso- 
trons observed at a higher altitude, thus indicating considerable production of low energy 
mesotrons at low altitudes. The number of electrons is much too small to allow the conclusion 
that all of the energy of decaying mesotrons goes into shower production. Unless the lifetime-to- 
mass ratio for mesotrons is larger than presently accepted values, the results favor the con- 
clusion that less than half of the energy of the mesotrons goes to the electron component. 


INTRODUCTION 


EAR the base of the atmosphere or under 

any very heavy absorber, almost the entire 
electron component of cosmic rays must have its 
origin from decay and collision processes of 
mesotrons, chiefly from the decay. Thus the 
relative intensity of mesotrons and electrons at 
low altitude (or under a heavy absorber) may be 
used to furnish information regarding mesotron 
decay. This has been pointed out by Euler and 
Heisenberg,! who found the existing data in 
accord with a rough calculation based on the 
assumption that the mesotron decays into an 
electron and a neutrino with a lifetime of 2 
microseconds. More recently, Nelson? was able to 
conclude, on the basis of data taken by Neher 
and Stever,’ that the relative intensities were con- 
sistent with a lifetime of 2.8 microseconds, again 


and W. Heisenberg, Ergeb. d. exakt. Naturwiss. 
?E. Nelson, Phys. Rev. 58, 771 (1940). 
oth V. Neher and H. G. Stever, Phys. Rev. 58, 766 


providing that only half of the mesotron energy 
goes to the electron component. The relative in- 
tensities reported by Nielsen, Ryerson, Nordheim, 
and Morgan‘ have led Nordheim® to conclude 
that probably Jess than half of the mesotron 
energy goes into shower production. But this 
conclusion was based on a very small value of the 
mesotron lifetime (1.25 microseconds), and would 
have been different if a lifetime of 2 microseconds 
or more had been assumed. 

The definiteness of the above conclusions has 
been limited both by experimental difficulties in 
determining the ratio of electron and mesotron 
intensities, and by inaccuracies in the theories 
used to predict the number of electrons which 
should be observed as a result of mesotron decay 
and collision processes. Recent additions to these 
theories, made by Rossi and Greisen,* Rossi and 


*W. M. Nielsen, C. M. Ryerson, L. W. Nordheim, and 
K. Z. Morgan, Phys. Rev. 59, 547 (1941). 

5 L. W. Nordheim, Phys. Rev. 59, 554 (1941). 

* B. Rossi and K. Greisen, Phys. Rev. 61, 121 (1942). 
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Klapman,’ and Richards and Nordheim,® have 
made it possible to calculate more accurately the 
number of electrons which should be observed 
above certain energy limits. The present experi- 
ments were designed so as to eliminate some of 
the experimental uncertainty in the ratio of the 
two intensities. 

Most of the experimental difficulties in de- 
termining the relative intensity of mesotrons and 
electrons (from counting rates in Geiger-Miiller 
counters below different amounts of absorber) 
arise from the comparatively small value of the 
electron intensity. Since one must rely on small 
differences between large counting rates, the 
statistical errors as well as all small systematic 
errors become important. The effect of meteoro- 
logical changes (pressure, temperature, distri- 
bution of air mass) is great, and the number of 
low energy mesotrons stopped by the absorbers 
must be determined. The effect of side showers 
incident upon a counter telescope creates a sig- 
nificant error, and the transition effects observed 
as the absorber is changed become important. 
Moreover, the fraction of the electrons recorded 
depends strongly on the energy limit set by the 
apparatus, which must therefore be known. 


EXPERIMENTAL ARRANGEMENT 


The arrangement of counters used for the 
principal measurements is shown in Fig. 1. The 
fourteen counters arranged in a semicylinder were 


Fic. 1. Counter arrangement used in 
principal measurements. 


60 cm long and were all connected in parallel. 
The central counter was 20 cm long. The princi- 
pal measurements were of the coincidence rate 
between the central counter and the outer ones, 
as different absorbers were placed between the 
counters. Since the outer counters cover almost 


7 B. Rossi and S. J. Klapman, Phys. Rev. 61, 414 (1942). 

8 J. A. Richards, Jr., and L. W. Nordheim, Phys. Rev. 61, 
735 (1942). We are indebted to the authors for communi- 
cating their numerical results to us. 


the entire solid angle above the central counter, 
the rates obtained give an intensity integrated 
over all angles. This serves a twofold purpose: 
(1) The error due to side showers is eliminated, 
since we wish to record the rays traversing the 
central counter in all directions, and (2) the 
counting rate is high, so that a high statistical] 
accuracy may be obtained in a comparatively 
short time. 

The absorbers used were cylindrical in shape. 
The energy limits set for electrons by the smallest 
absorbers were determined with fair precision by 
constructing these absorbers of a material of low 
atomic number, carbon (for which the “‘strag- 
gling’”’ in the energy loss is small). Throughout 
the measurements with other absorbers (iron and 
lead), a central absorber of carbon was kept in 
place, so that only those electrons were recorded 
which emerged from the iron or lead with an 
energy of about 10’ ev. 

The errors in the counting rates, due to chance 
coincidences and inefficiency of the counters, 
were very small (a few tenths of a percent) and 
were determined by auxiliary measurements of 
the resolving time of the apparatus and the 
inefficiency of the counters. Other more serious 
sources of error are discussed in detail below. 


DEVIATIONS DUE TO METEOROLOGICAL 
CHANGES 


A simple correction for pressure and tempera- — 


ture changes cannot be made, since the cosmic- 
ray intensity depends on the distribution of air 
mass and temperature above the apparatus, 
more than on the values of pressure and tempera- 
ture at any one point. Therefore an attempt was 
made to obtain an empirical correction simul- 
taneously for all meteorological fluctuations, by 
using the deviations from the average of two 
“standard” rates, which were taken repeatedly 
during the running of the experiment, on alter- 
nate nights, for periods of more than 12 hours. 
These were the rate with the maximum lead 
absorber and the rate with only a small thickness 
of carbon as absorber. The statistical errors in 
the individual runs were less than 0.5 percent. 
The ‘‘standard” rates were utilized as follows. 
Since the rates had been taken so often, it could 
be assumed that the simple averages were the 
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Percent Deviations 


Pressure in mm Hg 


----- Curve from which corrections were token 


Borometric record 


t+ +t 
Date (noon) 


Fic. 2. Percent deviations due to meteorological changes. 


correct rates for the average environmental con- 
ditions. The percent deviations of the individual 
runs from the average rates were plotted against 
the mean time at which the individual rates were 
obtained. In spite of the fact that the rates with 
the maximum lead absorber included only meso- 
trons, while the rates with the small carbon ab- 
sorber included also the electron component, 
there was nothing to indicate that the two graphs 


Fic. 3. Illustrating errors due to showers and secondary 
particles, and methods of correction. 


of percent deviations against time should not be 
superposed ; therefore, the points taken with both 
absorbers were plotted together, giving a graph 
of percent deviation from the mean which should 
be valid for all the intermediate absorbers as well. 
The corrections for the rates taken with the other 
absorbers were read from this graph, according 
to the mean time at which the rates were ob- 
tained. From these corrected rates, the average 
rate was computed for each absorber used. 

In Fig. 2 we have plotted the percent deviations 
from the average, including the individual runs 
not only with the “standard” absorbers but with 
the intermediate absorbers as well. For compari- 
son, we have also drawn the barometric record 
for the same period of time. It is obvious from 
this graph that changes in cosmic-ray intensity 
are not always exactly ‘in phase’’ with corre- 
sponding changes in barometric pressure, and 
that the same change in barometric pressure does 
not always accompany the same change in 
cosmic-ray intensity. 


EFFECT OF INCIDENT SHOWERS AND SECOND- 
ARIES GENERATED IN THE ABSORBER 


Figure 3a illustrates the error associated with 
showers incident upon the apparatus, which leads 
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T = Shower effect pius collision effect- 
C = Collision effect for all particles 

i M = Collision effect for mesotrons 

S= Shower effect 


Counts per minute 
w 


1 


g/em2 Pe above carbon 


Fic. 4a. “Shower” and “collision” effects in iron. 


6 
T = Shower effect plus collision effect 
C = Collision effect for all particies 
sr M = Collision effect for mesotrons 
S = Shower effect 
t 
4 
2 
3 re 
AM 
= 
g/cm? Pb above carbon 
Fic. 4b. “Shower” and “‘collision”’ effects in lead. 
to an underestimate of the number of electrons caused by collision electrons knocked out of the 13 
striking the apparatus. Except for very small ab- absorber by mesotrons.® 42 
sorbers, this effect is unimportant, because there In order to determine the error due to the a 
is no error unless both particles have enough shower and collision effects, an auxiliary experi- ; 
energy to penetrate the absorber. Even for small ment was performed with counters arranged as in 1 
absorber thicknesses the effect is not very large, Fig. 3c. The two upper counters were connected } 
because of the great spread of air showers. For in parallel, and threefold coincidences between a 
brevity we shall refer to this effect as the ‘‘shower these and the lower counters were recorded as a 1 
effect.” function of absorber thickness, with the same 1 
Figure 3b illustrates an error in the opposite absorbers as were used in the principal measure- t 
direction, due to secondary particles generated in ments. For comparison, we also recorded the pe 
the absorber. This is a sort of ‘‘transition effect’”» —————— — 
which may be described as an increase in effective _°* Admittedly the term “collision effect” is not entirely 
appropriate, since a large part of this effect is caused by tw 
size of the central counter due to the presence of owers generated in the absorber. The ¢ only real distinction ac 
; H between the “shower” and “‘collision” effects, as we have cot 
the absorber. We shall ref er briefly to this effect called them, lies in whether the multiplicity existed above per 
as the ‘‘collision effect,’’ because it is largely _ the absorber or was caused by the absorber. to 
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Mesotrons 


Counts per minute 
8 


6s 
64 — 

momentum in 10% ev/c 
“a0 
0.6 os 10 | 


Fic. 5. Absorption curves in iron and lead. 


twofold coincidence rate between the upper 
counters and one of the lower counters. It can 
easily be shown that the threefold coincidence 
rate obtained in this experiment is proportional 
to the sum of the errors due to the shower effect 
and the collision effect."° The proportionality 


TABLE I. Corrections to counting rates, and final 


corrected rates. 
Average Net correc- 
rate cor- Net correc- tions for 
rected for tion for i i 
meteoro- “‘collisions” and 
logical chance Corrected 
Absorber changesonly “showers coincidences rate 
None 90.3 min.-! +0.3 min.“ +0.1 min.~ 90.7 min.“ 

1.36 g/cm?C 87.8 —0.5 +0.1 87.4 
2.84 g/cm*C 85.7 —1.1 +0.1 84.7 
4.2 g/cm*C 83.7 —1.5 +0.1 82.3 
6.1 g/cm*C 81.6 —1.8 +0.1 79.9 
2.84 g/cm?C 

+2.5g/cm*Fe 82.4 —2.0 +0.1 80.5 

+7.8g/cm*Fe 80.1 —3.1 +0.1 771 

+14.0g/cm*Fe 77.8 —3.9 +0.1 74.0 

+21.4g/cm*Fe 76.3 —4.2 +0.1 72.2 

+38.2g/cm*Fe 74.2 —4.3 +0.1 70.0 

+62.8g/cm*Fe 70.8 —3.8 +0.1 67.1 
2.84 g/cm*C 

+3.8g/cm*Pb 81.8 -3.1 +0.1 78.8 

+11.3g/cm*Pb 78.9 —4.8 +0.1 74.2 

+22.5g/cm*Pb 76.0 —4.6 +0.1 71.5 

+33.8g/cm*Pb 72.6 —4.0 +0.1 68.7 

+67.5g/cm*Pb 68.6 —2.5 +0.1 66.2 

+101 g/cm*Pb 66.7 —2.3 +0.1 64.5 


1” This conclusion rests on two assumptions: (1) That for 
two particles striking the apparatus simultaneously in such 
a direction that one of them penetrates one of the lower 
counters, it is equally likely that the second particle should 
penetrate the same counter or an identical counter adjacent 
to the first ; and (2) that particles which pass by the central 


constant is the ratio of the coincidence rate in the 
principal experiment to the twofold coincidence 
rate in this auxiliary experiment, which ratio had 
the same value (3.42), within statistical errors, 
for all of the absorbers used. 

In Figs. 4a and 4b we have plotted the results 
of this auxiliary experiment. The uppermost 
curve (7) represents the sum of the errors due to 
shower and collision effects, determined as ex- 
plained in the above paragraph. The correction to 
be applied to the counting rates in the principal 
experiment, however, is the difference between 
these twoerrors, rather than the sum. Fortunately 
the shower effect is small and decreases rapidly as 
the absorber thickness is increased. From the 
results obtained with very small absorbers in the 
experiment described above, it was possible to 
estimate the magnitude of this effect for zero 
absorber, while the variation of the effect with 
absorber thickness could be roughly calculated ; 
such estimates (admittedly not exact) have led to 
the curves S, the lowest dashed curves in Figs. 4a 
and 4b. Curve C represents the difference be- 
tween the ordinates of T and S; i.e., the collision 
effect alone. The corrections applied to the 


counter at distances greater than a counter diameter 
(4 cm) do not produce secondaries which discharge the 
counter. The first assumption is justified by the large 
spread of air showers; the second by the small average 
spread of showers in iron or lead. 


64 
| 
80 
| 
the ; 
the 
eri- 
S in q 
ted 
sa 
me 
ire- 
the 
by 
tion 
ave 


counting rates are given by the differences be- 
tween the ordinates of C and S. The curves M in 
Figs. 4a and 4b represent the collision effect for 
mesotrons alone, as determined with the arrange- 
ment shown in Fig. 3d. 

In Table I we list the absorber thicknesses used 
in this experiment, and the corresponding count- 
ing rates, together with the corrections applied to 
the rates. 


-SEPARATION OF MESOTRON AND 
ELECTRON INTENSITIES 


The method of obtaining the separate intensi- 
ties of mesotrons and electrons from the counting 
rates resembles a method recently used by Auger" 
in analyzing measurements taken at a higher alti- 
tude. The momentum limits for mesotrons were 
calculated, corresponding to all of the absorbers 
used in the experiment, from the range-mo- 
mentum relations given by Rossi and Greisen.” 
(In determining the ranges, the average absorber 
thicknesses were used, taking into account the 
differences in absorber thickness for mesotrons 
arriving in different directions.) The counting 
rates obtained with the iron and lead absorbers 
were plotted as a function of these momentum 
limits, yielding the two upper curves in Fig. 5. 
At each abscissa on this graph, the ordinate 
representing the mesotron intensity must have 
been the same for both the iron and the lead 
absorber. The difference between the curves 
arises from the fact that more electrons penetrate 
the iron; i.e., thicknesses of iron and lead which 
are equivalent for mesotrons are not equivalent 
for electrons. In fact, for electrons of energy large 
compared with the critical energies (25 Mev in 
iron, 7 Mev in lead), the iron and lead are equiva- 
lent for electrons in radiation lengths (14.4 g/cm? 
in iron, 5.9 g/cm? in lead). Thus in order to pene- 
trate the maximum thickness of iron, the elec- 
trons must traverse, on the average, in addition 
to the carbon, only 4.75 radiation lengths of iron, 
while the corresponding thickness of lead is 15.8 
radiation lengths. 

It was assumed that no electrons penetrate the 

thickness of lead equivalent (for mesotrons) to 


1 P, Auger, Phys. Rev. 61, 684 (1942). 
a be B Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
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the maximum thickness of iron." The difference 
between the two absorption curves at this 
abscissa represents the number of electrons which 
can penetrate a smaller thickness of lead (the 
thickness which is equivalent for electrons to the 
maximum thickness of iron). By subtracting this 
difference from the lead absorption curve at the 
appropriate abscissa, one point is obtained on an 
extrapolated curve representing the total meso- 
tron intensity (see Fig. 5). The difference between 


TaBLe II. 

Av Total num- Electrons Electrons 
energy limit ber of from from 

for elec- electrons _ collision _ remaining 

Absorber trons (Mev) __ per min. processes sources 
None 3.4 20.8 7.6 13.2 
1.36 g/cm?C 6.4 17.6 5.6 12.0 
2.84 g/cm?C 10.0 15.1 4.2 10.9 
4.2 g/cm 13.6 12.8 3.4 9.4 
6.1 g/cm 18.9 10.8 2.7 8.1 
2.84 g/cm?C 

+62.8 g/cm*Fe 350 2.0 0.0 2.0 


this extrapolated curve and the iron absorption 
curve at another thickness of iron may be sub- 
tracted from the lead absorption curve at a still 
smaller abscissa. By proceeding in this way, the 
six points on the ‘‘Mesotron”’ curve in Fig. 5 have 


been obtained. We are assisted in drawing the — 


curve by the fact that it must be tangent to the 
lead absorption curve at the greatest thicknesses, 
and must become horizontal as the momentum 
approaches zero (not only because of the increase 
in the decay probability per cm as the velocity 
decreases, but also because of the increase in the 
energy loss). 

The only questionable part of the above pro- 
cedure is the determination of thicknesses of iron 
and lead which are exactly equivalent for elec- 
trons. Actually, different absorbers can never be 
equivalent for all electrons, because of the differ- 
ence in critical energies. However, the carbon 
absorber kept permanently under the iron and 
lead removes from consideration the mary elec- 
trons of energy below 10’ ev which emerge from 
the lead. If the energy limit set by the carbon 
absorber had been somewhat higher (several 
times the critical energy in iron, 25 Mev), we 
should have been able to say the iron and lead 


13 This assumption is justified by the graphs in Fig. 4b, 
which show that at this thickness the transition effect for 
all particles traversing the lead agrees with the transition 
effect for mesotrons alone. 
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absorbers were exactly equivalent in radiation 
lengths; as it is, this is still true after such thick- 
nesses of iron and lead that the energy limit has 
been so increased. Thus, the average energy loss 
of the electrons barely capable of penetrating the 
absorber, in the first radiation length of iron, is 
the same as that in the first 1.6 radiation lengths 
of lead; beyond these thicknesses, the iron and 
lead are equivalent in radiation lengths. Hence, 
for values of x greater than 1, we have assumed x 
radiation lengths of iron equivalent to x+0.6 
radiation lengths of lead. 

From the ‘‘Mesotron” curve in Fig. 5 we have 
obtained the counting rate due to mesotrons 
corresponding to each of the absorbers used in 
this experiment; in particular for the carbon ab- 
sorbers, for which the average energy limits set 
for electrons could be calculated. By subtracting 
the counting rates due to mesotrons from the 
total counting rates, we have obtained the count- 
ing rates due to electrons, which are given in 
Table II. The number of electrons penetrating 
the maximum iron absorber has also been given, 
for which the approximate energy limit has been 
calculated with the Snyder-Serber theory." 


MESOTRONS DESCENDED FROM A 
HIGHER ALTITUDE 


Under the assumption that practically no 
electrons penetrate the maximum lead absorber," 
the counting rate with this absorber gives the 
intensity of ‘‘fast’’ mesotrons (i.e., mesotrons of 
momentum above 2.5X10® ev/c, the average 
limit set by this absorber). If the number of 
mesotrons at Ithaca with momentum below this 
limit can be calculated, we have an alternative 
method of obtaining the total number of meso- 
trons and, hence (by subtraction from the ab- 
sorption curves), the number of electrons. Pro- 
vided that no mesotrons are produced between 
the two elevations, we may calculate the number 
of low energy mesotrons at Ithaca as the number 
which have descended from a higher elevation, 
where they had a greater energy. 

Mesotrons of momentum between 0 and 
2.5108 ev/c at Ithaca (275-meters elevation) 
would correspond to mesotrons with momentum 


4H. Snyder, Phys. Rev. 53, 960 (1938) ; R. Serber, Phys. 
Rev. 54, 317 (1938). 
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between 6.6 and 8.1X10*8 ev/c at Echo Lake 
(3240-meters elevation), where a differential mo- 
mentum spectrum of mesotrons, including this 
range, has been obtained by Rossi, Greisen, 
Stearns, Froman, and Koontz.'® By using this 
momentum spectrum, together with the proba- 
bility of survival of the mesotrons between the 
two elevations,'* we have calculated the number 
of slow mesotrons which have survived to reach 
Ithaca from the elevation of Echo Lake. By 
adding these to the number of fast mesotrons, we 
have obtained the dashed curves shown in Fig. 5, 
which should represent the total mesotron in- 
tensity. The two curves correspond to calcula- 
tions made with two different values of the 
lifetime-to-mass ratio for mesotrons, 8.4 10~‘ 
cm-c/ev (upper curve) and 6.7X10-* cm-c/ev 
(lower curve), between which the correct value 


probably lies. The value 8.4X10~* cm-c/ev is 


taken from the measurements of Rossi and 
collaborators,'® but may be somewhat too high if 
mesotron production is an important process at 
low elevations. A recent direct measurement by 
Rossi and Nereson"’ indicates that the mesotron 
lifetime is about 2.2 microseconds; this together 
with a mesotron mass of 160 electron masses 
would indicate 8.4X10-* cm-c/ev for the life- 
time-to-mass ratio, in agreement with the other 
measurements, but with a mesotron mass of 200 
electron masses indicates 6.7 X10-* cm-c/ev for 
the lifetime-to-mass ratio. 

The graphs in Fig. 5 show that this calculation 
of the number of slow mesotrons is incorrect. In 
the first place, it yields a graph of total mesotron 
intensity which is not tangent to the absorption 
curve in lead at the largest thicknesses. Secondly, 
it predicts that the number of slow mesotrons is 
only 3 percent of the total number of mesotrons, 
while the differences between the absorption 
curves in iron and lead indicate that the number 
of slow mesotrons is 8 percent of the total num- 
ber. From this discrepancy, we may conclude 
that more than half of the slow mesotrons qb- 
served at Ithaca have originated in the air below 
3240-meters elevation. 

It should be observed that this conclusion has 


16 B. Rossi, K. Greisen, J. C. Stearns, D. K. Froman, and 
P. G. Koontz, Phys. Rev. 61, bn (1942). 

16 See B. Rossi, Rev. Mod. 11, 

17 B. Rossi and N. Nereson, Phys. Rev. 62 (1942). 
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little bearing on the place of origin of the fast 
mesotrons. If one takes into account not only the 
decay but also the decrease in rate of mesotron 
production with increasing atmospheric depth, a 
rough calculation indicates that most of the high 
energy mesotrons observed near sea level should 
have originated near the top of the atmosphere, 
while most of the low energy mesotrons should 
have originated near the place of observation. 


ABSOLUTE INTENSITIES 


From the counting rates given above, it is 
possible to obtain, by a purely geometrical calcu- 
lation, the absolute integrated intensities of 
mesotrons and electrons at Ithaca. The inte- 
grated intensity to which the counting rates are 
most closely related is the one which is connected 
with the specific ionization, and which gives the 


number of particles per minute crossing a sphere. 


of unit cross section. In terms of the absolute 
directional intensity 7(@) and the zenith angle @, 
this integrated intensity is defined as 


f sin 61(6)dé. 
0 


Because of the cylindrical shape of the counters, 
the rates obtained in this experiment do not 
yield J directly, but yield an elliptic integral of 
the directional intensity, which may be written as 


R=2s gin 01(0) F(6)d8. 
0 


F(@) is a slowly varying factor, the total variation 
in which is less than 30 percent of its value at 
6=0. Hence we may write R= KJ, where K isa 
constant which depends slightly on the zenith- 
angle variation of the cosmic-ray intensity, and 
may be found from a tedious calculation in- 


TABLE III. Absolute integrated intensities of mesotrons 
and electrons. 


Particles per 

Component of cosmic rays cm? per min. 
Mesotrons above 2.5 X 108 ev/c 1.082 +0.007 
Mesotrons above 0.5 X 108 ev/c 1.175+0.012 
Electrons above 3.4 Mev 0.339+0.015 
Electrons above 6.4 Mev 0.287+0.014 
Electrons above 10.0 Mev 0.246+0.012 
Electrons above 13.6 Mev 0.208+0.014 
Electrons above 18.9 Mev 0.176+0.015 
‘Electrons above 350 Mev 0.033 +0.008 
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volving the dimensions of the apparatus. In the 
calculation of K, account is taken of the fact that 
the outer counters in the present experiment did 
not completely cover the solid angle above the 
central counter; however, only a few percent of 
the particles were missed in this way, because of 
the small number arriving at large zenith angles, 
Because the electron intensity has a somewhat 
steeper zenith-angle dependence than has the 
mesotron intensity, we have obtained slightly 
different values of K for electrons and for 
mesotrons (K=59.6 for mesotrons; 61.4 for 
electrons). The absolute integrated intensities 
thus calculated are listed in Table III. The errors 
given include not only the statistical errors but 
also an estimate of the likely errors which may 
have arisen in the method of computation. 


ELECTRONS ARISING FROM COLLISION 
PROCESSES OF MESOTRONS 


The number of electrons of energy above 10’ ey 
in equilibrium with the mesotrons at sea level, as 
a result of the collision processes, has been calcu- 
lated accurately by Rossi and Klapman,’ who 
obtained the result that the number of electrons 


was 6.7 percent of the number of mesotrons with © 


momentum above 3X10* ev/c. This would be 
about 6.6 percent of the mesotrons with mo- 
mentum above 2.5 X 108 ev/c, the limit set by the 
maximum absorber in the present experiment. 
Thus the counting rate for collision electrons 
above 107 ev (which corresponds to the energy 
limit set by one of our carbon absorbers) should 
be 4.2 per minute. 

An approximate method of calculation of the 
number of collision electrons with energy above 
other limits has been suggested to the author by 
Professor Rossi and rests on the conservation of 
energy. Since the collision electrons are in equilib- 
rium with the mesotrons, the energy given to 
electrons of energy above E in each g/cm* must 
be equal to the energy dissipated by them in each 
g/cm*. This energy is lost in two ways: (1) by 
electrons drifting through the boundary E to 
lower energies, each taking away an energy E£, 
and (2) by collision and radiation losses of elec- 
trons which remain above the boundary E. The 
last term is very closely approximated by the 
number of electrons of energy above E times the 
energy loss of electrons of energy E. The number 
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Fic. 6. Electrons arising from collision processes of mesotrons in air. 
N.(E) =number of electrons of energy above E in equilibrium with the 
. mesotrons, per 100 mesotrons of momentum above 2.5 X 108 ev/c. 


which drift through the boundary to lower ener- 
gies is equal (under the same approximations) to 
the number of new electrons added to the group 
by mesotron collisions. This term, as well as the 
energy given to the electrons of energy above E, 
may be calculated from the collision probabilities 
of mesotrons and integrated over the momentum 
spectrum of mesotrons given by Blackett.'* The 
results of such calculations are presented as a 
function of E in Fig. 6. 

The accuracy of these calculations may be 
tested by comparing the result for E=10" ev 
with the value obtained for this energy by the 
more accurate calculations of Rossi and Klapman. 
The present calculations give N.(10")=6.4 per 
100 fast mesotrons, while those of Rossi and 
Klapman give WN,(107)=6.6, in very close 
agreement. 

The counting rates due to collision electrons, 
corresponding to the carbon absorbers used in 
this experiment, are listed in Table II.’ By 
subtracting the collision electrons from the total 
electron intensity, we have obtained the electron 


18 P. M. S. Blackett, Proc. Roy. Soc. 159, 1 (1937). 

'* These have been obtained from the graph in Fig. 6 by 
multiplying the number of electrons per fast mesotron by 
the counting rate due to fast mesotrons in the present 


he 
“Oran rege The results have been corrected in the ratio 
6.6/6.4 to give exact agreement at E = 10? with the calcula- 


tions of Rossi and Klapman, which were considered more 
accurate than the present calculations. 


intensities listed in the last column of that table. 
The principal source of these remaining electrons 
must be the decay of mesotrons ; but regardless of 
whether or not some of the electrons have other 
origins (e.g., primary electrons or radiation by 
mesotrons), the numbers listed in this column 
must be at least as great as the numbers of elec- 
trons arising from mesotron decay. 


ELECTRONS ARISING FROM DECAY OF 
SLOW MESOTRONS 


The primary electrons arising from decay of 
mesotrons with momentum below 2.5 ev/c 
have energies of the same order as the critical 
energy in air; hence the shower theory cannot be 
applied to calculate the number of electrons to 
which they give rise. Therefore we have used 
instead the approximate method of calculation 
which was used above for the collision electrons. 
Because the energy spectrum of the decay elec- 
trons is not as steep as that of the collision 
electrons, the approximations are not as good, 
but because of the small number of slow meso- 
trons, a slight error in this calculation is of no 
significance. For the momentum spectrum of slow 
mesotrons, over which an integration is required, 
we have used the results deduced in the present 
paper from the differences between the absorp- 
tion curves in iron and lead. The results of the 
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TABLE IV. Electrons arising from mesotron decay. 


Energy limit Assumed decay electron, : 1 electron, 
(Mev) Ss g F T T T 
3.4 2.5 19.5 22.0 2.0 15.6 17.6 14.7 11.7 
64 22 17.9 20.1 1.8 14.3 16.1 13.4 10.7 
10.0 2:0 16.4 18.4 16 13.1 14.7 12.3 9.8 
13.6 18 15.3 17.1 1.45 12.2 13.7 11.4 91 
18.9 1.6 14.0 15.6 L. 11.2 12:5 10.4 8.3 
350 0.0 2:3 23 0. 18 18 1.5 12 


Legend: S = 


from decay of fast mesotrons; T =S +F =total predicted counting rate due to e' 


calculations are listed in Table IV, corresponding 
to the two values of the lifetime-to-mass ratio 
with which the calculations have been made (see 
discussion above in the section on mesotrons 
descended from a higher altitude). These calcula- 
tions have been made under the assumption that 
the decay products of the mesotron are an 
electron and a neutrino. If the decay products are 
an electron and a photon, the numbers should be 
multiplied by approximately 2; while if the 
mesotron decays into an electron and two 
neutrinos, the results should be multiplied by 3. 


ELECTRONS ARISING FROM DECAY OF 
FAST MESOTRONS 


Since most of the primary electrons arising 
from the decay of fast mesotrons have energies 
large compared with the critical energy in air, the 
shower theory can be applied to calculate the 
number of electrons to which they give rise. It 
has been shown by Rossi and Greisen” that if one 
neglects the decrease in the rate of production of 
the primary electrons over a distance equal to the 
maximum range of the showers, then the number 
of electrons of energy above E in equilibrium 
with the mesotrons is given by the integral track 
length of all electrons above energy E produced 


' per g/cm? by the mesotrons. In another article,® 


the same authors have shown that the variation 
in the rate of production of the primary decay 
electrons can be taken into account by evaluating 
the track length for the electrons produced per 
g/cm? at a distance of 130 g/cm? above the place 
of observation. 

The track length for E=0 is obtained simply 
from the conservation of energy as 43 Eo/e, where 
E> is the energy going to decay electrons per 
g/cm’, and ¢€/43 is the ionization loss per g/cm’, ¢ 


icted counting rate due to electrons arising from decay of slow mesotrons; F = predicted counting rate due to electrons arising 
arising from decay. 


mesotron 


being the critical energy, and 43 g/cm? being the 
radiation length in air. For other energies, we 
have used the accurately known track length for 
E=0 together with the relative track lengths (or 
energy distribution at the maximum of a shower) 
recently calculated by Richards and Nordheim.* 
The track length for E=10" ev, calculated by 
Rossi and Klapman,’ agrees very closely with the 
corresponding value obtained by Richards and 
Nordheim, the latter value being higher by about 
3 percent. 

The energy Eo going to decay electrons per 
g/cm?, due to the disintegration of fast mesotrons, 
is given by Ey=43(u/r0)(N/p) (where p is the 
density of air, N is the number of fast mesotrons, 
and 70/y is the lifetime-to-mass ratio), under the 
assumption that the disintegration products are 
an electron and a neutrino. If the mesotron 
decays into an electron and a photon, E» should 
be twice as large; if into an electron and two 
neutrinos, Ey should be ? as large. Ey must be 
evaluated at 130 g/cm? above the place of ob- 
servation; hence the values of N at 130 g/cm’ 
above Ithaca were read from a graph of mesotron 
intensity vs. elevation recently obtained by the 
author.?° Data were available at four different 
zenith angles, so that an integration over zenith 
angle could be performed, as was required for 
comparison with the data of this experiment. 
Thus we have obtained the numbers in Table IV, 
which represent the hypothetical counting rates 
due to electrons arising from decay of fast 
mesotrons. 


CONCLUSIONS 


The total numbers of electrons arising from 
decay, given in Table IV, should be comparable 


2° K. Greisen, Phys. Rev. 61, 212 (1942). 
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with the counting rates listed in the last column 
of Table II; at any rate, they should be no larger 
than those counting rates. Hence it is impossible 
that the mesotron should decay into a photon 
and an electron, for in that case the counting rate 
due to electrons arising from decay would have to 
be three times what was observed; in fact, it 
would have to be 50 percent greater than the 
total difference in counting rate between the 
minimum absorber (only the counter walls) and 
the maximum absorber (including 101 g/cm? of 
lead). 

Because of the indirect way in which the 
counting rates listed in the last column of 
Table II were obtained, they cannot be con- 
sidered very accurate. However, considering all 
the possible sources of error, it is unlikely that 
any of these numbers should be in error by as 
much as 3 counts per minute, or that the last of 
the numbers should be in error by as much as 0.5 
counts per minute. Hence we may conclude that 
the results also discourage the hypothesis that 
the decay products are an electron and a neutrino. 
Atany rate, if these should be the decay products, 
the results indicate that the lifetime-to-mass 
ratio for mesotrons is at least as large as 8.4 10-* 
cm-c/ev. Thus, if one accepts the value of the 
lifetime (2.2 microseconds) recently obtained by 


Rossi and Nereson,'’ the mass of the mesotron 
must be no larger than 160 electron masses. 

In general, however, the results favor the 
hypothesis that Jess than 4 of the energy of 
decaying mesotrons goes to the electron com- 
ponent. They would be in agreement with a 
decay into an electron and two neutrinos, for 
instance, in which case the spin of the mesotron 
might be } but not 0 or 1, and only 3 of the 
mesotron energy would go into shower pro- 
duction. 

It has been suggested!*"” that some of the 
mesotrons may disappear by another process 
other than radioactive decay and, hence, without 
contributing to the electron component. If this 
were an important process in the atmosphere, it 
might reconcile the disintegration into an electron 
and a neutrino with the small number of electrons 
observed. However, because of the small density 
of the atmosphere, and because most of the 
electrons should arise from disintegration of 
mesotrons of great momentum, we do not think 
that the suggested process can affect the con- 
clusions discussed above. 

The author wishes to express his gratitude to 
Professor Bruno Rossi for his invaluable aid in 
the design and analysis of this experiment. 


21S, pomenaee and G. Araki, Phys. Rev. 58, 90 (1940). 
=F. Rasetti, Phys. Rev. 60, 198 (1941). 
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Meson Theory of the Magnetic Moment of Proton and Neutron 
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The magnetic moment of proton and neutron is determined by a second-order perturbation 
calculation. A new limiting process secures the convergence of the integration in momentum 
space. The result for both pseudoscalar and vector meson is the same and leads to a positive 
magnetic moment for the neutron if the magnetic moment of the vector meson is assumed to be 
one meson magneton (eh/2uc). Agreement with the known values of the magnetic moment of 
the heavy particles can only be obtained if we give the vector meson an anomalous magnetic 


moment of about three meson magnetons. 


1. 


HE relativistic wave equation of particles 
with spin 3 leads in a natural way toa value 
of the magnetic moment of one magneton 
(eh/2mc), where m is the mass and e the charge of 
the particle. In the case of proton and neutron the 
magnetic moments differ appreciably from the 
values which they should have according to this 
theory. Measured in units of nuclear magnetons 
the experimental values for these moments are: 


= +2.785+0.02 


for the proton! and 
—1.935+0.02 


for the neutron.? The minus sign in the case of the 
neutron indicates that the spin points in the 
opposite direction of the magnetic moment. 

It was pointed out by Kemmer’ and others‘ 
that the meson theory of nuclear forces would 
also lead to anomalous magnetic moments of the 
nucleons, due to the virtual emission and re- 
absorption of mesons by the nuclear particles. 
Unfortunately a second-order perturbation calcu- 
lation leads to divergent integrals. This diver- 
gence difficulty was dodged by cutting off the 
integrals suitably, which is equivalent to intro- 
ducing a finite size of the nucleon. 

This procedure, however, is highly unsatis-. 
factory for two reasons: 


728 (1999). : , Rabi, Ramsey, and Zacharias, Phys. Rev. 56, 
2 L. W. Alvarez and F. Bloch, Phys. Rev. 57, 111 (1940). 

3 N. Kemmer, Proc. Roy. Soc. A166, 127 (1938). 

4 Fréhlich, Heitler, and Kemmer, Proc. Roy. Soc. A166, 
154 (1938). H. Fréhlich, Phys. Rev. 62, 1801. (1942). The 
result announced in this note, that the pseudoscalar meson 

ves no magnetic moment, is incorrect. I am indebted to 

ofessor Fréhlich for correspondence on this point. 


(1) It introduces an arbitrary constant into 
the theory, the physical significance of which is 
not quite clear. 

(2) It destroys the relativistic invariance of 
the theory. 

In spite of these grave objections one feels that 
the meson theory of the magnetic moment is 
essentially correct and that a theory of interaction 
of nucleons with the meson field, formulated in 
such a way that the divergence difficulties disap- 
pear, should give the correct value for these 
moments. 

Recently a theory was proposed by Dirac tal 
the interaction of an electron with the electro- 
magnetic field, which avoids all divergence diffi- 
culties for that case. The theory consists of two 
essentially different parts. The first part is the so- 
called \-limiting process, first investigated by 
Wentzel.® It is essentially a classical procedure 
and consists of using the difference of retarded 
and advanced potential at the position of the 
particle. Dirac showed’ that the \ process may 
best be introduced in quantum mechanics by a 
change in the commutators of the field quantities 
entering into the interaction energy.® In this form 
the procedure may easily be generalized to 
charged fields with rest mass different from zero. 

The integral which causes the divergence of the 
meson theory of the magnetic moment is the 
integral ° 

k*/(k?+?)*dk 


P. A. M. Dirac, Proc. Roy. Soc. A180, 1 (1942). 
® 4 Wentzel, Zeits. f. Physik 86, 479, 635 (1933) ; 87, 726 
P.A 


. M. Dirac, Am. de I'Inst. H. Poincare, 9, 13 (1939). 
8 For a detailed account see an article by W. Pauli, to be 
published in the Rev. Mod. Phys. 
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with »=meson mass. This integral will be made 
convergent with the A process. We must, how- 
ever, remark that this integral occurs only if a 
second-order perturbation theory is applied and 
if the mass m of the heavy particle is infinitely 
large. In higher approximations in the coupling 
constant and in »/m other integrals will appear in 
addition to this which no longer converge with 
the A process alone. In general only integrals 
which for large values of k behave like 


f . 


with m integer will be made to converge with this 
method. A limitation to a second-order pertur- 
bation treatment would not be justified if such 
higher order divergent terms remain. 

In order to make the other integrals converge, 
Dirac has introduced another kind of field oscil- 
lator corresponding to states with negative 
energy and quantized according to an indefinite 
metric in Hilbert space.* Dirac has shown that 
the two methods lead to the elimination of all 
divergencies in the quantum theory of fields. 
Serious interpretation difficulties still remain, 
however. 

It may easily be seen that the higher order 
terms lead always to integrals of the form® 


f "R(ko, dk 


with a rational function R(ko, k) for the integrand 
containing only even powers of k. The introduc- 
tion of the negative energy oscillators changes this 
integral into the form 


af “[R(ko, b)-+R(—ko, 


and this is sufficient to make all the remaining 
integrals occurring in this perturbation treatment 
converge. 

The calculation of the magnetic moment due 
to the meson field with this method leads to a 
positive value for the neutron both in the 
pseudoscalar and the vector theory. This is in 
disagreement with the experimental result. Only 
by making use of the possibility that the vector 


* This is not true for the theory of holes, where integrals 
the form /dk/(k*+ occur, which are logarithmically 
divergent, even in this theory. 
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meson may have a magnetic moment different 
from one meson magneton can agreement with 
the experiment be obtained. The contribution of 
the meson field to the magnetic moment has then 
a certain value M which is positive for the proton 
and negative for the neutron, so that we shall 
have for the magnetic moment of these particles 


Mp=1+M, 


By choosing M~+1.9, agreement with ex- 
periment will be obtained within the limits of our 
approximation. 

The objection may be raised that this assump- 
tion of an anomalous magnetic moment of the 
vector meson introduces a new constant into the 
theory and we stand where we were before. The 
situation is, however, not quite so bad. We must 
bear in mind that we have introduced the new 
parameter in a way which does not destroy the 
relativistic invariance of the theory. Moreover, 
the constant has a definite physical meaning, 
namely, the magnetic moment of the: vector 
meson, which may be checked with an experi- 
ment, deciding for or against the theory here 
presented. 


2. THE PSEUDOSCALAR THEORY 


In the following we often use one single letter 
“a’’ for a four vector with components a,. Here 
=a* is the space part of the vector. The invari- 
ant scalar product of two four vectors shall be 
denoted by (a, b) = (a, b) —a°b°=a,a*. (a9= —a’.) 
Charged particles of spin zero are described by 

a complex field ¥(x) satisfying the field equation 


LW 
where uw is the meson mass in units A/c and 
[ ]=0?/dx*dx,. It is quantized according to the 
commutation rules 
iLy(x), ¥*(x") ]=D(x—x’), (1) 
where D(x) is the invariant D function defined by 


exp [7(k, x) ]sin k®°x®, (2) 


D(x, x®)= 


1 
+ (k?+y*) 
The operator for the current density véctor is 


- (3) 
Ox* 
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and the total magnetic moment of the field 
e 


where the positive number ‘‘e’”’ is the magnitude 
of the charge of the mesons in units (hc)}. 

For our calculations it is somewhat more con- 
venient to introduce two real fields y; and y2: 


1 1 
—tp2), Wits), 


for which we have the commutation rules 
iLWa(x), Va(x’) (x —2’). (5) 


The spin dependent interaction of the heavy 
particle with the coordinates z may be written" 


H' = grad)ya(2). (6) 


The ra represent the matrices for the isotopic 
spin and o=(01,¢2,¢3) is the ordinary spin 
vector. They satisfy the equation 


T1T2>= —T2T1=173, °° oa» 
(7) 


together with all the other equations which can 
be obtained from them by cyclic permutation of 
the indices. The summation over a is extended 
over the two indices a= 1, 2. In the symmetrical 
theory there is a third real field ys present repre- 
senting the neutral mesons, and the summation 
goes then over all three indices. Since neutral 
particles do not contribute anything to the cur- 
rent and to the magnetic moment, we can restrict 
ourselves to the first two indices. 

The interaction energy contains only the non- 
relativistic part of the heavy particles. This ap- 
proximation can be made since we are only 
interested in the effects which are of highest order 
in »/m, and they are independent of m. This 
amounts to setting the mass of the heavy particle 
infinitely large, and then relativistic effects of the 
heavy particle are of course unimportant. The 
mesons, however, are treated relativistically all 
through the calculations. The interaction term in 


10 In the complex notation this term corresponds to the 
expression 
= g(8x)'[r_(@, grad)y*+r+(0, grad)y] 


with ra =1/2(riir2). Our g is therefore v2 times the 
coupling constant which is usually used in the literature. 
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the Hamiltonian leads to an additional term in 
the current so that the total current now is" 


Su= (Opi /Ox")P2— 
— ri2)5(x—z).- (8) 


We develop the field operators into eigenstates 
in momentum space 


f { 2) 
(9) 


This definition of the ga(k) and g*(k) is useful 
since they are constant operators if there is no 
interaction between heavy particles and field. In 
case of an interaction they are slowly varying 
according to the law” 


ga}. (10) 

From (1), (2), and (9) it follows that the 
commutation rules for the ga(k) are 

Cea(k), s*(k’) (11) 

from which it follows in the well-known way that 
Na(k) = (1/k°) ¢a*(K) ¢a(k) 


has the eigenvalues 0, 1, 2, - - - and represents the 
number of mesons in the state (a, k). 
The interaction energy is 


d*k 
H'=ifre f k) 
X (12) 


f=g/2r. 


The \ process may be introduced, as Dirac” 
has shown, by a change in the commutation rules 
of the field operators which enter into the inter- 
action energy. In x space these commutation 
rules may be written 


iLVo(x), Vox") (13) 
Dy(x) (14) 
A= {2,A°} isa time-like four vector, (A, 4) <0, (15) 


Wa(x) 


with 


where 


"See W. Pauli and S. M. Dancoff, Phys. Rev. 62, 105 
(1942). There would be a third term in the total current due 
to the proton which, however, in our approximation is zero, 
since the velocity (although not the momentum) of the 
nucleons is zero. 

12 See W. Pauli, Rev. Mod. Phys. 13, 210 (1942). 

13 P. A. M. Dirac, Ann. de I'Inst. Henri Poincaré 9, 13 


(1939). 
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which ultimately goes to zero in such a way that 
the inequality (15) is always satisfied. If A 
transforms like a four vector, then the scheme is 
evidently relativistically invariant. In k space the 
introduction of the \ process shows itself in the 
commutation rules 


[ea(k), cos (16) 


which follow from (13). 

A possible representation of the G.(k) in terms 
of the ga(k) is given by 
A), 


—sin ——}, 


= a cos 


k,X 
Ga* (Ik) = ) 


The quantities which depend on the field alone 
remain unchanged." For instance, the magnetic 
moment may be. written 


(17) 


Ct ¢1* * l, x) } 
{ z) + I’, z) } 


—sym. (12)]—sym. (ws), (18) 


where sym. (uv) stands for the preceding expres- 
sion with the indices » and v interchanged. 


“v2 (20)! 


x f 


+ 9} —sym. (1, 2)]. (19) 
The first of these two expressions may be 
simplified by replacing x, by a differentiation of 


“This definition of the adjoint ators of the ¢ is 


a generalization of the Hermitian conjugate, similar to the 
one described by Pauli in a paper to be published in the 
Rev. Mod. Phys. 
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the exponential and subsequent integration by 
parts 


ue fe dl’ 
lo 


dgill) | d¢gi*(1) 


—sym. (12) ]—sym. (ut). (20) 


Using 
1 34 — 
we find 
) 


¢1(1) +compl. con}. ) 


—sym. (1, 2} —sym. (w) | (21) 


The Schrédinger functional Q satisfies the 
Schrédinger equation 


1(02/dx*) = H’Q. (22) 


We look for stationary solutions of this equation 
which are correct up to the second order in g. For 
that purpose we write 


where Q, 91, 22, -- + represent those parts of the 
functional with 0, 1, 2, --- mesons present. The 
leading terms of these functionals will then be of 
zeroth, first, second, --- order in the coupling 
constant g. 

Let u represent the state of the nucleon which 
we choose in such a way that o3u=u. The spin 
then points in the positive z direction. Further- 
more, let w represent that state of the field in 
which no mesons are present, so that 


Ga(k)w=0. (23) 


Then we have for Q)= tw. 
We then write the equations for the Q’s 


i(d2,/dx®) = H’Q, 


i(dQ2/dx®) = H’Q,. 
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If we can integrate these equations, we can use 

the ’s to calculate the expectation value of 

M. 3- 

2* M3Q= 22 M322 
(25) 

since M; is Hermitian. 

It is seen from the last equation that only 
terms of zero order in g which occur in Q will 
contribute anything to the magnetic moment in 
our approximation. With this remark in mind we 
can integrate the equations (24) approximately. 
The first one with the help of (22) and (23) reads 


d d*k k i(k, 2) 
~ifre (0, k) ga*(k)e*™ 
or 


and from this and the second equation of (22) we 
get 
ap! 


d 
i—%=f*rats | 


dx® Re’ 


(bot by!) ba? k’)(o, k) 
X Ga*(k’) +* uw. 


These integrations are only approximately 
correct since the operators g, and ¢,* are slowly 
varying according to (10). But this variation 
would lead to a correction term to Q of the 
second order and a similar term in Q, of the third 
order. Another term of the second order has been 


neglected by including in the expression H’Q; 


only those terms which correspond to two mesons 
present. In view of the above remark, we may, 
for our approximation, neglect all these terms. 
With these wave functions we can now calculate 
the expectation value of the magnetic moment. 
First we prove that the last term in (25) does 
not contribute anything in this approximation. 
This term is proportional to 


. 
Fo *)—sym. (12)} 


XTa k’) uw, 
0 
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which is equal to 
3 
—iu* (2102 — 2201) Ts f k)u cos (R, d) 


when we use the commutation relations. The 
integral vanishes and therefore 


M32, =0. 


Of the remaining terms let us first calculate 
2,*M;2,=M;'. We notice at once that only the 
part in M; which contains an emission and an 
absorption operator can give a contribution to 
this expression. The remaining part may be 
simplified by adding to it the expression 


0 0 
k 
( { g1* ¢192*} 


which vanishes when integrated over k space. In 
this way we obtain for the part M;! of M;'5 


d2*(k) 
Ok, 


dk d*k! 
Mit f ss 


d* 


a 
+ ¢(l) —sym. | 


X uw. 
This expression may be simplified by using the 
commutation rules (11). We have 


¢2(1) +e ()) 


cy) 


(ora 
~sym. (xy) 
= 
al, 


46 In order to avoid ambiguity we have put u=x, v=y. 
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so that after a partial integration'® 
My = —ifru*rs (0, k)u cos (k, d). 


If we choose a special coordinate system for 
which the space part of \ is zero, we may first 
carry out the integration over the angles and 
obtain 


The integral may be evaluated in the following 


f (=) cos koedk = f cos kododk 


Ro! 
The first integral is zero in the limit \>=0, as is 
shown in the appendix, and for the second we get 


—(3x/4)u, by going to the limit before 
carrying out the integration so that finally 


M3 = feu =F (gu)*(e/2u). (28) 


The upper sign is for the proton and the lower 
for the neutron. 

For the calculation of 22o*M;Q2 we need only 
to consider that part in M3; which contains two 
absorption operators 


cos kododk. 


sym. (12) 


—sym. 


X — sym. —sym. cy) | 


dk’ dk 
x f k’)(e, k) 


16 We have coment terms which arise from the differenti- 
ation of the cos-factor and which tend to zero as ) goes to 
zero and terms coming from the differentiation of the 
exponential function which leads to an expression of the 
form o1%:—02%; and which is zero after averaging over s. 
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Again applying the commutation rules we have 


bf 


(29) 


— 6(1+k’) 
|cos r), 


dk’ 
Ms" = 2iefru nf Ro(Ro-+ho’) Ro? 


X (oike’ — (o, k)5(kK+k’)u 
f (=) cos (k, A\)dku*r3u 
= (gu). (29) 
2u 


The total magnetic moment due to the 
pseudoscalar meson field is 


Ma= ¥2(gu)*(¢/2n) (30) 
for the proton and neutron, respectively. 
3. THE VECTOR THEORY 


The mesons of spin 1 are described by a com- 
plex vector field y, satisfying the equations 


dp,/dx,=0. (31) 
It is quantized according to 


2 

XD(x—x’). (32) 
The current is defined in terms of the quantities 
and Wyo = — (dp,/dx*). 


and the magnetic moment is then again given by 
Eq. (4) of the preceding section. 

We decompose the field again in the real fields 
by 


1 


1 1 


ate 
fix wf cos | 
to 
be 
way 
| 
| 
| 


We put ¢={y,}(k=1, 2, 3) and write for the 
interaction energy 


H’' = g(4r)*ra(e, curl) ¢(z), (34) 

with this interaction present the current is 

1 2 
Su= — 

(q'r2— },5(x—z). (35) 
The second part corresponds to the second 
part in (8) in the scalar theory. It will give rise to 
an additional term in the magnetic moment. A 
consideration similar to the one given above will 
show that this term does not contribute anything 


to the expectation value of the magnetic moment 
within the limits of our-approximation. 


We develop y,* and Viv in eigenstates in 
momentum space 


1 | d®k 
Vn" (x) = anid 

X + (36) 
and 

1 d®k 
= 


(Rugo®(k) — ky gy 
— (ky — (37) 


Because of the subsidiary condition dy,/dx"=0 
the ¢,*(k) satisfy the equations 


(¢%, k) k) =0. 


From the commutation relations (32) it follows 
that 


1 
=k° RyRy —k’ ap- 


The ¢,*, ¢,** again are defined in such a way 
that they are slowly varying according to 


=iLH’, (39) 


dx® 
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For H’ we get from (34) in momentum space 


d*k 


—(o, [kx (40) 
f=g/2r. 


The A process is again introduced by defining 
the operators g,* and ¢,** as was done in the 
scalar theory. As before this is equivalent to 
replacing the D function in (32) by the D, 
function. This alters the commutation relations 
(38) in momentum space to 


X 6(k—k’) cos (k, A). (41) 
The magnetic moment operator 
(42) 


may be transformed with a calculation exactly 
analogous to that which led to the expression (21) 
in the scalar theory to give 


X 
+(. )e 


+compl. con}. |—sym. —sym. (ut). (43) 


The solutions of the equations 
(44) 
corresponding to the solutions (26) in the scalar 
theory are then 
Qo = Uw, 
dk ~ 
(45) 


x [k’x ]) uw. 


—f*rats 


an 


| 
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| = 
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The part M;! corresponding to (27) is 


X(o, [kx 


+compl. 


Writing 
(o, [kX$¢°]) 


we can simplify the expression (46) considerably 
by using the commutation rules (37) 


d¢,*? 
l nf 1p 
th 


rok 
+compl. con}. |-sym. (29) 

= —k,k, al. 


1 
—sym. 


r# 
cos (kd). 
t#p 
After multiplying with 
(d*k’/ko’*) exp [t(k—k’, z) 
and integrating over k’ we get by using 
o,ok-k,=0 
r¥e 
—sym. (xy) cos (k, 
t#p 


so that (46) becomes 


d*k d*l 


ko 


ei T3 


—sym. (ay) *) cos (k,X). (47) 
r#p,t#p 


‘als 
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Integrating by parts over | first we get a term 
from the differentiation of 1/1)’ which, however, 
is zero just as in the scalar theory because it con- 
tains a factor since t¥p. 

The other part gives a term 


— 
+ (ks02— eos) 0103 
— +o2k2) — kiss) 0203. 


Most of these terms will give zero after the 
integration over the angles. The only term 
different from zero is 


introducing this in (47) we get 


8x k\4 
f (;) cos = F 2 
0 


The term M;"! is 


dk dk’ 
f 
3 freu*w* rats — he”? 


X sym. 


X (0, [kX [k’X$**]). 


A calculation similar to the one above leads to 


d®kd*] 


él, 


—sym. cos (k, d). 


In comparing this with (47) we find that again 
we have 


Ms =F 


so that the total magnetic moment in the vector 
theory is 


M3=M3'+ M3" = ¥2(gu)*(e/2u). 
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The upper sign is for the proton and the lower for 
the neutron. This is the same result as for the 
pseudoscalar meson. 


4. A VECTOR MESON WITH AN ARBITRARY 
MAGNETIC MOMENT 


The Lagrangian for a vector meson field is not 
unique insofar as it is possible to introduce terms 
in the Lagrange-function which are proportional 
to the external electromagnetic field.'7 Such a 
term makes itself felt in an altered expression for 
the current density and consequently a change in 
the magnetic moment. The additional current 
density is 


a 
(48) 
Ox, 
where K is a dimensionless and arbitrary con- 
stant. The change in the magnetic moment 
operator is then from Eq. (4) 


Myo =1eK f Vo", 
=eK 


=}eK fs 


+ ¢,'(K) ¢.**(k)) —compl. conj. ] 


—sym. (uv)}. (49) 
We have, therefore, 


d*k d*k’ 
[kx$*]) 


X[o, [k’X 


= —2eK r303— = f cos (RoAo) dk 


= +K(gu)*(e/2u). (50) 


( ne H. C. Corben and J. Schwinger, Phys. Rev. 58, 953 
1 4 


A similar calculation for = leads to 
the same result so that 


My! = = 42K (gu)*(e/2u). (51) 


The total magnetic moment in the vector 
theory is then 


Ms=+2(K—1)(gu)*(e/2u) (52) 


for the proton and neutron, respectively. 
Taking the two theories together and assuming 

for the vector meson twice the mass of the scalar 

meson, but the same coupling constant, we find 


M=+2(gu)?{2(K—1)—1}(e/2u) (53) 
for (gu)?~0.1, u~zom, K~2, 
M~ +2, 


for the magnetic moment of proton and neutron, 
respectively, due to the meson field only. If we 
add to that the magnetic moment of the nucleons 
(0 for the neutron and +1 for the proton), we 
find agreement with the experiment within the 
limits of our calculation. 

This problem was suggested to me by Professor 
W. Pauli of the Institute for Advanced Study. I 
am very grateful to him for his generosity in 
taking time for discussion and helpful advice. 


APPENDIX 


The integral = a cos kedodk may be evaluated in the 


following way. With the substitution k =i sin vAqu =a, we 
have 


| cos cos £08 ¥) dy, 
For the second integral we substitute v by r—v. 
cos ve cos ve cot 


The paths of integration have to be taken in the regions for 
which Im(cos v) >0, in order to secure the convergence of 
the integrals. They may be chosen in such a way that the 
contribution from the branches leading to infinity cancel, 
and we are then left with the expression 


cos —>-Ji(a) lim I(a) =0. 
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The mathematical problem of the perfectly conducting cylindrical cavity resonator containing 
two and three dielectrics is set up, and solved in terms of the electric and magnetic Hertzian 
vector potentials. The relations between the various natural frequencies and resonant lengths 
are determined and discussed, in particular for a cylindrical resonator containing a movable slab 
of dielectric. It is shown that modes of the form Expo and Hypo are not excitable if there is more 
than one dielectric present in the resonator. The band transmission characteristics of such 
resonators are considered briefly, and a simple and convenient absolute method for the measure- 
ment of dielectric constants at microwave frequencies is suggested. Experimental measurements 
were made at 3070 megacycles per sec., for the Hi:: mode, and the results were found to agree 
well with the theory. Experimental and theoretical curves illustrating various points of the 


discussion are included. 


I. INTRODUCTION 


T has been known for a long time that, if the 
frequency is sufficiently high, electromagnetic 
waves can be propagated inside conducting tubes 
of various cross sections, rectangular, circular, 
elliptical, etc.' Any closed, highly conducting 
chamber, practically not more than a few wave- 
lengths long, may serve as a cavity resonator, 
inside of which standing wave patterns may be 
set up. Such resonators may oscillate in a large 
number of different modes,?~‘ the mode of oscil- 
lation depending on the manner of excitation.® 
By varying the dimensions of the cavity, one can 
tune the resonator to different frequencies; die- 
lectrics (other than air), also, have a considerable 
effect on the resonant dimensions. In earlier 
work, King® and Lamont’ have discussed this 
property of a dielectric of shifting resonant posi- 
tions, with particular emphasis on the measure- 
ment of dielectric constants at ultra-high fre- 
quencies. King’s paper describes a parallel wire 
method for measuring the dielectric constant of 
fluids and solids at wave-lengths around 2 meters, 
based on the solution of the telegraphist’s equa- 
tions for the parallel wire system. The measure- 
ment depends on a simple mathematical expres- 


1 Lord Rayleigh, Phil. 125-132 (1897). 

*W. L. Barrow, Proc. I.R.E. 24, 1298-1328 (1936). 

*J. R. Carson and S. P. Schelkunoff, Bell Sys. Tech. J. 
15, 310-333 (1936). 

4S. P. Schelkunoff, Proc. I.R.E. 25, 1457-1492 (1937). 

5G. C. Southworth, Proc. I.R.E. 25, 807-832 (1937). 

*R. King, Rev. Sci. Inst. 8, 201-209 (1937). 

TH.R. L Lamont, Phil. Mag. 29, 521-540 (1940). 


sion involving the dielectric constant, the fre- 
quency at which the measurement is made, the 
dielectric thickness for both liquids and solids, 
and a bridge shift measured along the parallel 
wires in air. Lamont extends the theory to include 
the hollow pipe or waveguide.* Both works em- 
phasize the importance of the change in resonant 
length of the particular system used due to the 
presence of a dielectric medium, particularly for 
the measurement of the dielectric constant of the 
sample under test. 

This paper, however, is concerned primarily 
with the problem of determining the relations 
between the natural modes and resonant lengths 
of a cylindrical cavity resonator containing 
several dielectrics. Some of the properties of such 
resonators as tunable circuit elements, and as 
wave filters, are discussed and finally, as a special 
case of the general theory, a simple method of 
measuring dielectric constants is suggested. This 
method differs from those described by King* and 
Lamont’ in that a single resonance position is 
determined with the sample under test at one end 
of the resonator; the advantages are that the 
calculations are somewhat simplified and less 


8 A distinction should be made between a cavity resonator 
and a waveguide. The waveguide is considerably longer, in 
fact, is treated mathematically as having infinite length; 
infinitely long waveguides are equivalent in the practical 
case to a finite guide terminated in its characteristic 
impedance. The wa ide may in general be treated 
analytically as a two-dimensional problem, whereas for 
cavity resonators the entire configuration, including the 
ends, must be considered. 
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experimental data are required, and in the case of 
fluids the data may be more conveniently taken, 
as no container is necessary. The disadvantage, 
however, is that the method is less sensitive be- 
cause the maximum shift conditions are not 
necessarily satisfied. 

It is the purpose, then, of this work to examine 
some of the resonant properties of cylindrical 
cavity resonators containing two and three die- 
lectric media. The results indicate a new and 
convenient absolute method for the measurement 
of dielectric constants of good dielectrics, for 
wave-lengths around 10 cm. They show also that 
the resonator behaves like a band transmission 
filter for the particular mode excited, and further 
the analysis yields formulae involving the reso- 
nant lengths and frequencies, which, it is hoped, 
may be of use in future experimental measure- 
ments and in the design of tunable ultra-high 
frequency resonators. The theory applies gener- 
ally for all modes and resonant dimensions of a 
tube containing three different dielectrics. Spe- 
cifically, we will discuss the Eo; and Hj,;; modes 
and corresponding resonant lengths of 


(a) a cylindrical resonator of fixed over-all 
length containing a movable slab of 
dielectric, 

(b) a cylindrical resonator of variable over-all 
length containing a movable dielectric slab, 
and 

(c) a resonator of variable total length con- 
taining a slab of dielectric at one end. 


Solutions for the electric and magnetic fields, 
E and H, will be obtained in terms of Hertzian 
vector potentials,® and application of the various 
boundary conditions will then give us expressions 
involving the resonant lengths and frequencies. 


Il. MATHEMATICAL ANALYSIS 


If we assume that the media considered are 
homogeneous and isotropic and if the disturbance 
is periodic with angular frequency w= 2zf, we may 
express the wave equations in the following form: 


=0, s=1,2, (1) 


where M® and M® denote, respectively, the 
electric and magnetic Hertzian vector potentials, 


*J. A. Stratton, Electromagnetic Theory (McGraw-Hill, 
New York, 1941), p 28. 
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Fic. 1. Cylinder containing three dielectric media with 
constants eg, ua, d= 1, 2, 3. 


which are related to the field vectors by the 
formulas: 


H=vv +9 XT, (3) 


and as usual, ¢, uw, and o are the dielectric con- 
stant, permeability, and conductivity, and ¢’ is 
given by 

=€—ia/w. (4) 


Equations (1) may be converted into scalar wave 
equations,’ and the solutions may be easily ob- 
tained. For cylindrical coordinates it can be 
shown that for our problem a complete set of 
solutions of Eq. (6) is made up of 


n=0,1,2--- 
s) [COS 
) cos m=0,1,2--- (5) 
sin 
p=1,2---, 


in various combinations to satisfy the boundary 
conditions. Choosing potentials having a z com- 
ponent only, we have for the EZ, or transverse 
magnetic modes alone (s=1): 


z—1 1 Vn 
am, am,” 
Apia 
r 00 or 


and for the H, or transverse electric modes 


10 See reference 9, page 349 et seq. 
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alone, we have, (s=2): 


i 
H,=i (9) 


The quantities i;, ie, is are unit vectors in the 
r, 6 and z directions, respectively; see Fig. 1 
above; the values of d=1, 2, 3 designate the 
various media. Henceforth we will omit the sub- 
script z on the field vectors, etc. 

We shall assume that the cylinder is perfectly 
conducting and that the dielectrics are perfect. 
Then «= for the cylinder, and o=0 for the 
dielectrics. 

For both the electric and magnetic modes, the 
boundary conditions" require that the B( =H) 
and D(=eE) vectors vanish everywhere outside 
the inner boundaries of the cylinder since it is a 
perfect conductor, and that the normal compo- 
nents of theB field and the tangential components 
of the E field also vanish there, as well as at the 
interfaces of the dielectric. The surface charge 
density » and the surface current density I are 
determined from the remaining boundary condi- 
tions: 

n-D= , (10) 


n X H=I, (11) 


where n is a unit vector normal to the surface 
under consideration. Since the dielectrics are 
assumed to be perfect, homogeneous, and iso- 
tropic, the surface current density I is zero at the 
dielectric interfaces, and the charge density 7 
must be also. 

In obtaining solutions of the scalar wave equa- 


tions we will need one arbitrary constant Pod 
for the first region, 0O=z=/,, two constants, Fe 
and for the second region, 


and finally one more constant _ for the third 
medium, /;+/:,=zS/. In the separation of vari- 
ables the following relations may be shown to 


hold: 
= (wun) (You) (12) 


and it can also be shown that ow p is the same for 


1 See reference 9, pages 34-38, for a general discussion of 
boundary conditions. 


each dielectric, i.e., the system as a whole has 
but one resonant frequency for the particular 
mode excited. 

The complete solutions then are, for medium 


one: 
a cos m 
sin Ym12 


Xexp (13) 
for medium two: 
(s) COS Ym2% sin 


II, => (45. n2 
SIN Ym22 COS Ym22 


XJn(Ynp?) COS nO EXP (14) 
and for medium three: 


(s) COS Ym3(/—2) 


sin Ym3(/—2) 


Xexp (iwmnst). (15) 


The upper or lower line in the cos, sin applies 
according as s=1 or s=2. 

The possible natural modes are described by 
the parameters n, p, m, in that order; where 


n=0,1, 2, --+ is the order of Tnlvoor), p is the 
particular root for a given order m, and m 
depends on the length of the resonator; the 
parameters p and m may also take on any posi- 
tive integral values. Thus, the Eo1: mode may be 
obtained from Eqs. (13)—(15) above by letting 
n=0, p=1, and m=1. 


(8) 
The expressions v,, are determined from the 


equations : 


=0; Jn'(vnp@)=0, (16) 


where a is the radius of the cylindrical resonator.” 
The following equations, obtained from the 
boundary conditions, suffice to determine the 7’s 


and 13 
Mi€1 Me€2 
_ (9)? 


Ca=1/(uaea)', 


12 It is interesting to note that Eqs. (16) have p-infinite, 
discrete, positive roots, all of which are mutually distinct. 

13 Henceforth we shall omit subscripts wherever there is 
no ambiguity in so doing. 
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where the subscripts refer to the medium, and the 
superscripts distinguish between the electric and 
magnetic modes, as before, with ca as the velocity 
of electromagnetic waves in medium d. The 
fourth equation is 

tan yi%,-+tan 


a; tan ml, tan 7212 — 1 


s=1, 2, (18) 


where the following abbreviations are used, for 
s=1: 


and for s=2: 
ag = psy2™/poys™; (20) 


An important special case is that of two dielec- 
tric media enclosed by the perfectly conducting 
cylinder. To obtain the boundary condition 
equations, allow /; (or /;) to approach zero. Then 
(17 and 18) become 


Me€2 
tan —ay tan yi). - (22) 


Expressions for the fields may be obtained by 
substituting (13)-(15) into (6)-(9), and the 
charge and current densities are then determined 
by (10) and (11). The relative field amplitudes in 
the various media may also be obtained without 
difficulty from the boundary conditions. The 
above analysis is easily generalized to include 
cross sections of the resonator other than circular, 
the only modification occurring for the terms 
involving the transverse dimensions. 


=(w)*, (21) 


and 


Ill. CALCULATION OF THE RESONANT LENGTHS 
AND NATURAL MODES 
Equations (18) are the expressions relating the 
resonant lengths /;, /2, 13; for the E and H modes. 
It is immediately obvious that one possible set of 
solutions is 


d=1,2,3. (23) 


This will be so only if 

cr (v)2] = + 

7 (24) 
hold simultaneously. 
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For m>0O we have modes higher than the 
fundamental.* (In the calculations and discussion 
of the experimental work that follows we shall 
be concerned primarily with the modes for which 
m=1, and hence the subscripts for m on the 7's, 
may be omitted.) If m were equal to zero, then 
the propagation constants ya for the three 
media would also be zero, which is immediately 
obvious from Eq. (23). Now w does not depend 
on any particular medium, and generally the 
various velocities, 1/(uzea)*, are different, so we 
see at once from Eq. (17) that ya“ must be 
different from zero, and hence we conclude that 
tf there is more than one dielectric present in the 
resonator, it is not possible to excite modes of the 
type Enpo, Hnpo. This is physically plausible since 
we would not expect to be able to excite modes 
independent of the thickness and position of the 
dielectrics. Such modes may be excited, of course, 
if there is only one dielectric in the cylinder. 

The set (23) is special in that (24) must hold 
also. But in general, the restrictions on /;, ls, J; 
are such that (24) will not be satisfied. The pro- 
cedure is to choose any two of the lengths /,, J, 
l; and a frequency, whence the third resonant 
length is determined from Eq. (18) for s=1, 2, 
depending on the type of mode excited. The 
quantities ya“ are determined from (17) and are 
the same for the three media as well as for the 


two. The upper limit on we is given by rm/la, 
and the lower cut-off limit occurs when 


—(Yap) =0, d=1, 2,3. (25) 


The lower cut-off frequency is determined by 
the medium having the lowest dielectric con- 


stant, i.e., for the value of serfs which is 
smallest and still satisfies (25), which becomes 


* Note: In the simple case of a resonant hollow pipe con- 
taining only one dielectric medium, the wave-length is 
given by \=2x[(ym™)?+(vap™)*]-#, and the propagation 
constant ym is determined from ym“mzx/l, where m=0, 
1, 2, ---, and 1 is the over-all length of the resonator. When 
the pipe contains several different dielectrics, however, we 
can no longer use the term wave-length in the pipe, and 
‘ym is then not in general equal to mx/I. This is use 
there are now as many different wave-lengths, and fractions 
thereof, in the tube, as there are different dielectric media 
and dielectric interfaces. In order to carry over symbolically 
the dependence of the higher modes on y,, in analogy 
with the one-dielectric case, we have used Eq. (23) to set 
the upper and lower limits on the propagation constant; 
for example, for the mode in which we have m=3, we have 
2x/la<y2a3/la. Additional half-wave-lengths in the 
various media now no longer define the modes, as far as 
dependence on the axial, or 2 direction is concerned. 
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Fic. 2. Resonant length of one air , 4, as a func- 
tion of /,), a chosen resonant length the other air space, 


od a given thickness of dielectric ]2, for one and two half- 
ve-lengths in the tube, g=0, 1, respectively. 


simply 
(s) (s) 
Wmn p(Cut-off ) = Vnpla- (26) 


(a). Cavity Resonator of Constant Over-All 
Length Containing a Movable Slab 
of Dielectric 


For three dielectric media an important case to 
be considered is that for which medium one is the 
same as medium three, with a thickness /, of 
dielectric separating them. The total length of the 


resonator is kept constant, and the position of the 


dielectric plug is varied. 
For aconstant tube length, /, and since medium 
one is the same as three, we have 


(27) 
a2 =a, =a: (28) 


and 


let 


where /, and 6 are constant. Substituting these 
values into (18) and solving for /,“, we obtain 
finally 


ean-( Bo + P—4Ac 
2A 


+rq]_ (30) 


q=0, 1, 2, aes 


where integral values of g>0O correspond to 
additional half-wave-lengths, g, of course, being 
such that /,;“+/,+/;@=1, which is fixed. We 


have also 
AM=a™ tan y16+(a™)? tan 
BO =[1—(a™)?] tan yi tan (31) 
tan yes. 


Thus there are two values of the resonant 
length /,, as might be expected, one for a position 
near one end of the tube, the other near the other 
end. The sum of these two values is constant and 
equal to 6, and the variation of 1,“ with fre- 
quency is symmetrical about /, equal to /;“, as 
can easily be shown (see Fig. 13). This indicates 
that the greatest effect of the dielectric on 
changing the resonant frequency occurs, as one 
would expect, at an anti-node of the E field, or in 
the middle of the tube, as shown in Fig. 13, if the 
resonator is a half-wave-length long (¢g=0). When 
we have (B®)?=44C, and this 
condition becomes 


Equation (32) is a complicated function of fre- 
quency, but the frequency at which /,“ equals 
1; can be obtained fairly readily by graphical 
methods. An upper limit to resonance for the 
particular modes excited is given when /;“ (or 
l,) is zero; then the problem reduces to that of 
two dielectric media. 

It was not possible to test (30) directly by 
experiment because an oscillator variable over a 
wide enough frequency range was not available. 


(b). Cavity Resonator of Variable Over-All 
Length Containing a Movable Slab 
of Dielectric 


It was found most convenient for the experi- 
mental work to select J,“ as the independent 
variable in (18), Js as the dependent, for a given 
dielectric thickness /,, and a given constant fre- 
quency. The over-all resonant length varies as 
1, and J, are changed, while the resonant 
frequency remains at some constant value; in the 
case of the theoretical curves (Fig. 14), this 
frequency was chosen as 3000 mc/sec., or a 10-cm 
wave-length; in the actual experimental work it 
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Fic. 3. Amplitude of the standing wave pattern vs. position 
in the resonator. 


was 3070 mc/sec., or about 9.78 cm. Let 
$=tan (a tan yh); (33) 
then from (18) we have 


(8) 4. 
tan-*( 


34) 


These equations (34) are general for three 
different media. 

The important case is as before, that media one 
and three be the same, a condition which 
simplified the experimental work. Then y,“ 
=73, etc., and (34) is symmetrical about 
1, =1,;. (To test this, replace /, by /;“ in the 
above.) Again the greatest effect of the dielectric 
occurs when the slab is symmetrically placed 
with respect to the ends of the tube, if it is half 
a wave-length long, as shown in Fig. 14. When 
the resonator is more than a half-wave-length 
long, the maximum effect is obtained as before at 
antinodes of the E field. Figure 2 shows this for a 
full wave-length in the resonator. 


(c). Cavity Resonator of Variable Total Length 
Containing a Slab of Dielectric at One End 


This is the case of two dielectric media; if we 
let 3=0, then Eq. (22) is easily solved for 1, ; J. 
is as before the thickness of the dielectric slab. 

When the resonator is \/2 cm long and when 
1, is zero, or equivalently when 


then this sets an upper limit to the frequency 


1000} ) Hy, mode 

B00 Theoretical 
Experimental 

‘0600 at 3070 


« -OjJoo +0200 "0300 


Fic. 4. Resonant length J; as a function of dielectric 
thickness /, at 3070 mc/sec., experimental and theoretical 
curves for the Hi; mode, with the dielectric at one end of 
the resonator. 


obtainable for a given size of dielectric and for a 
given mode (see Fig. 9). When this upper limit is 
equal to the lower cut-off frequency, we have the 
maximum value of /;“ possible and still excite (in 
this instance) the Eo.; or Hi, mode. This means 
that yi“ =0 at cut-off, and c,>c2; putting this 
into (22) we obtain ~ 


1, (max.) = (q+1)x/y2 
(36) 


and g>0 for more than one-half wave-length in 


the tube. The lower cut-off ‘requency is the same 
for all thicknesses of dielectric equal to or less 
than 7,“ (max.). If we combine Eqs. (26) and 
(35), letting Af be the band of frequencies for 
which a given mode may be excited, we find that 

which depends on the dielectric thickness; thus, 
by inserting different lengths of dielectric, the 
transmission band for a particular mode may be 
adjusted between the upper limit, given by Eq. 
(35), and the lower cut-off limit, given by 
Eq. (26). Under these conditions the resonator 
acts like a band transmission filter, for the given 
mode excited. 


IV. RESULTS 


The general theory has been developed in 
Section II, and the analysis for three special cases 
has been given in Section III. A brief discussion 
of the experimental and theoretical results 
follows. See Figs. 3-14. 
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(a). Experimental Procedure and Results 


With a dielectric plug of thickness /, at one end 
of a hollow cylinder and a movable piston at the 
other, the resonant air length /,“ was determined 
for several values of J,. A klystron oscillator, 
loosely coupled to avoid forced coupled-circuit 
effects, was used to excite the resonator. Because 
of the presence of the antenna and detector, as 
well as a small longitudinal slot down which the 
antenna and detector could be moved, and be- 
cause of the fact that the resonator was not a 
perfect conductor, the fields were slightly dis- 
torted, and the maximum and minimum values 
of the E and H fields did not occur quite at the 
ends of the cylinder, as predicted by theory. The 
effect of these factors was to decrease the over-all 
resonant length by 0.48 cm, the correction which 
was added to /, to obtain the first calculated 
half-wave-length value of 7.89 cm. The calcula- 
tions for this case were made from Eq. (22). The 


H,, mode 
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Fic. 5. Resonant length /;° as a function of a chosen 
resonant air length /,@ for two different thicknesses of 
—s" l:, at 3070 mc/sec., for a tube length greater than 
Atube/2- 
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Fic. 6. Propagation constants, y:, v2 vs. frequency, with 
(air), €2=2.60eo, Hin. mode. 


average experimental value of \/2 in the tube 
was found to be 7.90 cm, good agreement with 
the theory. The dielectric used was a cylindrical, 
tightly-fitting slab of polystyrene, and the meas- 
urements were made at a wave-length of 9.78 cm, 
corresponding to a frequency of 3070 megacycles/ 
sec. for the Hy; mode. 

Figure 3'* shows part of the resonance curve 
obtained for the Hi; mode in the hollow cylinder, 
of radius a=3.65 cm, at a frequency of 3070 
mc/sec. Figure 4 compares the experimental and 
theoretical values of 1, vs. lz, for various values 
of lz, at 3070 mc/sec.; the theoretical curve was 
calculated from Eq. (22) for ¢:=2.41¢, which 
was the experimentally determined value for the 
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Fic. 7. Resonant air-length 1, as a function of frequency 
for different thicknesses of dielectric /2, for a half-wave- 
length in the tube, Hin mode. 
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Fic. 8. Resonant air length /,® as a function of dielectric 
thickness /, for various resonant frequencies, for a half- 
wave-length in the tube, H,; mode. 


“4 Note: the bar, (—), appearing over any symbol in the 
figures, denotes the value of that quantity for the Hi or 
magnetic mode, s = 2. 
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Fic. 9. Thickness of dielectric as a function of the upper 
cut-off frequency for one-half wave-length in the tube. 
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Fic. 10. Propagation constants 7:, y2 as a function of 
frequency, with ¢: (air) ¢:=2.60e0, Eo: mode. 


dielectric used. The experimental work is seen to 
be in good agreement with the theory. 

Figure 5 below shows typical experimental 
curves for the three dielectric case (medium one 
the same as medium three), where the air length 
1; is plotted vs. 1;®, for a dielectric of thickness 
I, at various distances /,” from one end. These 
curves are similar to the theoretical ones shown in 
Figs. 2 and 14. Figure 14 was obtained by 
calculation from Eq. (42). 


(b). Properties of The Cylindrical Resonator 


Equations involving the resonant lengths and 
frequencies of cylindrical cavity resonators con- 
taining two and three dielectric media have been 
developed in Sections II and III. The properties 
of such resonators indicate their usefulness as 
ultra-high frequency resonant circuits which are 


tunable by varying the length of the air space or 


by shifting a dielectric plug. Figure 13, calculated 
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Fic. 11. The same as Fig. 7, Zo: mode. 
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Fic. 12. The same as Fig. 8, Eo.; mode. 


from (32) for the Hi: mode, shows the variation 
of resonant length /,® with frequency, a cylinder 
of dielectric being moved about within the tube, 
whose total length is constant. As we have seen, 
a half-wave-length resonator containing a die- 
lectric at one end or at some position about the 
middle of the cylinder acts like a filter, for a given 
mode, whose frequency transmission band width 
is adjustable by varying the dielectric thickness. 
Figure 9, calculated from Eq. (35), shows the 
variation of dielectric thickness /, as a function of 
the upper cut-off ‘frequency, for one-half wave- 
length in a tube of radius a = 3.65 cm, for both the 
Eou and Hi, modes; from such a curve it is a 
simple matter to obtain the thickness of die- 
lectric needed for a given band width, and vice 
versa. Resonators containing a greater number 
of half-wave-lengths may also serve as band 
transmission filters, but may be impractical be- 
cause of their size and the length of dielectric 
required for a given passband. 
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(c), Measurement of Dielectric Constants 
at 10 cm 


The tube containing two different dielectrics, 
for example, polystyrene at one end, the rest of 
the interior filled with air, suggests a simple 
absolute method for measuring the dielectric 
constants at microwave frequencies, in this 
case approximately 10 cm. 

Figure 8 shows calculated curves of the reso- 
nant air length /, as a function of the dielectric 
thickness /., for one-half wave-length in the tube, 
with frequency as a parameter, and ¢.=2.60¢; 
these calculations were made from Eq. (22). 
Figure 7 is obtained from this curve, where the 
parameter is now /, and the resonant air length 
1, varies as a function of the natural frequency 
of the resonator. 

To find the dielectric constant of a sample of 
the unknown, thickness /2, all that is required is 
to place the sample at one end of the tube, and to 
determine the new resonant air length /,®, then 
determine the resonant length with air only, and 
plot tan and —tan as a 
function of 72. The point of intersection of these 
two curves determines y2”, while y:® is readily 
calculated for air from (21) and the frequency at 
which the data are taken. The H1,, mode was the 
one used experimentally, though any physically 
possible mode will do. Equation (21) becomes 


(v2)? + (rrr)? J/ (wir) = €2 = 
meters/sec., (38) 


me 


Fic. 13. Resonant air length 1, as a function of fre- 
quency, /;= 1.0 cm, in a resonator of constant total length 
1=.074 m, for one-half wave-length of the Hi., mode. 


Fic. 14. Same as Fig. 5: calculated curves for one-half 
wave-length in the tube, at 3000 mc/sec. 


from which é2 is readily found. The method is 
absolute, as it depends only on the dimensions of 
the resonator, the sample, and on the wave- 
length in the tube. By this method, with data 
taken from Fig. 4, ¢2 for polystyrene was found 
to have the value 2.41¢9, which is in fairly close 
agreement with Lamont’s value'® of 2.47¢, the 
difference being due to the fact that different 
samples of polystyrene have different dielectric 
constants, because of variation in purity and 
polymerization. 

This method of measurement is convenient for 
solids, as only a thin sample, 1.0</,<8.0 cm, is 
needed in most cases, and the larger the dielectric 
constant the thinner the sample. As may be seen 
from Fig. 8, the sample should be of the order of 
thickness of a quarter-wave-length or greater, so 
that an appreciable resonance shift will occur. 
With a little adaptation, « of liquids can also be 
measured, one possibility being to place the 
resonator in a vertical position, with the bottom 
end sealed to prevent the escape of any fluid; this 
does away with the need of a separate fluid con- 
tainer ; here, too, a few centimeters or so of liquid 
should be sufficient. . 

The author wishes to express his sincere ap- 
preciation to Professors Ronold King and J. H. 
Van Vleck for their generous encouragement and 
many helpful suggestions in the preparation of 
this manuscript, and to Professor King also, for 
the use of his apparatus in the experimental work. 


1H. R. L. Lamont, Phil. Mag. 30, 15 (1940). 
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The Effect of High Electrostatic Fields upon the Vaporization of Molybdenum* 


G. B. EstaBroox** 
University of Pittsburgh, Pittsburgh, Pennsylvania 
(Received March 6, 1943) 


It has been found for wires of varying diameters, (1) that, for field strengths of about 0.52 < 10° 
volts/cm or less at the wire surface, there was no change in the rate of change of resistance with 
time and hence the rate of vaporization upon application of the field; (2) that, for field strengths 
above approximately 0.52 X 10® volts/cm, as successively higher fields were applied, the rate of 
vaporization first decreased slowly, then more rapidly, and finally asymptotically toward zero 
at field strengths over 1.8 X 10* volts/cm; (3) that, with increasing temperatures, the same field 
strength applied at the surface of a wire produced progressively smaller decreases in the rate of 
vaporization; (4) that, when there was an effect upon the rate of vaporization, sudden increases 
in resistance occurred when the field was applied, and sudden decreases when it was removed. 
Worthing has reported similar results for tungsten except that the sudden increases and de- 
creases were opposite to those of platinum and molybdenum. 


INTRODUCTION 


HE filament of a tungsten lamp operated at 
its normal: temperature of about 2450°K is 
known to decrease in cross-sectional area with 
time. Evidence of this change may be seen by 
the deposition of metal upon the sides of the 
lamp. This change, because of vaporization of 


L_] 


Fic. 1. Tube used for mounting molybdenum wires. 
A, copper tube, 2.8 cm I.D., walls 0.2 cm thick; B, Nonax 
glass tubing; C, molybdenum wire; D, 20-mil tungsten lead 
wires; E, small spring of 5-mil tungsten wire to provide 
slight tension; F, liquid-air tube; G, side tubes containing 
tungsten electrodes, on either side of main tube, for 
“cleaning up”’ residual gases. 


* This work was carried out during the year 1931-32. 
Publication has been delayed because of an apparent dis- 
crepancy with data on the vaporization of molybdenum 
already published. Now that an explanation of this 
discrepancy has been found (see following paper), it is 
offered for publication. 

** Now at the Department of Physics, University of 
Maryland, Baltimore, Maryland. 


the metal, is accompanied by an increase in 
resistance of the wire. 

In an investigation of the effects of high elec- 
trostatic fields upon the vaporization of small 
heated tungsten wires Worthing! has found that 
for a wire 0.020 mm in diameter, operated at 
2500°K, a radial field of strength 2.7 X 10° volts 
per cm at the wire surface reduced the rate of 
increase of resistance, and hence the rate of 
vaporization, to 45 percent of the value without 
the field. He further observed that, upon apply- 
ing the field, there was a sudden decrease in 
resistance of the wire and that, upon removing 


Fic. 2. Electrical arrangement for measuring the rates of 
vaporization of molybdenum wires. A, tube containing 
molybdenum wire; B, standard resistance; C, milliam- 
meter ; D, variable rheostat ; E, storage batteries; F, switch; 
G, potentiometer; H, voltage multiplier; Z, Wimshurst 
machine; J, electrostatic voltmeter; K, single pole double 
throw switch for grounding the circuit containing the 
molybdenum wire. 


1A. G. Worthing, Phys. Rev. 17, 418 (1921). 
* E. H. Greibach, Phys. Rev. 33, 844 (1929). 
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.Fic. 3. Change in the rate of vaporization of a molybdenum wire (r = 0.00243 cm, 
T= 1462°K) upon application of four different electrostatic potentials. (Ordinate 


scales differ by a constant.) 


the field, there was a sudden increase. Greibach?* 
has predicted, on the basis of classical theory, 
that the effect of an electric field, such as used 
by Worthing, should be to decrease the rate of 
vaporization by less than one percent. 

The present investigation was undertaken to 
obtain information on the behavior of small 
molybdenum wires under conditions somewhat 
similar to those described for tungsten. Prior to 
using molybdenum wires, a similar study of small 
platinum wires was attempted. A decrease in the 
rate of vaporization was observed when they 
were under the influence of the electric field. 
Contrary to the case for tungsten, a sudden 
increase in the resistance of the wire took place 
at the instant of applying the field and a decrease 
when the field was removed. However, the data 
were not consistent, probably on account of the 
difficulty in out-gassing the platinum; and so the 
work on platinum was dropped temporarily. 


APPARATUS 


The molybdenum used was obtained from the 
Cleveland Wire Works of the General Electric 
Company at Cleveland, Ohio, in the form of 
wires 0.00486 cm, 0.00684 cm, and 0.00762 cm 
in diameter. While no analyses of the particular 


samples used were available, the Cleveland Wire 
Works reported that analysis of powdered metal, 
similar to that from which the wires were made, 
ran 99.9 percent molybdenum. A fourth size 
wire used was obtained by electrolytically etch- 
ing the 0.00486-cm wire to 0.00458 cm. The wires 
were examined under a microscope and initially 
found to have polished surfaces with no pits. 
Their diameters were determined by weighing a 
measured length on a microbalance. Approxi- 
mately 10-cm lengths of these wires were mounted 
in highly evacuated tubes, as indicated in Fig. 1. 
During the evacuation, which took about 18 
hours, the tubes were heated in a furnace at 
about 425°C. 

Figure 2 represents the electrical arrangement 


for measuring the change in resistance of the - 


wire while keeping the current constant. For 
convenience of measuring potential differences, 
a Leeds and Northrup potentiometer, modified 
to operate as a limited double potentiometer, as 
described by Worthing and Forsythe,’ was used. 
The electrostatic potential difference was meas- 
ured by means of an electrostatic voltmeter as 
described by Shrader.* 

3A. G. Worthing and W. E. Forsythe, J. Opt. Soc. Am. 


10, 599 (1925). 
‘J. E. Shrader, J. Opt. Soc. Am. 6, 273 (1922). 
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Fic. 4. Variation in the rate of vaporization of molybdenum wires with electrostatic 
field strength applied at the surface of the wires. 


PROCEDURE AND RESULTS 


Before making the initial run on any sample, 
the wire in the tube was first aged for a period 
of several hours, at a temperature above 2000°K, 
and then run at the desired temperature for 
approximately two hours, or until conditions 
became steady. Following this, by repeated ad- 
justments, the current in the wire was kept 
constant to about one part in 200,000 and differ- 
ences in potential between the ends of the wire 
measured every five minutes. By use of the 
limited double potentiometer mentioned above, 
it was possible, once adjustments were made, by 
merely throwing a single switch, to change 
directly from the current balance to the potential 
balance without disturbing the galvanometer. 
Thus the two measurements could be made in 
rapid succession, a condition essential in this 
work. 

Figure 3 contains a typical set of curves show- 
ing the changes in the rates of increase of voltage 
with time, for a wire 0.00486 cm in diameter 
operated at an average temperature of about 
1462°K, upon the application of several high 
potentials between the cylinder and the wire. 
To obtain a single curve, three steps are neces- 
sary. Consider the curve for 17,000 volts. After 
the preliminary heating, the wire was operated 
at constant current for a period of 50 min., with 


a potential of 45 volts only between the wire and 
the cylinder, to eliminate thermionic emission 
effects. During this time, the potential difference 
between the ends of the wire was found to in- 
crease regularly, so that when voltage was plotted 
against time, the slope of the curve, that is, 
the rate of increase of potential difference, was 
constant, and in this instance had a value 8.00 


; 


Fic. 5. Relative change of vaporization rate of a 
molybdenum wire (r=0.00243 cm) with temperature. 
Electrostatic field 1.03 x 10* volts/cm. 
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Fic. 6. Sudden changes in potential difference between the ends of molybdenum wires 
upon application and removal of high electrostatic fields. 


X10-* volt per min. The 45-volt potential 
difference then was changed suddenly to 17,000 
volts and the changes in potential drop along the 
wire followed for one hour and 15 min., the rate 
being 4.00 10-* volt per min. Finally the high 
potential was removed, the 45 volts reapplied, 
and again changes in potential drop along the 
wire were followed for 50 min. The slope of the 
curve returned to 8.00X10-* volt per min. 
In many cases, this procedure was continued 
through two more steps, making five in all. 
Since in each run the current was kept constant, 
the slope of the curve is proportional to the rate 
of change of resistance of the wire. Now as a 
wire evaporates, its cross-sectional area dimin- 
ishes; and, since this varies inversely with re- 
sistance, a change in the rate of increase of 
resistance is interpreted as a change in the rate 
of vaporization. As appears later in the dis- 
cussion, there is a disturbing factor in this 
interpretation. 

The slopes of the curves before and after the 
application of the high potentials were always the 
same for any one run. For the runs graphed in 
Fig. 3, the application of a potential of 8000 volts 
resulted in no change of slope. A slight change 
of slope was observed for 10,000 volts, more of a 
change for 17,000 volts, and for 20,000 volts the 
curve became almost horizontal. In addition, for 
all potentials above 8000 volts, a sudden increase 


in potential difference between the ends of the 
wire was observed when the high potential was 
applied, and a decrease, though not quite so 
great in magnitude, when it was removed. 
Similar runs to those just described were made 
on the same wire with other high potentials and 
at other temperatures, and on other wires of 
different diameters. Curve 1 of Fig. 4 is a com- 
posite curve showing the results obtained from 
four sizes of wire operated at about 1462°K. In 
this curve, the percentage variation in the rate 
of vaporization is plotted against field strength 
at the wire surface. For reasons appearing later 
under “Discussion,” field strength has been 
chosen as abscissas rather than potentials. As is 
evident from the curve, for a given field strength, 
the change in rate of vaporization is independent 
of the size of the wire. No change of rate of 
vaporization was found for field strengths of 
0.52X10* volts per cm or less. However, for 
increasing field strengths above this value, the 
rates of vaporization, upon applying the field, 
decreased, at first slowly, then more rapidly, 
and finally asymptotically toward zero at field 
strengths above 1.8X10® volts per cm. Curve 2 
shows similar results for an average temperature 
of 1644°K. As for 1462°K, no change was found 
for field strengths of 0.52X10® volts per cm or 
less. At higher field strengths, the vaporization 
rates are higher for the higher temperature. 


| 
fa 


356 G. B. ESTABROOK 


Figure 5 represents the variation of the rate of 
vaporization for the 0.00486-cm wire with tem- 
perature, the field in all cases being the same. 
Again, the higher the temperature of the wire, 
the less the field decreases the vaporization rate. 

In Fig. 6, curve 1, the sudden increase in re- 
sistance of the wire upon application of a field 
is plotted against field strength at the wire 
surface, and in curve 2, the sudden decrease upon 
removing the field. The values in the former case 
are somewhat larger than those in the latter. 
Both curves, however, have similar shapes. In 
neither case was an effect observed for field 
strengths of 0.52 10° volts per cm or less. 

To determine whether or not the change in re- 
sistance could be ascribed to vaporization, ap- 
proximately 20-cm lengths of the wires were 
mounted in tubes and operated for 50 hours or 
more at an average temperature of 1462°K, the 
average temperature used in most of the runs 
just described. From the weights of a measured 
length of wire before and after the prolonged 
heating, the amount of metal vaporized was cal- 
culated and compared with that expected from 
the resistance changes in the runs on the wire 
from which the data of Fig. 3 were obtained. For 
a 0.00486-cm wire, which had first been aged at 
a temperature above 2000°K for about an hour 
and then weighed and put in operation, the ratio 
of the loss in weight of the wire to its original 
weight was 0.0253, while for the same sized wire 
for which resistance variations were followed the 
rate of change of resistance (dR/di) was 114.3 
X10-* ohm per min. In this latter case, the 
average value of (1/R)(dR/di) was 6.16X10-® 
min.—!, and the value of (1/R)(dR/dt)t was 
0.0211, where ¢ is the total time of operation of 
the wire. On the assumption of no pitting or 
roughening of the surface during operation, this 
value is subject to a correction for change in 
temperature. This is necessary to give the frac- 
tional change in resistance had the initial and 
final resistance measurements been made at the 
same temperature. When thus corrected, a value 
of 0.0147 is obtained. (See Appendix.) If one 
considers the uncertainties resulting from (1) the 
difficulty in obtaining small differences in weight 
with a microbalance (the difference was 0.068 
mg), (2) the end losses (the vaporization rate was 
less at the cooled ends of the wire than at the 


center), (3) the slight roughening of the surface 
of the wire with operation (the roughening of a 
surface by increasing the radiating area cuts 
down the increase in temperature which normally 
follows operation of such a wire on constant 
current), and (4) the possible decrease in re- 
sistivity, or in radiancy, or both (this follows a 
decrease in V*J, which Worrell® has shown to be 
true for a molybdenum wire mounted in a highly 
evacuated tube), the ratio 0.0253 seems not in- 
consistent with the ratio 0.0147, and it seems 
reasonable to attribute the change in resistance 
to vaporization. 

It is interesting to compare the rate of vapor- 
ization which may be obtained from the above 
with that reported by Langmuir and Mackay.*® 
They give for molybdenum the following ex- 
pression : 


log m= 17.11 —38600/T —1.76 log T, 


where m is the rate of vaporization in g/(cm? 
sec.), and 7 the temperature in °K. On the 
scale of temperature used by Langmuir and 
Mackay, 1462°K becomes 1443°K. Substituting 
this temperature in the formula, one obtains a 
value of 6.31X10-" g/(cm?sec.) for m. The 
value obtained from the above run is 1.36 10-® 
g/(cm? sec.). This difference is discussed in the 
following paper by Reid. 


DISCUSSION 


To account for the decrease in the rate of 
vaporization of the various sizes of wire with the 
application of a strong electrostatic field, it may 
be assumed that the atoms which have left the 
surface of the wire in a high electrostatic field 
become dipoles and that many of them are 
forced back by the field to the surface of the wire 
where the field is strongest. Consider a dipole of 
length / and charge e in a radial field of strength 
R between two concentric cylinders, the inner 
cylinder having a radius a, and the outer cylinder 
a radius 


«-_r—to positively charged wire. 


5F. T. Worrell, Phys. Rev. 61, 520 (1942). 
( —— and G. M. J. Mackay, Phys. Rev. 4, 377 
1914). 
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The force on the dipole is accordingly 
F=(R+41(9R/dr) Je—[R—41(0R/dr) Je 
=le(dR/dr). 


Let us assume, as is usually done, that the dipole 
moment /e is related to the field strength R, thus 


aR=le, 


where a@ is the polarizibility. Now the work done 
by a vaporized atom at the expense of its kinetic 
energy against the field in moving from the wire 


to the cylinder is 
W= f Fdr, 


or substituting the above values, we have 
b b 
w= f le(@R/ar)dr =a f R(@R/ar)dr. 


Substituting for R and dR/dr the values 
= — V,/r In (b/a) 
(b/a) J, 


where V, is the potential at the outer cylinder, 
the inner cylinder being grounded, we obtain 


—aV,? aV,? 1 1 
w-——— f 
b\*Lb? a? 
( In -) 2{ In -) 
a a 
To produce equal decreases in the rates of 
vaporization for two wires of different radii, at 
the same temperature and placed concentrically 
in cylinders of the same size, requires the same 
work per dipole for the transfer from a to b. 
Thus, where W, is the work in the case of a wire 
of radius a1, and W, the work for the same 
effect upon the rate of vaporization in case of a 
wire of radius do, 


Wi = Vox[In (b2/a2) — 


W: In (61/a) — (1/2) 


Neglecting the 5? term since } is large compared 
with a, it follows that 


Wi) W2= Rai/Raz=1, 


where Ra; is the field strength at the surface of 
the wire of radius a:, and Rez that for a wire of 
radius a2. Hence for equal field strengths at the 
surfaces of the wires, equal decreases in the rates 
of vaporization are expected. This is in accord 
with Fig. 4. It is to be noted that though the 
rate of vaporization is proportional to the field 
strengths, the existence of a dR/dr is the charac- 
teristic of the field which produces the change. 
In the case of the sudden increase in resistance 
observed when the field was applied, or decrease 
when it was removed, it was first thought that 
this might result from a current between the 
wire and the surrounding cylinder due to either 
a positive ion current from the wire or a leakage 
current over the surface of the glass. To test for 
such a current, the effect of switching the ground 
connection from one terminal of the standard 
resistance B, Fig. 2, to the other was investi- 
gated. No indication of such a current was found. 
If a current had been present, due to its inclusion 
or exclusion from that passing through the 
standard resistance, it would have been detected 
in the potential measurements. While potential 
differences could be read to about one part in 
200,000 and a sudden increase in resistance 
caused a change of as much as one part in 10,000, 
the change in the ground connection never pro- 
duced any detectable effect. Further, since the 
introduction of grounded guard rings did not 
alter the sudden increase and decrease of re- 
sistance, the possibility of leakage currents 
seemed excluded. Again since a positive ion 
current, if present, would be expected at fields 
lower than that which just produced the sudden 
changes in resistance, the possibility of such a 
current being present seemed likewise excluded. 
As mentioned above, there is a disturbing 
factor in the interpretation of the rate of in- 
crease of resistance as proportional to the rate of 
vaporization. It is the difference in the sudden 
changes of resistance when applying and re- 
moving the electric field. During a run of three 
steps, as described earlier, the total change in 
resistance is not equal to the sum of those 
occurring during the individual steps, the differ- 
ence being that between the two sudden changes. 
Increases and decreases of resistance in the 
same direction were observed for platinum, but 
for tungsten Worthing reported opposite changes. 
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These changes are further discussed in the fol- 
lowing paper. 


APPENDIX. CORRECTION OF OBSERVED 
(1/R)(dR/dt) FOR CHANGE 
IN TEMPERATURE 


Consider the resistance of an electrically heated wire in 
vacuum as a function of its temperature 7, and its radius r, 
thus 

R=f(T,r). 
fo tes wire evaporates, its resistance changes with time ¢, 
an 
dR/dt =(8R/dT),(dT (1) 


[(@R/R)/(8T/T),] = (T/R)(AR/dT), =B. 
Multiplying (1) by (1/R) and substituting, we obtain 


Ra di di (2) 


Now the energy dissipated from a wire of length /, carrying 
a current J, is 


R (3) 


where p is the resistivity, and ® the radiancy expressed in 
watts/cm*. Then 

p/ R= 2x*r?/I? = (p0/Ro)(T/T) (4) 
where p=po(T/To?; R= Ro(T/To)* and po and Ro are 
fixed values at JT». The terms 8 and x are taken as constant 
values for the small interval of temperature considered. 
From (4) 

= Ro) (T/T)? (5) 

and 


d(In r) /dt = (8—x/3)[d(In T) /dt]. (6) 
Substitution in (2) gives 
1/R(dR/dt) =[B— 3(8—x)}(1/T) (dT /dt). (7) 


For (7/R)(@R/dT),=8, and for (T/R)(@R/aT),=x, at 
1466°K, Worthing’ gives, for molybdenum, 1.145 and 5.14, 
respectively. Substitution in (7) and multiplication by ¢ 
gives 
(1/R) (dR/dt)t = 3.80(1/T) (dT /dt)t. (8) 

Since for the region considered (dR/dt) and (d7/dt) are 
constant, except for minor fluctuations, we may replace 
(1/R)(dR/dt)t and (1/T)(dT/dt)t by AR/R and AT/T 
where AR and AT are the changes in R and T during the 
interval ¢. 

For the run described in ‘Procedure and Results,” 


(t/R)(dR/dt) = (SR/R)o=0.0211, (9) 


the observed value which is to be corrected for variation in 
temperature during the run to give the fractional changes in 
resistance had the initial and final resistance measurements 
been made at the same temperatures. 

Substitution of the value 0.0211 in (8) gives AT/T 
=0.0056. From this it follows that the average tem- 
perature is equal to the final observed temperature times 
[1—4(A7/T)], and in this case the final observed tem- 
perature of 1466°K yields an average temperature of 
1462°K. 

Upon applying this correction to (AR/R)o, the corrected 
value (AR/R) becomes 


Ry(1 —0.00288) — Ro(1 +0.00288) 
4(Ro+ Ro) 


where R, is the final resistance and Ro the initial resistance 
of the wire. Substitution in (10) yields for AR/R (corrected 
for temperature variation during the run) 0.0147 in place 
of the uncorrected value 0.0211. : 

The writer wishes to express his indebtedness to Dr. A. 
G. Worthing, who suggested this problem, for his help and 
encouragement during the course of the work. 


(AR/R)-= 


(10) 


7A. G. Worthing, Phys. Rev. 28, 190 (1926). 
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This paper endeavors to show that equations can be ob- 
tained to fit Estabrook’s observed variation of vaporization 
rate with electric field strength at wire surface if we con- 
sider the vaporization decrease as due to the formation of 
an adsorbed gas layer on the filament. It is mentioned that 
electric field strengths of the magnitudes used can be shown 
to have no appreciable effect upon the motion of the gas 
particles or metallic atoms in the region between filament 
and cylinder, and that time lags of the size observed can 
be approximately accounted for by the times required for 
formation and removal of the gas layer when the field is 
applied or withdrawn. The following postulates are pro- 


posed: (1) the resistance changes were due to temperature 
changes resulting from small changes in emissivity; (2) re- 
sistance changes in different directions were due to spectral 
variation in the behavior of the emissivity of a filament 
when it is coated with gas; (3) Estabrook’s vaporization 
rate was higher than that calculated from an equation 
given by Jones, Langmuir, and Mackay because the wires 
used by the latter were more coated with gas; and (4) Jones 
did not observe any effect of an electric field upon resistance 
or vaporization because his tube did not contain enough 
oxygen. These hypotheses remain of course to be proved or 


disproved by appropriate experiments. 


INTRODUCTION 


UPPOSE that we had a fine wire supported 

along the axis of a metallic cylinder, the 
whole being in a vacuum. If the wire were 
heated to incandescence by means of an electric 
current, its resistance would gradually increase 
because of vaporization of its atoms. This much 
is well known. In addition however, Worthing! 
and Estabrook? found that the application of a 
sufficiently high potential difference between 
filament and the concentric cylinder reduced the 
rate of increase of resistance of the filament. 
Their cylinders were made negative with respect 
to the filaments in order that there would be no 
emission of electrons. In addition to the change in 
slope of the resistance versus time curve, there 
was also a change in the actual resistance itself 
at the instant the field was applied, and a change 
in the opposite direction when the field was re- 
moved. These two effects were found with Pt, 
W, and Mo filaments, and presumably would 
exist also with other metals. Estabrook found 
that for Mo and Pt there was an increase in re- 
sistance when the field was applied, and a de- 
crease when it was removed. Worthing found 


.* Submitted to the Graduate School of the University of 
a in partial fulfillment of the requirements for the 
of Doctor of Philosophy. 
eaten at the United States Rubber Company, Detroit, 
igan. 
1A. G. Worthing, Phys. Rev. 17, 418 (1921). 
* See preceding paper. 


changes in the opposite direction with W. Both 
observers, however, found a decrease in the rate 
of resistance change when a radial field was 
applied. Figure 1 of this article, and Fig. 3 of 
Estabrook’s article illustrate the above. In the 
graphs it is seen that there is a rounding of the 
corners indicating that the resistance and vapor- 
ization changes do not take place instanta- 
neously. 

It seems reasonable to attribute the steady 
increase in resistance of a wire with time to a 
gradual vaporization of its atoms. If such is the 
case, then the applied radial field must decrease 
the rate of vaporization. The ratio of the rate 
of vaporization after the field is applied to that 


t WITHORAWAL OF 
APPLICATION OF 
ok 15,300 VOLTS. VOLTAGE. 
i i i i 
1 2 3 


Fic. 1. This curve shows how the resistance of a tungsten 
wire was affected by the application of a high potential 
difference between it and a concentric cylinder. The 
temperature was 2500°K, and the wire radius 0.010 mm. 
[Taken from A. G. Worthing, Phys. Rev. 17, 419 (1921). ] 
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before is given by the ratio of the two correspond- 
ing slopes of the resistance versus time curve. 
Figure 2 shows how this ratio varies with the 
field at the surface of the filament. 

In Fig. 3 we see how the resistance change is 
affected by the field strength. This figure differs 
from Estabrook’s Fig. 6 only in that the ordinates 
have been doubled, and a different curve has 
been chosen. Note that in it, as in Fig. 2, there 
is no apparent effect for fields of less than 0.5 


RA 


VAPORIZATION RATE WITHOUT FIELD 


FIELD STRENGTH AT WIRE SURFACE=€ 


Fic. 2. These curves show the variation in rate of 
vaporization of molybdenum wires with the electrostatic 
field at the wire surface. The circles are from Estabrook’s 
data. The solid lines are calculated values. 


X10® volts/em. Both Worthing and Estabrook 
found the resistance change to go more or less 
hand in hand with the vaporization change, 
although the former change has also been found 
to occur at temperatures where there is no appre- 
ciable vaporization at all. Vissat* has studied 
such resistance changes at room temperature and 
above. 

In this paper we wish to propose a theory to 
account for the phenomena mentioned. Any such 
theory should provide an explanation for the 
following features of the problem as well as those 
cited above: (1) The vaporization rate observed 
by Estabrook was much larger than that pre- 
dicted by the work of Jones, Langmuir, and 
Mackay,‘ and (2) Lloyd Jones (results unpub- 
lished), working with a very good vacuum, did 
not observe either the resistance or the vapor- 
ization change. Before undertaking our treat- 


3P. L. Vissat, Ph.D. Thesis, University of Pittsburgh 


(1941). 
I. Langmuir and G. M. Mackay, Phys. Rev. 4, 377 
Mackay, 


1914). H. A. Jones, I. Langmuir, and G. M. 
hys. Rev. 30, 201 (1927). 


ment, let us consider the dipole theory. This is 
done because the latter is apparently the explana- 
tion which occurs first to a person confronted 
with the problem. 


THE DIPOLE THEORY 


The metallic atoms leaving the filament were 
found to be uncharged, and so it is not easy to 
see why an electric field should have: any effect 
upon them unless they become dipoles under its 
influence. If they do become dipoles, then each 
atom upon evaporation will experience a force 
tending to return it to the filament. This force 
arises, of course, from the fact that we have a 
radial rather than a uniform field. Greibach’ 
studied thermodynamically the problem of evapo- 
ration of dipoles into an electric field, but was 
able thereby to account for but one percent of 
the effect observed. 

At first it appears plausible to consider the 
decrease in vaporization rate as being due to 
the influence of the electric field upon evaporat- 
ing metallic dipoles. A closer investigation, how- 
ever, reveals that an unreasonably large value 
must be taken for the dipole moment of the 
metallic atom at the wire surface in order to 
produce decreases in vaporization of the magni- 
tudes observed. It is also found that many 
phases of the problem cannot be explained suit- 
ably with such a physical picture, as, for ex- 
ample, the fact that quite large fields must be 
applied before there is any appreciable effect on 
the vaporization rate. These and other features 
make the dipole theory untenable. 


THE ADSORBED GAS THEORY 


It has been quite well established® that at 
least some of the gas particles striking a surface 
stick to it for a while before leaving again. The 
process of reflection of gas from a surface con- 
sists, therefore, at least in part, of a continual 
condensation and subsequent reevaporation of 
the particles. Since some of the molecules linger 
for a while, it follows that the surface will at 
all times be partially covered with an adsorbed 
gas film. The extent to which it is covered will 


5 E. H. Greibach, el Rev. 33, 844 (1929). 
® See, for example, L. B b, The Kinetic Theory of 
— (New York, McGraw Hill, 1934), second edition, 
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of course depend upon the conditions. In general 
gas films are driven off at temperatures much 
lower than that of 1452°K at which Estabrook 
worked, but chemisorbed films can adhere quite 
tenaciously even at temperatures considerably 
higher. The layer of atomic oxygen on the sur- 
face of tungsten filaments has, for example, been 
studied quite extensively, and is known to persist 
up to temperatures in the neighborhood of 
2000°K. Since Mo is quite similar to W, we can 
perhaps safely assume that an atomic layer of 
oxygen could adhere to it too at temperatures 
such as those used by Estabrook. 

Let us get an expression for the fraction of the 
wire surface covered by gas.’ We shall do this 
first for the case where there is no electric field, 
and on the basis of Langmuir’s physical picture. 
We denote by fo the fraction of the total number 
of adsorbable gas particles striking the filament 
that stick to it. These molecules upon condensing 
will dissociate. The separate atoms will remain 
upon the surface for a while, but will subse- 
quently leave singly, and with a heat of desorp- 
tion of ¢@o per atom. Let vp be the number of 
adsorbed atoms per unit area; 7 the number of 
Mo atoms per unit area; 0 the fraction of the 
surface covered by gas; No the rate at which gas 
molecules strike per unit area; and No’ the rate 
at which molecules of the type being adsorbed 
strike the surface per unit area. Then No’(1— 4) 
will be the rate per unit area at which the latter 
strike that portion not covered by gas, and 
No'fo(1—) will be the rate per unit area at 
which they stick. Since each gas molecule con- 
sists of two atoms, it follows that the number of 
atoms being adsorbed per unit area and per unit 
time will be given by 


2No'fo(1— (1) 


To get the number of atoms leaving® per unit 
area of the surface per unit time, let us first 
consider the case of a gas escaping through a 
small opening into a vacuum. Assuming a Max- 
wellian distribution of velocities for the gas 
particles, we get as the total number leaving the 


™See Neil K. Adam, The Physics and Chemistry of 
~~ (Oxford, Clarendon Press, 1938), second edition, 
Pp 


§ From . K. Roberts, Trans. Faraday Soc. 31, 2, 
1710 1938), 
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hole per unit area and per unit time with 
velocities greater than vo: 


Xexp (—myv?/2kTo). (2) 


If, as in our case, we have gas particles escaping 
from a surface rather than through an opening 
in a container, then we can use the above ex- 
pression for the number leaving provided that 
we let ¢9= 4$mv,? be the amount of work required 
to remove a gas particle, and replace o, the 
number of molecules per unit volume in the 
container by something proportional to the sur- 
face concentration of gas atoms, or let c—Cvo. 
This gives 


[Cvo(do+kT 0’) /(2emkTy’)*] 
Xexp (—do/kTv’) (3) 


as the number of gas atoms leaving per unit area 
and unit time. The prime is used on the 79 be- 
cause the effective temperature of the gas depart- 
ing from the filament may not<beJthe same as 
that of the surface. 


# 


FIELO STRENGTH AT WIRE SURFACE 


5 10 


Fic. 3. This graph shows the percentage changes in 
resistance obtained by Estabrook for the case when the 
field is removed. The data were obtained with molybdenum 
wires of radii 0.00229, 0.00243, and 0.00342 cm at 1452°K. 
Similar changes occur upon application of the field. 


It has been found that the number of atoms of 
oxygen in a complete chemisorbed layer on W is 
about equal to the number of W atoms on that 
surface. If the same is true for the gas adsorbed 
on Mo, then 7 will be a maximum value for vo, 
and we shall have 


00 = vo/n. (4) 


Since the rate at which atoms leave must, in the 
steady state, be the same as the rate at which 
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they stick, we set expressions (1) and (3) equal 
to each other. Upon doing so, and using (4), 
we get 


2N, o fo(2rmk 
[ 14 
exp (5) 


as the fraction of the surface not covered by gas 
before the field is applied. The fraction 1— 0, of 
the surface not covered after the field is applied 
is obtained from the last equation by changing 
the subscripts from 0 to 1 throughout. 

If we let V be the rate of vaporization of Mo 
per unit area of the filament, and assume that 
it is proportional to the fraction of the surface 
not covered by gas, we have 


Vi/Vo= (1 6:)/(1 80) (6) 


as the ratio of vaporization after field is applied 
to that before. Let us apply this equation to 
Fig. 2. 6; is equal to 4) when the field E at the 
filament surface is zero. It seems reasonable to 
imagine that 6, steadily increases as E increases, 
but yet in Fig. 2 we see that V:/V» remains 
equal to unity until E is about 0.5 X 10° volts/cm. 
Hence 6) must be essentially equal to zero if 
Eq. (6) is to apply to Fig. 2, and so Eq. (6) 
becomes 


(7) 


We next consider the number N’ of adsorbable 
molecules striking the filament per unit area and 
per unit time after the field is applied. Because 
of the fact that the gas molecules probably be- 
come dipoles under the influence of the field, it 
follows that they would experience forces tending 
to move them toward the filament. The motion 
of such dipoles in a radial field was studied 
mathematically, and led to the result 


Ni =N;'[1 +yE?/2kT3+ (8) 


where N;’ is the number of the molecules striking 
the cylinder wall per unit area and time after 
field is applied. y is the polarizability of a gas 
molecule. Let us consider it to be about 10-* cm’, 
or 1.1X10-** erg cm*/volt?. The largest value of 
E used by Estabrook was about 1.8 X 10° volts/ 
cm. 7; is the temperature of the cylinder, and 


was probably just a shade above room tempera- 
ture, or about 310°K. Using these values, to- 
gether with the Boltzmann constant k, we obtain 


/ 2k Ts) maximum = 0.00042. (9) 


Hence 
Ny =Ny’. (10) 


This means that when the field is on, the number 
of gas molecules striking the filament per unit 
area and time is essentially the same as the 
number striking the cylinder, or that field 
strengths of the magnitude used have no appre- 
ciable effect upon the motion of the gas particles, 
Since N;’ is known from the kinetic theory of 
gases, we have 


Ni’ = (11) 


where 2m is the mass of a gas molecule and p,’ 
is the partial pressure at the cylinder of which- 
ever gas is being adsorbed. 

If we evaluate Eq. (7) by means of Eqs. (5) 
and (11), we get 


[ bs fil2 exp (¢1/ (12) 


In order to use this equation, we should know 
how all of the various terms in it vary with the 
electric field strength at the wire surface. A study 
of Fig. 2, however, reveals that the data can 
most easily be fitted by a linear increase of 
¢:1/k7,' with the applied electric field. 7)’ is 
known from the experiment to have remained 
about constant, and so it is a linear increase in 
¢: which best fits Eq. (12) to the data. f would 
probably also increase with E if ¢, increased, 
but for simplicity let us assume it to be constant. 

We shall next put numerical values into Eq. 
(12). 71’ is the average temperature upon leaving 
the filament of those gas atoms which have 
adhered to it for a while, and hence is probably 
about equal to the temperature of the filament, 
or 1452°K for the first curve. As before, let us 
take 7;=310°K. » is the number of Mo atoms 
per unit area of the surface. Assuming that the 
110 plane is the one most likely to be along the 
wire surface, we get »=(14.4)10" atoms/cm’. 
The constant C, which was originally defined as 
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is® 
FIELD STRENGTH AT WIRE SURFACES 


Fic. 4. This graph shows the curves obtained when the 
heat of desorption of the adsorbed gas is assumed to vary 
linearly with the field strength at the wire surface. The 
circles are from Estabrook’s data. 


the proportionality constant in expression (3), 
can be found from an equation for the rate of 
evaporation of adsorbed atoms and molecules 
given by J. K. Roberts.’ The constant in his 
relationship was obtained from some of Lang- 
muir’s data on the rate of evaporation of oxygen 
atoms from a tungsten surface and leads to the 
value 


C= (2.56) X10’? dynes/calorie. 


According to Roberts, this value for C should be 
the same for all substances. By comparison with 
experiment, he found his equation to apply 
accurately to a considerable range of atomic 
films. Hence we shall use the above value for C. 
For f we can use the following equation given by 
Langmuir 


logie fo= 1.76—5940K°/T (13) 


for 1000°K < T < 2000°K. This expression was ob- 
tained for oxygen atoms striking a W wire, but 
might perhaps hold for our case too. For 1452°K, 
we get fo=0.0047(=f;). 

Experiments dealing with the interaction be- 
tween gas molecules and solid surfaces have led 
to many conflicting theories and contradictory 
experimental results. This is due primarily to the 
fact that the results are quite sensitive to par- 


ticular conditions existing at the surface, and to . 


the practical impossibility of knowing the nature 
of the surfaces.'® Hence the above value for f is 
quite likely to be in error. 


*T. Langmuir, J. Am. Chem. Soc. 35, 105ff (1913). 
wL. B. b, reference 6, p. 347f. 


The partial pressure p;’ at the cylinder of 
whichever gas was coating the wire is not known. 
We must guess at it, and because of the uncer- 
tainty in fi, let us consider the product p;'f; as 
being unknown, and perhaps within the range 
covered below: 


Case 1. Assume p;'f; = (10-* mm Hg) (0.0047), 
Case 2. Assume ?;’f:= (10-7 mm Hg)(0.0047), 
Case 3. Assume mm Hg) (0.0047). 


Upon putting the values above into Eq. (12), 
and applying the latter to the data at 1452°K 
in Fig. 2, we find that the equation will fit the 
data if 


0+ (4.02)10-°°E cal. cm/atom volt, (14) 


where ¢ is about (18.3)10-*° cal./atom for 
Case 1, 17.2 for Case 2, and 16.1 for Case 3. 
The coefficient of E is approximately the same 
for the three cases. In curve 1 of Fig. 4, we see 
the graph of Eq. (12) when the values given 
above are used for the various quantities in it. 
,The graph turns out to be about the same for 
each of the three cases and is essentially the same 
for an even wider range of values chosen for 
ps'fi. Thus the coefficient of E in Eq. (14) and 
the graph of Eq. (12) are not greatly influenced 
by errors in the choices for constants in Eq. (12). 
This is the case because the changes in $,+7;' 
with ¢; are small compared to the changes in 
exp (¢1/k7T;’) in Eq. (12). 

When we consider the data at 1644°K, we 
again find that a linear variation of ¢, with E 
will give fair agreement between Eq. (12) and the 
experimental results. The slope of the ¢; versus E 
line turns out to be about the same as at 1452°K, 
and as before the value for ¢) depends upon our 
choice of the constants. If we take p;’=10-’ mm 
Hg, 7,’ =1644°K, f from Eq. (13) as 0.014, and 
C, », and T; as before, we get 


$1 (at 1644°K) = [19.0+ (4.02)10-*E cm/volt 
X10-*° cal./atom. 


Using this last equation in Eq. (12), we get the 
second curve in Fig. 4. 

From above we see that ¢ appears to be in 
the neighborhood of 16 to 19X10-*° cal. per 
atom. By way of comparison we find that 
Roberts" gives a mean value of 128 kilogram 


uJ. K. Roberts, Proc. Roy. Soc. A152, 464 (1935). 
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calories per gram atom, or 21X10-*° cal. per 
atom for the heat of desorption of oxygen from 
W filaments. Other values in the literature are 
of the same magnitude. Our values for ¢o seem, 
therefore, to be reasonable ones and are of a 
size which seems to indicate that the adsorbed 
gas on Estabrook’s filaments was oxygen. 


AN EXTENSION OF THE ADSORBED GAS THEORY 


In Fig. 4 we see that the previous treatment 
does not give perfect agreement with Estabrook’s 
data. We wish now to show how better correla- 
tion can be obtained. 

We have assumed that only those Mo atoms 
which were not covered by the gas layer were 
capable of evaporating. It is reasonable however 
to suppose that a Mo atom might be able to 
leave the surface along with the oxygen atom on 
top of it. In fact, Langmuir has found that the 
combination of.a W and an oxygen atom does 
indeed evaporate from the surface of a W fila- 
ment, and has used this fact to explain the clean- 
up of oxygen by an incandescent W wire. He 
showed it to be WO; which evaporates, but the 
molecule is formed by the combination of an 
impinging O, molecule and a WO group from the 
wire surface. Hence as far as the surface is 
concerned, it is W plus O which leaves. By 
analogy, therefore, we shall assume that Mo 
plus O is capable of leaving the wire surface, 
perhaps combining with an O2 molecule to form 
Mo0O;, although it does not matter for our con- 
siderations how it leaves. 

Let us take the rate at which Mo plus O leaves 
the wire to be directly proportional to the frac- 
tion of the surface covered by gas. Since we 
have previously assumed that Mo alone leaves 
at a rate proportional to the fraction of the 
surface not coated, it follows that the total 
number of Mo atoms leaving the surface per 
unit area and time after the field is applied will 
be given by 

Vir=A(1—0:)+BA, 
and 
Vir/ Vor=[A(1—61) + BA, ]/ 
[A (1— 60) ] 
=[A(1—6,:)+B0,]/A, 


or 
Vir/ (16) 


where A, B, and s are assumed to be constants 
depending upon the temperature of the filament, 
but not upon the electric field applied. Substi- 
tuting for @, from Eq. (12), we get 


bs'fi(2T exp (¢:/kT1’) 
Vir Ca (Ts) (G1 + RT 1’) 


exp 


As E gets larger, 6: approaches unity, and 
hence we see by Eq. (16) that Vir/Vor ap- 
proaches s as a limiting value at high fields, 
Thus Vir/Vor=s is an asymptote for our curve, 
If we are to apply Eq. (17) to Estabrook’s data, 
we will have to know s at temperatures of 1452°K 
and 1644°K. From the data at 1452°K, it seems 
to be about 0.01. Since the data at 1644°K do 
not extend to high enough fields, it is difficult to 
determine a value of s for that temperature. By 
trial, it was found that a curve which fit the data 
best would be obtained if s were about 0.12. 
This seems too large in comparison with the 
asymptote at 1452°K. The choice of 0.12 was 
needed because of the last point on the curve. 
Upon throwing out the latter, it becomes quite 
impossible to determine an asymptote from 
Estabrook’s data. By developing an approximate 
expression for s, and substituting into it the best 
numerical values we could find, we obtained 
s+=0.017 at 1644°K. This seems more reasonable, 
and so we shall use it. 

When we use the above values for s in Eq. (17) 
and apply the latter to the data, we again find 
that ¢; must vary approximately linearly with E 
in order to give agreement. In order to get com- 
pletely satisfactory correlation between theory 
and experiment, however, it is found to be 
desirable to have ¢; increase linearly with the 
fraction of the surface covered as well as with 
the electric field. The equation, therefore, is of 
the general form 


(17) 


dot (18) 
and if we substitute for 6; from Eq. (16), we get 
$1= do+gE+j(Vir/Vor—1)/(s—1). (19) 


When this relationship is put into Eq. (17), we 
obtain an expression in terms of Vir/Vor and E 
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which we shall use as the equation for the curves 
in Fig. 2. Values for the constants ¢o, g, and j can 
be determined by making the equation fit the 
data. If we choose p;’ = 10-7 mm of Hg, we get for 
$o, g, and j, respectively, the values (17.92)10-*° 
cal./atom, (2.86)10-°* cal. cm/atom volt, and 
(0.90) 10-*° cal./atom at 1452°K, and (19.42)10-*° 
cal./atom, (3.40)10-°* cal. cm/atom volt and 
(0.70)10-*° cal./atom at 1644°K. The curves 
shown in Fig. 2 were obtained by using these 
values in Eq. (19), and the latter in Eq. (17). 
They are seen to fit the data quite well except 
for the last point on the curve at 1644°K, which 
we decided to ignore, as mentioned before. 

The fact that it is necessary to include the 
term 76, in Eq. (18) is not at all surprising. The 
heats of desorption of gases are known to change 
markedly with the fraction of the surface coated. 
Roberts,” for example, found the heat of de- 
sorption of H from W wires to decrease from 45 
kilocalories per mole for a bare surface to 18 for 
a surface covered with gas. This change is pre- 
sumably due to interaction energy between the 
adsorbed atoms. 

The question then is not one of why a variation 
of ¢: with the fraction of the surface covered was 
found, but instead—why was the change in our 
case so small? The answer to this is found in a 
theoretical treatment by Miller™ of the problem 
of adsorption of dipoles. Quoting from his paper, 
we have “. . . the electrostatic and van der 
Waals forces give contributions to the variation 
of the heat of adsorption which are of opposite 
sign, and almost counterbalance one another, so 
that the resultant variation in the heat of ad- 
sorption is very much less than would be expected 
from consideration of forces of one type only.” 


CONSIDERATION OF THE TIME LAGS, 
RESISTANCE CHANGES, ETC. 

As has been mentioned before, the resistance 
and vaporization changes do not take place in- 
stantaneously. This is apparent because of the 
rounding of the corners in Fig. 1. On the basis 
of the theory proposed, the explanation for the 
time lags becomes quite simple. Upon applying 
an electric field, a gas layer begins to form on the 
filament. It is not until the fraction of the surface 


2 J. K. Roberts, Proc. Roy. Soc. A152, 445 (1935). 
#8 A. R. Miller, Proc. Camb. Phil. Soc. 36, 69 (1940), 


covered reaches a maximum, constant value that 
the resistance versus time curve can have a 
steady slope. In like manner, we must wait for 
the gas layer to leave the wire before we can ex- 
pect the slope of the curve to return to the value 
that it had before the field was applied. 

When we calculate the lags to be expected 
on the assumption that they are due to the time 
required for gas layers to form and leave, we 
find that the magnitudes observed can be ac- 
counted for quite easily. It must be assumed, 
however, that f; increases with the electric field 
if we are to explain the lags with exactly the 
same numerical values that were used in the 
treatment of the vaporization change. 

By comparing Figs. 2 and 3, we see that the 
resistance and vaporization changes apparently 
go hand in hand. Note, for example, that both 
effects begin at about the same field strength, 
and that if the curve in Fig. 3 were inverted, it 
would have the same general appearance as that 
in Fig. 2. In Worthing’s work at 2500°K he also 
found the two effects to begin at about the same 
field strength and to increase together. It would 
seem, therefore, that the theory which has been 
presented to explain the decrease in vaporization 
rates should also be able to account for the re- 
sistance changes. 

Any such explanation must be consistent with 
the following experimental facts. (a) Worthing 
and Estabrook found resistance changes in oppo- 
site directions. These findings were not contra- 
dictory, because different metals were used, and 
at different temperatures. It means, however, 
that any theory must be able to account for a 
change in resistance in either direction depending 
upon the conditions. (b) Vissat observed re- 
sistance changes at room temperature. Thus the 
changes can occur even when there is no appre- 
ciable vaporization. (c) Vissat found that his 
resistance changes were influenced by vacuum 
conditions, and were accompanied by a leakage 
current. (d) Lloyd Jones, working with a very 
good vacuum, did not observe any resistance or 
vaporization changes. 

Let us attribute the resistance changes to 
variations in the temperature of the wire re- 
sulting from its being coated with gas when a 
field is applied, and uncoated when field is 
withdrawn. If the temperature changes are pro- 
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duced by changes in emissivity, and if the 
emissivity in turn varies linearly with the frac- 
tion of the surface covered with gas, then it is 
not hard to show that the percentage change in 
resistance would be directly proportional to the 
change in the fraction of the wire surface covered 
with gas, or 


AR/R= (constant) (6; — 4) 
= (constant) (22) 


With this relationship, the curve in Fig. 3 would 
look the same as the 1452°K curve in Fig. 2 
turned upside down. This is seen to be approxi- 
mately the case. 

The fact that Worthing observed a decrease in 
resistance of his W filament when a field was 
applied at 2500°K, while Estabrook and Vissat, 
using Mo and W wires, found increases at lower 
temperatures can be understood if we imagine 
that there is a spectral variation in the behavior 
of the emissivity when a filament is coated with 
gas. The results are explainable if the presence 
of a gas layer decreases the long wave-length 
emissivity below that for a bare wire, but in- 
creases the short wave-length emissivity so that 
at high temperatures there is a resultant increase 
in the total emissivity, but at lower temperatures 
a resultant decrease. 

Another feature of this work that must be con- 
sidered is the fact that Estabrook obtained a 
vaporization rate for Mo at 1452°K which was 
very much greater than that predicted by an 
equation set forth by Jones, Langmuir, and 
Mackay.‘ Estabrook’s value was 1.4 g/cm? 
sec., while from J., L., and M., we get 8.0X10-" 
g/cm? sec. Because of this discrepancy, Norris! 
also studied the vaporization of Mo. His work 
gives 1.9X10-" g/cm? sec. as the vaporization 
rate for Mo at 1452°K, and hence does not agree 
with either of the above values. Thus the ques- 
tion of what the correct vaporization rate is 
remains unsettled. 


“LL. A. Norris, M.S. Thesis, University of Pittsburgh 
(1933). 
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It may be that the variation in values obtained 
was due to differences in extent to which the 
filaments were covered with gas. If such is the 
case, then Estabrook’s value should be the best 
because it is the highest. This is in agreement 
with our treatment above, in which we found his 
filament to be essentially uncoated before the 
field was applied (i.e., +0). The vaporization 
rate of Mo as given by Jones, Langmuir, and 
Mackay would then seem to be in error. As a 
partial corroboration of this, we note that 
Wahlen and Whitney" were forced to discard the 
vaporization rates of W as predicted by equations 
of J., L., and M. In so doing, they say “‘. . . but 
their (J., L., and M.) calculations are based on 
evaporation data obtained at a time when 
vacuum conditions had not been developed to 
the present extent; so this value is somewhat 
questionable.” 

At the beginning of this paper we mentioned 
that Lloyd Jones, working with a very good 
vacuum, had failed to observe either the resist- 
ance or the vaporization change. Such a fact 
would of itself lead one to suspect that the 
phenomena were in some way due to the gas in 
the tube. Since we have attributed the effects 
observed to an adsorbed gas layer, presumably 
oxygen, on the wire, we can easily account for 
Jones’ failure to find anything by assuming that 
he had removed too much of the gas from 
his tube. 
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The consideration of the phases of the fractional-parentage coefficients allows the extension 
of the matrix methods to configurations with more than two equivalent electrons. Tables are 
given for the parentages of the terms of ~* and d*. Applications are made to the spin-orbit 
interaction of the d* terms and to the electrostatic interaction between the configurations @", 
d"™—1s5, and d*~*s?, Errata in Part II are indicated. ™p 323 belo CM 


§1. INTRODUCTION 


HIS paper deals chiefly with the application 

of matrix methods to calculations within 
configurations with more than two equivalent 
electrons. 

It is known that the eigenfunctions built up 
with the usual vector-coupling formulas! are not 
antisymmetrical as required from the exclusion 
principle and they must be antisymmetrized 
afterwards. But if certain of the electrons are 
equivalent, these antisymmetrized states are no 
longer normalized and some of them are linearly 
dependent, so that the calculations become very 
complicated. 

An escape from these difficulties was proposed 
by Gray and Wills? who started from the nlm,m, 
scheme with antisymmetrized eigenfunctions and 
computed the SZ eigenfunctions using angular- 
momentum operators and orthogonality con- 
siderations. This method leads to an ortho- 
normal system of eigenfunctions, but since it 
gives up the vector-coupling formulas, the 
matrix of each operator must at first be calcu- 
lated in the nlm,m, scheme and then transformed 
to the SL scheme, and no use may be made of 
the powerful matrix methods developed in 
Chapter III of TAS' and also extended in a 
previous paper of the author.* 

In order to make full use of the above-men- 
tioned methods, we shall calculate the eigen- 
functions of the configuration /" as linear com- 


1E. U. Condon and G. H. Shortley, Theory of Atomic 
S; a (Cambridge, 1935) (which we shall denote by TAS), 


maniet Gray and L. A. Wills, Phys. Rev. 38, 248 (1931); 
*G. Racah, Phys. Rev. 62, 438 (1942) (which we shall 

denote by II). We refer to this paper and to TAS for 

definitions, notations, and bibliographical indications. 


binations of the eigenfunctions obtained by the 
addition of a further electron / to the con- 
figuration /*-'. This possibility was already 
indicated by Goudsmit and Bacher,‘ who intro- 
duced the concept of fractional parentage; but 
they were interested only in the squares of the 
coefficients of these linear combinations and 
calculated them with a procedure which, being 
based on a diagonal-sum method, did not permit 
them to separate the fractional parentages of 
duplicated terms.5 The consideration of the 
phases of the coefficients of fractional parentage 
will enable us to calculate them separately also 
for terms of the same kind occurring in a given 
configuration and to calculate the matrix ele- 
ments of every symmetrical operator between 
configurations containing equivalent electrons. 
The fractional parentages of the configurations 

p" and d"* are calculated in §3 and §4, whilst §2 © 
contains a lemma on which these calculations 
are based and §5 deals with the matrices of sym- 
metrical operators. In §6 an analysis is made of 
the structure of the configurations /" in connec- 
tion with the appearance of more terms of the 
same kind, and §7 contains an application to 
configuration interactions. 


§2. TRANSFORMATIONS BETWEEN THE DIFFER- 
ENT COUPLING SCHEMES OF THREE 
ANGULAR MOMENTA 


If we add two angular momenta j; and js, the 
magnitude J of the resulting vector and its 
s-component M characterize completely the 
states of the system ; but if we add three angular 
momenta, several states may occur with the 


as ene and R. F. Bacher, Phys. Rev. 46, 948 
- Don. Menzel and L. Goldberg, Astrophys. J. 84, 1 
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TABLE I. (p*SL[p*(S’L’)pSL). The different rows are 
normalized separately, and N is the normalization 
factor of each linear combination. 


N 1§ 3p iD 
«S$ 1 0 1 0 
op 18-4 2 —3 
2D 2-4 0 1 -1 


same J and M and a complete characterization 
of the states needs the specification of the type 
of coupling of the vectors. 

We may for instance couple at first j; and je 
and then add j; to their resultant J’: In this case 
the eigenfunctions are 
= V(jrjoJ’M’) (jams) 

m3M’ 
*(J’jsM'ms| J’jsJM) 
= L 


mimym3M’ 
jij2J'M’) 
(1) 
but we may also couple at first j2 and j; and then 


add their resultant J” to j;, and in this case the 
eigenfunctions are 


jejis(J”), JM) 
= 


(jojsmems| 
| jiJ”JM). (2) 


The unitary transformation which connects 
these two representations of the same system is 


(jrja(J") js | J) 
= 
mymom3M’ M’’ 
| 
| IM) ; (3) 


introducing the expression (16’)II for the trans- 
formation coefficients for the addition of two 
angular momenta and using Eqs. (19)II and 


(37)II we obtain ro 
| jojs(J”), J) 
=[(2J’+1)(2I" +1) js; JJ"), (4) 


where W is the function defined by (36’)II. 

It is sometimes useful to consider the changing 
of the coupling together with a change in the 
order of the vectors; the same way as before 
yields 


| 
=[(2I'+1)(2I" +1) ; (5) 


If we have three electrons or groups of elec- 
trons, the transformations between the different 
parentages in SZ coupling are obvious extensions 
of (4) and (5): For instance, 


| sili, SL) 
+1) 
(6) 


a particular case of this transformation was con- 
sidered in TAS 6° 14. 


§3. THE EIGENFUNCTIONS OF GROUPS OF 
EQUIVALENT ELECTRONS 


If we couple two equivalent electrons with the 
usual vector-coupling formulas, we obtain anti- 
symmetric or symmetric eigenfunctions accord- 
ing to whether S+L is even or odd (TAS, p. 231) ; 
the eigenfunctions of the states with S+L even 
are therefore the normalized eigenfunctions of 
the allowed states of P. 

If we add in the same way to the allowed states 
of F a third / electron, the obtained eigenfunc- 
tions are in general antisymmetric only with 


respect to the first two electrons, but not with’ 


TABLE II. 


N oS 2P 1D 1G 
30-4 0 74 154 —8 0 
15-4 0 — 8 0 —7 0 
2D 60-+ 4 —3 —3 
#D 140-4 0 454 21% 
°F 0 284 — 104 7% 
5-4 0 0 0 
#G 42-4 0 0 — 104 21% 
24 0 0 0 1 


h 


6) 


| © 


respect to the third. If we apply in effect to 
y(2(S’L’)ISL) the transformation® 


= W(S"L"), SL) 
(1, (7) 


where the transformation matrix is given by (6), 

we obtain in general in the sum (7) allowed and 

forbidden values of S’’L” and, therefore, 

y(?2(S’L’)ISL) cannot be an eigenfunction of F. 
Only such a linear combination 


¥(aSL) = 
s’ 
(8) 


may be the eigenfunction of ? for which the 
coefficients of (I, SL) vanish for every 
forbidden value of S’L” after the application of 
the transformation (7) ; the ‘coefficients of frac- 
tional parentage” must 
therefore satisfy the equation system 


U(S"L"), SL|P(S'L’)ISL) 


(S"+L" odd). (9) 


Since a function antisymmetric with respect 
to the electrons 1 and 2 and also with respect to 
the electrons 2 and 3 is antisymmetric with 
respect to all three electrons, the condition (9) 
is necessary and sufficient for the determination 
of the coefficients of fractional parentage of the 
terms of /*, and the number of independent non- 
vanishing solutions of (9) for a given SL equals 
the number of allowed terms of this kind in FP; 
if this number is greater than one, the different 
terms may be distinguished by a parameter a. 

As an illustration of this method let us cal- 
culate the eigenfunction of the term p*?D. It 
follows from (6) that 


2D) = (3/16) W(p, pp(?D), *D) 
+(3/4)¥(p, pp(*D), 2D) 
—(1/4)¥(p, pp(*P), *D), 


Since we shall mostly consider transformations and 
operators which are diagonal with respect to Ms and M, 
and independent of them, we shall in general neglect these 
quantum numbers. 
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TABLE III. 


oS 1 0 o 1 0 0 0 0 0 
as 1 0 o 0 1 0 0 0 0 
2P 3604 1354 -—35' 0 0 
5 10 4 0 0 
2D —424 0 1054 454 o 0 
140% 424 0 20 0 154 0 
207 7 +o oo 45t 0 
10% 0 3 0 7 0 0 
560% 1204 0 0 200! —105! o 
1680-4 | —200! 0 —1754 —1124 —<05! 
20 504-4 0 1898 —5 70! 0 —154! 
eG 168074 0 0 o 200! 3154 297! 30st 
eH 0 0 2 264 
10+ 0 o oO 0 0 0 7 


¥(p?('D) p *D) = (3/16) pp(*D), *D) 
— pp('P), *D) 
—(1/4)¥(, pp(*D), *D) 
+(3/4)¥(p, pp(?P), *D); 
since in the development of 
V(p* *D) (p?(*P)p *D) *D) 
the coefficients of 


pp@D),*D) and ¥(p, pp(P),? D) 


must vanish, the only possibility, apart from a 
phase factor, is 


*D) = (1/2)'y(p?(*P)p *D) 
— (1/2) *D). 
The same method may also be extended to the 


configurations /", if the fractional parentages of 
I’ are known. In this case 


V(l"aSL)= 


a’S'L’ 
S'L')ISLW"a SL) 
a’ 
SY S'L’) 
(10) 


and the coefficients of fractional parentage 
(I—"(a'S’L')ISLW"aSL) must satisfy the equa- 


ing 
the ‘ 
ore 
ec- 
ent 

he 
ti- 
; 
en 
of 
es 
th’ 


370 GIULIO RACAH 
TABLE IV. 
ds 

JeS ¢S 27D #D “PF oF eG eH “ 
ss s4io 0 0 0 o 33 0 0 0 0 0 0 0 0 
&S 1 |0o 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 
1504/0 0 144 5 30! 10! 0 -ist -5 0 0 0 0 0 
*#P 3007/0 o -8 144 0 0 354 o 56} 0 0 0 0 0 

50-4 | 0 .-3 0 —st 0 0 0 o -3 0 0 0 
iD 3504/0 -144 -7 60+ 0 354 4st 0 
#D 7007/0 —sot 0 1264 0 908 60! 0 0 13894 0 994 4s} 0 0 
“D 7004/0 0 0 1264 0 o 0 o 0 0 180+ 0 0 
#F 2800-4/ 0 -—200! -160! 180! 1208 0 ~-175! -20 220! 0 
28004) 0 0 0 3608 0 -10 #600! 0 o —315¢ 0 49st —3 —396! 0 
*F 7004/0 -4 0 o —175% 2248 9 0 —110! 0 
2G 84007+/ 0 0 0 0 -10 6008 0 —7+ 1680! 945! sso 84s? soit 924th 7288 
8G 184804] 0 0 oO 14524 968! —2541! 0 4235! -1215' —ss774 -—308! —2184) 
4204/0 0 0 0 0 5 0 —708 0 0 —66' 1544 0 
2H 11007] 0 0 0 0 0 0 0 334 —220! 2208 —st 286% 1824 
0 0 0 0 0 0 0 0 0 0 0 2314 
tion system not exhaust all states of /"—'l, but only those 


(S"L", MSL"), 
S'L)ISLW"aSL) =0 
odd). (11) 


The systems (9) and (11) do not fix the phases 
of the eigenfunctions of the different terms, nor 
the scheme in the case of more terms of the same 
kind, but give the fractional parentages in any 
arbitrary orthonormal scheme; the convenience 
of a particular choice of the scheme will be con- 
sidered in §6. 

The fractional parentages of the terms of #’, 
d*, and d® calculated with this method are 
given in Tables I-IV. The phases of the eigen- 
' functions of p* and d* are in agreement with 
those of TAS 4° 6j and 5°6 with the exception 
of p* ?P; it must however be remarked that these 
phases differ from those of Ufford’ for the terms 
‘P, *F, *F, and °G of a’. 

The coefficients of fractional parentage con- 
sidered by Goudsmit and Bacher and by Menzel 
and Goldberg are nm times the squares of our 
coefficients. 

It must be pointed out that the matrix 
is not an ordinary uni- 
tary matrix, but only a rectangular matrix which 
is a part of a unitary one, since its columns do 


7C. W. Ufford, Phys. Rev. 44, 732 (1933). 


which are allowed in /"; the hermitian conjugate 


(Pa 


(12) 
does therefore satisfy the relation 
2 S'L" 


but a matrix multiplication in the opposite order 
has no sense, if the sum is limited to the antisym- 
metric states of /". In calculations with only sym- 
metrical operators we may however, treat 
formally the matrix (/"—(a’S’L’)ISLW"aSL) as 
a common unitary matrix without weakening the 
general laws of matrix calculations, since sym- 
metrical operators do not connect states of 
different symmetry and, therefore, the sum over 
the neglected states vanishes. 


§4. FRACTIONAL PARENTAGES IN ALMOST 
CLOSED SHELLS 


We shall determine in this section a relation 
between the fractional parentages of the terms 
of an almost closed shell /*‘+?-" and those of the 
terms of /**!, This relation will not only avoid 
long numerical calculations, but will also give 
us the eigenfunctions of the terms of /*'+?-* with 
the phases fixed by the convention of §6 of II. | 


1 
| | 
{ 


I 


8 


2) 
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According to §6 of II two terms of /* and of and 


]##+#2-* will be called conjugated® if their eigen- 
functions appear multiplied with each other in 
the relation 


aSLMsML 


[(2S+1) 


*(2L+1) e(l"aSLM 5M) 

SL —Ms— M1) 

*(SSMs—Ms|SS00) 
(14) 


where £ denotes the group of the first » electrons 
of the shell and ® the group of the remaining 
4l+2—n; owing to (16’)II we have 


1S) = ( 


aSlMsM L 
(15a) 


In the same way, if we consider the group { 
of the first +1 electrons of the shell and the 
group ®’ of the remaining 4/+1—~n, we may also 


write 
15) = 
n+1 
a’ 
1’) 


(15b) 
It follows from (10) that 
Ve s'M 1’) 


Ve(l"aSLM sM 1) 


aSLMsM 


sm, | S$S'M 5’) 
(LIM Lm} | LIL'’M 
(16) 


® To the correlation defined in §6 of II we shall reserve 
the word “conjugation,” since other t 
between terms 
§6 of this paper. 


of correlations 
different configurations will be found in 


a’ 
— M s'— Mx’) 
Ms'm,|S'3S— Ms) 

here is the group of the electrons 
n+2, ---4l+1. Since Yx is antisymmetric, the 
substitution of the electron 41+-2 by the electron 
n+1 and of the group R” by the group 8’ mul- 
plies ¥x by (—1)**', and then 
Vx M1) 

a’S'L'’Ms'’ ML’ memi 
S'L'’ — Ms'— Mz’) 
M3'm,|S'4S— Ms) 
-(L'l— M'm,| L'IL— Mz) 


If we introduce (17) in (15a) and (16) in (15b), 
we may equate the coefficients of each product 
WV ¢n41¥m' separately, since for different quan- 
tum numbers these products are orthogonal, and 
obtain 


41+2\-! 
-(S’'$— Ms'm,|S'$S— Ms) 
-(L'l— 'm,| L'IL— Mz) 
SL) 
Ca 
-(S$M gm, | S$S'Ms')(L1M im, | LIL'’M 1’) 
(18) 
Owing to (16’) IT and (19a) II, and to the fact that 
n+1 has the same parity as 2(S’+L’), we get 
(a! SL) 
= 
(n+1)(2S'+1)(2L'+1) 7 
—n)(2S-+1)(2L+1) 
(19) 
which is the requested relation. 
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TABLE V. (p*SL|| U®||p*S’L’). entages of the terms of 7 it must, however, 
be observed that (19) gives the parentages of 
ad with respect to and that the 

eigenfunctions of /?'+! determined by the methods 
“ap 36 0 14 (21)! of the preceding section are Ve(/?'*'); since it 
was shown in §6 of II that the terms of ?'**' split 

TABLE VI. (p*SL||U®||p*SL’). in two classes, according to the two possibilities 

of (76)II, we must change the sign in the relation 

‘S 2p (19) if belongs to the class for 

4S 0 0 0 which the minus sign holds in (76)II. The clas- 

= : a —@) sification of the terms from this point of view 
(3) will be considered in subsection (5) of §6. 

From (19) and (13) we obtain also le: le 
~ (2S+1)(2L+1)(-'(a'S’L')ISLW"aSL) We are at first interested in the matrix com- 
a I II i 
ponents (A!| F| A!) of quantity 

(20) F=2i fi (21) 


For the determination of the fractional par- where f; is an operator which operates on the 


TABLE VII. 


#P #D vF 7G #H 
#P —2(21)8 0 —216(10)! 16(210)! —4(21) 0 0 0 
0 0 ‘ 0 0 — 14(6)3 0 0 
216(10) 0 356 —7(15)+ 0 —15(7)* 0 
¥D — (210) 0 1 0 —5(3)+ 0 
—4(21) 0 (15)3 9(35)+ 7(6) 0 2(210) — (2310)! 
0 — 14(6)! 0 0 0 —7(6)* 0 0 
0 —15(7)! —5(3)% —2(210) 0 3(22)! — (462)! 
vH 0 0 0 0 — (2310)+ 0 (462)! (3003) 
TABLE VIIa. (d* 
oS as 2D aD 1G “G 
oS 0 0 7(30) 0 0 0 0 0 
aS 0 0 3(70)4 4(35)+ 0 0 0 0 
2D 7(30)+ 3(70)3 -5 — 30(2) 0 4(5)+ 8(55)* 0 
aD 0 4(35)+ — 30(2)! —15 10(14)8 10(10)+ 2(110)* 0 
aF 0 0 0 —10(14)+ —7(70)4 — 0 
1G 0 0 4(5)4 10(10)+ 5(22)! 5(2)4 —2(455)! 
2G 0 0 8(55)4 2(110) 14(770)+ 5(2)+ —12564(22)4 — & (5005) 
al 0 0 0 0 0 —2(455)+ — ,(5005)* 35, (143)! 


TaBLeE VIIIb. (dvSL||35 U® |\dv’S’L’) for S=1, 2. 


#D #D oF @H 
iP 144 (6)! — 284(15)3 0 144 (6)? 284 (6)+ 0 0 
28¢(15) 44(210)! 5 0 —4(35)+ 4(35) 20(3) 0 
SP 435)! 49606)! 3(210)! 
oF 284 84(21)3 —4(35)4 0 496 (6)! 7% — 14(210)! 44(2310)8 
&G 0 0 20(3) 0 3(210)+ 46(210)4 — 36(22) — 2(462) 
@H 0 0 0 0 2¢(2310)# (2310)! 2(462)$ 14 (3003)! 


372 
= 
* 
= 
3} 
1] 
2] 
2] 
fi 
0 
tl 
( 
_ 


THEORY OF COMPLEX SPECTRA wy , 373 
gS 0 0 0 70) 0 0 0 0 0 0 
0 —7(30)4 0 5(105)+ 0 4(105)+ 0 0 0 
O 7(30)4 0 0 —14(5)4 0 —10(21)4 0 0 0 
#D 4(70)3 0 15(7)* 0. 5(2)8 0 10(7)4 0 —6(55)# 0 0 
2D 0 —5(105)4 0 5(2)8 0 0 0 20(6)* 0 0 0 
0 0 14(5)4 0 0 —7(30)4 0 0 0 
0 4(105)* 0 —10(7)4 0 —7(30)4 0 4(42)4 0 = —2(462)' 0 
0 —10(21)4 0 (6)+ 0 —4(42)4 0 9(30)* 0 4(273)4 
2G 0 0 —6(55)* 0 0 0 9(30)# 6(70)* 0 
0 0 0 0 —2(462)4 0 —6(70)4 0 7(13)¢ 
|) 0 0 0 0 0 0 4(273)% 0 —7(13)! 


TABLE IXb. ,AL’). 


0 7(15)8 0 
—7(15)4 8(35)! 0 
0 —8(35)4 0 15(14)4 
0 —15(14)+ 0 
TABLE X. (p*SL\|64V || p2S’L’). 
1s ap ip 
0 (30)# 0 
—14(30)4 (6)4 
1D | ~ 0 
TABLE XI. 
4S 2P 2D 
‘Ss 0 2(3)8 0 
2p 2(3)8 0 (15)4 
0 —(15)4 0 


electron 7, and \! and X" are states of the con- 
figurations I and II; owing to the antisymmetry 
of ¥(A') and W(A!) we have 


f (22) - 


If putting i= and assuming for ¥ 
the expression (10), we obtain 
("aSLMsM F | Ina’ S'L'’M s'M 1’) 
=n > 
aiSili 


,SLM sM fa| S’L'’Ms'M 11’) 
("ay (23) 


where (SiLi,SLM sM1 | | SiL J, 
may now be calculated with the ordinary matrix 
methods of Chapter III of TAS and of II. 

As application of this formula we calculated 
the matrix components of the tensors U® and 
VD, defined by (102)II, for the configurations 
p*, p®, and d°; the results are given in 
Tables V—XIV. For d? the matrices were already 
given by (103)II; it must, however, be noted 
that an error occurred in the final form of the 
manuscript, and all the elements of (103c)II and 
(103d)II must be multiplied by (3/2)!. 

From the elements of V“ the matrix com- 
ponents of the spin-orbit interaction may easily 
be obtained : it follows in effect from the relation 


t= (1+-1)(2/+1) (24) 
and from (38)II and (102)II that 
(IPaSLJIM | ; 
= }! 
|| VO» J1). (25) 


If in (22) [=/" and II=/"~/’, the terms of II 
are characterized by S and L and by the quantum 
numbers of the parent ion /"-'; W(A") has in 
this case the expression 


s'M 1’) 
= (1/n)' (—1)" 
(26) 


where in the right side we consider the group 
as composed by the electrons 1, 2, ---, 7—1, 
i+1, ---, m, and P; is the parity of the per- 
mutation which exchanges 7 with n. Introducing 


\ 
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TaBLe XII. ||dv'S’L’). 


#P #D “F #G #H 

#P 2 —2(14)# 9 (2)4 0 0 0 0 

. 2(14)4 (10)4 —4(3) 0 0 0 0 0 

2D —4(3)4 34(5)4 — 34(105)+ 42)4 — (42)4 0 0 

— 0 — 34(105)4 — 2) 5(2) 0 0 

0 0 — —(2)4 —14(14)4 —(14)# 3g—(10)4 0 

0 0 — 5(2)4 (14)4 2(35)4 — 3(10)4 0 
0 0 0 0 3¢(10)4 —3(10)4 (55)4 
0 0 0 0 0 —9§(S5)# (55)4 


TABLE XIII. 


oS 0 0 3(3)3 0 0 0 0 0 0 0 0 i) 0 0 0 0 
os 0 o --( 2(2)4 0 0 0 0 0 0 0 0 0 0 0 0 
2P 3(3)t 1 0 0 0 0 0 0 0 0 0 
0 2 2 (2)4 0 0 0 0 0 0 0 0 
2D 0 0 —ila)) 2 0 0 2(10)4 0 0 2 4 0 0 0 0 0 
“D 0 0 a7 0 o 0 0 (2) 0 0 0 0 0 
0 0 0 0 154 0 (7 0 0 0 0 
oF 0 0 0 0 0 0 —2(5) 0 0 (35) 3g(35) -% 0 0 
oF 0 0 0 0 a (70) —(14)) —3g(14)# — 34083)! (10)! 0 0 
2G 0 0 0 0 0 0 0 0 o -@ 53) 0 0 3 —(66)* 0 
0 0 0 0 0 0 0 (33) —34(33)h 0 0 
0 0 0 0 0 0 0 -% —3(10)' —3 —34(11)# —*/6(55)4 0 
eH 0 0 0 0 0 0 0 0 0 0 0 —3¢(26) 
“0 0 0 0 0 0 0 0 i) 0 0 0 0 0 —39(26)% 0 


| 


(10) and (26) in (22) and putting i=, we have 
(MaSLM gM s'Mz') 
=n} 
5M 1’). (27) 


This formula, which is the extension of TAS 6°17, 
gives a rigorous demonstration to the method of 
Menzel and Goldberg® and also fixes for such 
transitions the phases of the matrix components, 
which are necessary for transformations to other 
types of vector coupling (TAS, p. 252). 

If in (22) and the 
terms of each configuration are characterized by 
S and L and by the quantum numbers of the 
groups of equivalent electrons; in the same way 
as for the precedent case we obtain 


l’?(a2S2L2), 
SLMsM F | 


*D. H. Menzel and L. Goldberg, Phys. Rev. 47, 424 
(1935) and reference 5. 


3 
S2Le, 
SLM 5M1|fa—p| 
Sole, S'L'’M s'M 1’) 
Sole, S’L’| Sy'Ly', S2L2(S2'L2’), S’L’) 
In connection with this result it must be 
observed that (I, SLW"aSL) is not 
(In (a’S’L’)ISLW"aSL), but it is easy to see 
that the two coefficients are connected by the 
relation 


(1, SLW"aSL) 
= 
S’L')ISLW"aSL). (29) 


It is unnecessary to consider the matrix ele- 
ments of F for transitions between more com- 
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TABLE XIV. (d5vSL||30! 


| 


ws 0 o -4 (4p 0 0 0 0 0 0 0 0 0 0 0 0 
es 0 0 3(10)+ 0 0 0 0 0 0 0 0 0 0 0 
4 0 0 0 1 2(2)8 0 0 0 0 0 0 0 
o -8 0 — (70) 0 0 0 0 0 0 0 0 
2D 0 0 —8 0 —(35 0 0 0 0 0 0 0 0 
#D 0 0 0 0 (10)# —4(5) o 0 0 0 0 0 0 
#D 0 o -1 —2(14)+ 0 aot 0 8 o -2 0 0 0 0 0 
0 0 (70 0 0 0 —4(5)) 0 0 0 0 0 
0 0 0 o Oo -8 0 0 —34(66)! 5(6)4 0 0 
oF 0 0 0 0 200% 0 —35(70) 0 —(70)' 0 0 0 0 
0 0 0 o +o -2 4(5)4 0 (7 0 0 
«(0 0 0 0 0 0 0 0 0 —3¢(2) 0 0 —3¢(22)! —3(2)t 0 0 
0 0 0 0 0 0 0 44 (66)? 0 —(66)# 0 o 0 
0 0 0 0 0 0 0 0 0 0 0 0 0 
vl (0 0 0 0 0 0 0 0 0 0 0 0 3(13) 0 


plicated configurations, since all other cases may 

be reduced to these three by means of TAS I*16. 
The calculation by the same method of the 

matrix components of the scalar operator 


(30) 


i<j 


needs in some cases the knowledge of V(/"aSL) 
as linear combination of the eigenfunctions of 
PP: 


W(maSL)= 


a1S1L1S2L2 


P(S2L2), 
P(S2L2), SLW"aSL); (31) 


the coefficients of this expression are given by the 
formula 
?(S2L2), 
= P(S2L2), 


(32) 
The following results are easily derived: 
(InaSL|G| SL) 


a1S1L1S2L2 


SL) g|?S2L2) 


P(S2L2), SLY"a’SL), (33a) 


SL) 
(n—1)n* 


a1S81L1S2L2 


B(S2L2), SL) (bd S2L2) 
’(S2L2), SL| SL) 
l’*(S2L2), SL) 
P(S2L2), (33c) 


Since the actual application of (32) needs 
generally very long calculations, the formulas 
(33) are of practical use only in a few particular 
cases; in other cases it is simpler to express g;; 
as a sum of scalar products of tensors and to 
reduce the problem to the calculation of tensors 
of the type F. Applications of both methods will 


be shown in the next sections. 


§6. THE STRUCTURE OF THE CON- 
FIGUATIONS 


In this section we shall classify the terms of 
the configuration /" according to the eigenvalues 


of 
Q= gis, (34) 
i<j 


where q;; is a scalar operator which operates on 
the two equivalent electrons i and j and is 
defined by the relation © 


(PLM 0). (35) 


H 
i) 
) 
) 
) 
0 
0 
0 
0 
0 
0 
) = 
) 
) 
3) 
e 
| 
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It will be shown that to every term of /" with 
non-vanishing Q a term of the same kind cor- 
responds in /"-*, and this fact will allow us to 
assign to each term a “seniority number”’ ac- 
cording to the value of m for which the term 


- appeared for the first time. Some useful relation 


between the fractional parentages of correspond- 
ing terms will be obtained and it will also be 
shown that the classification of the terms of 
P+ according to the two possibilities of (76)II 
depends only on the seniority of the term. 


(1) The Eigenvalues of Q 
It follows from (42)II and (40a)II that 


(2r+1) Or) W (lll; Lr) =6(L,0); (36) 


0 
expressing W/(Illl; Or) by (36’)II and using also 
(38)II and (58)II we get for g;; the expression 


0 
Since u is a scalar and, owing to (33)II, 


(38) 
we have also 
(37) 
and 
(39) 
We shall henceforth consider only schemes for 
which also Q is. diagonal, i.e., schemes for which 
SL) =Q(l"aSL)i(aa’). (40) 


It follows from (39) and (74)II that if Q is 
diagonal in a given scheme of /”, it is also diagonal 
in the conjugate scheme of /‘'*?-", and that 


SL) = Q(l"aSL)+2i+1—n. (41) 


In order to calculate the possible values of 
Q(l"aSL), we express it by means of (33a) and 
(35): 

Q(l"aSL)5(aa’) = 
24S), SLW"a’SL). (42) 
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Multiplying the two sides by 
SL{l"-*(BSL), SL) 
and adding with respect to a’ we have 
Q(l"aSL) (4S), SL) 
=> l?(4S), SL) 
(4S), 
SL{I"-*(BSL), (4S), SL); 


the summation with respect to a’ may be made 
by means of (42) after transforming the last two 
factors by the relation 


(in—*(BSL), (4S), SLY"aSL) 
(41+3—n)(41+4—n)}! 
" n(n—1) 
SL), 1?(4S), (43) 


which is analogous to (19) and may be obtained 
in the same way; we get 


Q(*aSL) SL) 
= P'S), 


Owing to (41) (/"aSL{ln-*(@SL), P('S), SL) 
may be different from zero only if 


Q(l"aSL) (45) 


and according to (42) the only non-vanishing 
values of Q(/"aSL) are those which are con- 
nected to a Q(/"*BSL) by (45). 


(2) The “Seniority Number” 
Putting for 
vap(QSL) = SL), (46) 


where a@ may assume all values for which 
Q(l"aSL) =Q, and £ all values which satisfy (45), 
we have from (42) that 


3n(n —1)(2/+-1) 

Sp (47) 
and also from (43) and again (42) that 
}n(n—1)(21-+1) 

=05(88"); (47’) 
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it follows that for two given values of Q which 
are connected by (45) the number of independent 
states of given S and L is the same in /* and /"~, 
and that the matrix 


uap(QSL) (48) 


is a unitary one. If we apply to the eigenfunctions 
of /* the transformation u, and consider the 
states with eigenfunctions 


= Ya V("aSL) tag, 
we obtain 
P(1S), SLY"BSL) 
=(Q("BSL) — 1) (22+ 1) (49) 


i.e., it is possible to find a scheme of /” in which 
not only Q is diagonal, but also each term of /* 
with Q0 corresponds to a well-defined term of 
In-? whose Q is connected to Q(/"8SL) by (45). 
If also Q(/"*BSL) #0, this term corresponds to 
a term of /*-* and so forth; each chain of cor- 
responding terms begins with a term /*8SL which 
has Q=0. 

We may thus assign to each term in the QSL 
scheme a “seniority number”’ v, which indicates 
the number of electrons of the first member of 
its chain; it follows immediately from (45) that 
Q depends only on m and v and that its values 
are given by 


O(n, v)=}(n—v)(4l+4—n—v). (50) 


Confronting (41) and (50) we see that conjugate 
terms have the same seniority. 

The seniority number suffices for distinguish- 
ing the different terms of the same kind in the 
configurations d" but not in f*, since there are 
in f* terms of the same kind which have also the 
same seniority. For such configurations an un- 
specified parameter a must be maintained besides 
v; terms corresponding according to (49) will 
have the same values of v and of a. 

With this convention Eq. (49) which defines 
the correspondence between terms of the same 
chain becomes 


(I"-*(a’v'SL), 12(!S), SLYnovSL) 
= (O(n, }'L4n(n—1) 
5(aa’). 


In this paper all tables of matrix elements are 
given in the QSL scheme and the seniority 
number is indicated by a prefix under the mul- 
tiplicity number of each term: for instance the 
two 2D terms of d*, which were indicated in TAS 
(p. 228) by a?D and &D, are therefore, respec- 
tively, denoted by D and ;*D. 


(3) The High Degeneracies 


Majorana’s operator of position exchange may 
be defined for equivalent electrons by the relation 


(PLM (51) 


and may also be expressed, according to Dirac’s 
vector model, by 


Mi;= (52) 
From (43)II we have 
(—1)"(2r+1) Or) ; Lr) 
=(—1)4W(lll; LO); (53) 


expressing Or) and W(llll; LO) by (36)II 
and using also (38)II and (58)II we get for M,,; 
the expression 

21 


M;;= (54) 
0 


Adding this equation to (37) and introducing 
(52) we have 


-u,29) 
(55) 


It follows therefore from §4 of II and par- 
ticularly from (51)II that if Slater’s integrals F* 
are proportional to (2k+1)/Cin, the electro- 
static interaction between two equivalent elec- 
trons is proportional to (55) and then the 
electrostatic-energy matrix is diagonal in the 
QSL scheme and its eigenvalues are only func- 
tions of n, v, and S. This fact explains the high 
degeneracies observed by Laporte and Platt'® for 
these particular ratios of the parameters; unfor- 
tunately these ratios are only hypothetical, 
since they are excluded by the property of F* of 
being a decreasing function of k (TAS, p. 177). 


1°. Laporte and J. R. Platt, Phys. Rev. 61, 305 (1942). 
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(4) Relations Between Parentages of 
Corresponding Terms 
If we express ¥(/") as linear combination of 
with as linear combination of 
express on the other hand W(J") as 
combination of with as com- 
bination of y(/"-*l), and compare the two 
developments, we obtain 
12(4S), 
= 


SL), 12S), SL"ewSL), (56) 
and owing to (49’) and (50) 


* 


SL) 
= SL) 
(57) 


It is easy to deduce from this recursion formula 
that 


=0 

(v'¥v+1), (58a) 

(58b) 

=[(m—v)(v+2)/2n]}! 

(58c) 


Comparing (58b) with (13) we obtain the more 
accurate orthogonality relations 


a’S’L’ 

]8(aa’) (59a) 
and 

a’ 
+ 
=[(n—v)(41+4—v)/ 
2n(21+2—) ]8(aa”); (59b) 


comparing (58c) with (20) we obtain also 
(2S+1)(2L+1) 


aSL 
= (2S’+1)(2L’+1) 
(60a) 
and 


(2S+1)(2L+1) 
aSL 


—1S'L')ISL) 
= 
*(41+-5 —v)/2n(2/+3—v) (60b) 
Another useful relation is the following : 
= 
(2.S’+1)(2L’+1)(v+1) ] 
(2S+1)(2ZL+1)(v+2)(2/+1—v) 
(61) 
since this relation is verified for v=0 (S=L=0, 
S’ =}, L’=1), it suffices to prove that if it holds 
for v=v'—1 it holds also for v=v’. 
We use for this purpose the expressions (32) 
and (49’) of (l*’(av’SL), P(4S), 
owing to (6), (58), and (50) we have 


(- 


(2S:+1)(2Z:+1) 7! 
1)(2S+1)(2L+1 | 
(2l+1—v')v’ 
(v'+2)(21-+2 on 
+ 
a2SoL2 
(2.S2+1)(2L2+1) 
1)(2S+1)(2L+ 
+ SL) 
2(21+1—v’) 


= | . (62) 
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If (61) holds for v=v’—1 we may calculate the 
first sum with the aid of (59a) and obtain 


v’ 1 2i+1—v’ 4 


the second sum must then have the value 
21+1—v’ 
U+2—v'L2(v’ +1) +2) 


and owing to (59b), to (60b) and to the well- 
known corollary of Schwarz’s inequality, this 
fact is possible only if (61) holds also for v=v’. 

By the use of the formulas (58) and (61) the 
calculation of the fractional parentages is con- 
siderably simplified: Only the parentages of the 
“new” terms (v=) must really be calculated by 
the methods of §3; all others may be quickly 
deduced from them. 


(5) Relations Between Correspondence 
and Conjugation 


It follows from (43) that if two eigenfunctions 
Ve(l"avSL) and e(l"**avSL) correspond accord- 
ing to (49’), also the eigenfunctions of their con- 
jugate states and 
correspond in the same way. But if, in order to 
make full use of (74)II, we assume as standard 
scheme the scheme of the We for n& 2/+1 and 
that of the Yx for m 2 2/+2, we cannot use (49’) 
for the determination of (P(a’vSL), P('S), 
SL}P?'**avSL), nor can we use (19) for n=2/ 
without knowing which of the two possibilities 
of (76)II holds for each term of /?'+'. 

In order to solve these questions we consider 
provisorily the system of functions WV e(/?'**avSL) 
defined by means of (49’) and the system of 
functions Vx(/?'+!a’v'S’L’) defined by means of 
(14), and seek the relation between the parentages 
of with respect to 
and the parentages of ¥x(/?'+2avSL) with respect 
to Vx(l?'+e’y'S’L’). Using (19), (58), and (61), 
and owing to the fact that 2S is even and 25’ is 
odd, we get 


= (24 (63a) 


and 
(21 x 
= — (63b) 
If we assume 
YL), (64) 


it follows by the alternate use of (63a) and (63b) 
that 


Wx = (65) 


If we had assumed a minus sign in (64), the 
relations (65) would have also the opposite sign ; 
the choice between these two possibilities depends 
on the phase of ¥(/*'*?1,S), and it may be shown 
that (64) is in agreement with the convention 
of §5 of TAS for the eigenfunctions of closed 
shells. 

The relations (65) must be taken in account if 
we use (19) for n=21 or (49’) for n=2/+2 and 
n= 21+-3, and also if we calculate the coefficients 
of fractional parentage for n22/+2 by means 
of (58) instead of (19). 


(6) Relations Between Matrix Components 
of Tensors 


It follows immediately from (23) and (58) that 
the matrix components of every operator F 
between two states of /" may be different from 
zero only if 


and that 
(mavSLM gM F\l"a'v—2S'L'M s'M 1’) 
= 
-(ltavSLM gM | (67) 


Av=0, +2, (66) 


owing to (65) a minus sign, however, must be 
introduced in this formula for m 2 2/+2. 

If Av=0 the sum (23) splits in two sums ac- 
cording to the two possibilities y—1 and v+1 of 
the seniority numbers of /*-, but only the first 
sum may be immediately expressed as in the 
preceding case by means of the matrix com- 
ponents for /*; the second sum is to be expressed 
by means of the matrix components for /* and 
those for another arbitrary configuration /*’. If 
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we assume n’=4/-+2—z, the final result is found 
to be 
(navSLM gM _| F\la’vS’L'M s'M 
] 
(l*avSLM gM F\l*a'vS'L'M s'M 
+[(n—v)/2(2/+-1—») ] 
(68) 
If F is an irreducible tensor, it follows from 
(68) and (74)II that 
T™ 
= (l*avSL||T™ (x+kodd), (69a) 
(lmavSL||T 


2i+1-—n 
= 
2/+1-—v 


(«+k even). (69b) 
From (67), (65), and (74)II we get also 
(ew SL|| —28'L’) =0 

(x-+k odd). (70) 


The remarkable result that a tensor of odd 
degree is diagonal with respect to v and that its 
submatrices are independent of m may be ob- 
tained also in a more direct way. It follows from 
i the triangular conditions and from the fact that 
in 2 only states with even S+L are allowed, that 
for x+k odd 


(? 1$| + |2SLMsM1z) 
= (2SLM sM_| +4,() 1§)=0 


TABLE XV. (p°SL||2V“|| p°S’L’). 


1S ap 1D 
15 0 0 0 
| 1D 0 me 0 
| ‘s 
4S 0 0 —2(2)4 
2p 0 (6)3 0 
2D 2(2)4 0 —(14)4 
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and, therefore, 
+t) = +4) 

(k+k odd); (71) 
since all other ¢,“*) commute with q;;, * 
commutes with g;; and also with Q, and is there- 


fore diagonal with respect to v. From (71) we 
have also 


QOT=TQ= bigis: 
i<j 
calculating the matrix of this operator with the 
methods of §5 and owing to (49’) we obtain 


(ltavSL||T 
= 
which is equivalent to (69a). 

The matrices of the tensor V“®) defined by 
(102)II were calculated for the configurations p°, 
p', d‘ and using also (69a) ; the results are 
given in Tables XV—XIX. The matrices given 
in Tables V and XIX are sufficient for the cal- 


culation of the spectra of the configurations p”/ 
and d*p with the methods of §8 of IT. 


§7. THE ELECTROSTATIC INTERACTION 
BETWEEN 4d", d*"'s AND d*-*s? 


The electrostatic interaction between 
and s?'!S is given by 


(d? 1S|e2/r|s? 1S) 
= R*(dd, ss)(d? 'S| P2(cos w)|s?'S), (72) 


where R? is defined by TAS 8°8 and w is the angle 
between the radii vectors of the two electrons. 
From (51)II we have 


(2||C||0) =1, (73) 
and hence from (45)II and (38)II we get 
1S| P2(cos w) | s* 1S) =1/5!; 
since 
R*(dd, ss) = R*(ds, sd) = G?(ds) = 5G,(ds), 
(72) becomes 
(d? e?/r|s? 1S) (74) 


Introducing this result in (33c) and owing to 


if 
a 
a 
a 
oO" 
a 
a 
m 
bi 
| 


he 
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TABLE XVII. (d*vSL||70V% |\dv’S'L’). 


#P —19(14) —28 0 8(35) —8(14)+ —22(14) “0 0 
28 14(35)# 0 —28(10)4 56 —28(10)4 
YD 0 0 35(6)* 0 0 0 0 0 
—8(35) — 28(10) 0 —5(6)* —4(210) 4(210) —60(2) 0 
—8(14)* 0 4(210)! -—77 —98 3(35)* —4(385) 
22(14) — 28(10)+ 0 4(210)3 98 — 14(10) — 18(35) 4(385)4 
2G 0 0 0 —60(2)+ —3(35)* — 18(35)* 12(77) 
2H 0 0 0 0 —4(385)? —4(385)* —12(77) — (2002)! 
TABLE XVIII. 
as 0 0 0 4(210)4 0 0 0 0 0 0 
0 6(14)t —15(35)* —3(35) —105 —12(35)4 12(14)* 0 0 
0 —15(35)! 0 5(6)* 0 0 0 0 60(2)* 0 0 
—4(210)4 3(35) —s(6)t 15¢(210) —10(21)4 9(210)4 6(165)* 15(2)4 0 0 
0 —105 0 —15¢(210)t —%5¢(30)¢ 0 35(6)4 0 15(70)* 0 0 
oF 0 12(35)* —10(21)4 0 0 —21(10)* 0 12(14) —6(154)4 0 
oF 0 12(14)$ 0 —9(210)? 35(6)$ 21410) —42 —12(35)* 6(385)! 
0 0 0 —6(165)* 0 0 0 0 27(10)* 0 
eG 0 0 —60(2)# 15(2)¢ —15(70)* 12(14)¢ 12(35) —27(10)+ —6(33) —12(91)4 
eH 0 0 0 0 0 6(154)4 6(385)* 6(210)* 6(77) —3(2002) —7(39)¢ 
-0 0 0 0 0 0 0 12(91) —7(39)4 0 
TABLE XIX. (d®5SL||70V% 

as &S “D oF 8G 
FAY 0 0 —4(210)! —6(105)! 0 0 0 0 
oS 0 0 0 70(3) 0 0 0 0 
—4(210)! 0 15(6)* 60(3)+ —20(21) -4 — 20(15)8 0 
6(105)* —70(3)+ —60(3)! 10(15)+ 0 8(330)* —40(3) 0 
oF 0 0 20(21) 0 —105 (1155)? 14(105)! 0 
2G 0 0 —4(165) —8(330)! —(1155)# 1254,(33)4 —10(3)+ 84 , (30030)! 
s'G 0 0 20(15)! —40(3)4 14(105) 10(3)! — 10(330)+ *3(2730)8 
gl 0 0 0 0 0 , (30030)! 2(2730) —%54,(858)* 


(49’) we obtain 
(d"vSL | e?/r;;|d"~*s*v'SL) 
=[Q(n, (75) 


owing to the conventions of subsection (5) of §6, 
a minus sign must be introduced for n=6, v=2 
and for n=7, v=3. 

The calculation of the interaction between the 
configurations d" and d*~'s by means of (33b) is 
easy only for n=3, since in this case 


(d, SL}d*vSL) 


may be obtained from Table II by means of (29) ; 
but, as it was already mentioned at the end of 
§5, the calculations for 24 become very long, 


and it appears more convenient to use the fol- 
lowing method. 
The interaction in question is given by 


(d"vSL| 
= R*(dd, ds) 
-(d*vSL | (76) 


owing to (45)II and to the fact that >>; C;®” is 
a scalar and its non-diagonal matrix components 
vanish, we have 


(d"vSL|¥ P2(cos 
i<j 
CP (77) 


/1) 
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re- 
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TABLE XX. (d*vSL| 


TapLe XXI. (d'vSL| 


a dts A: a ds He da d's 
#P)?P 3(35)4 2P (70)4 —15 
0 2P 0 2F —2(5)4 
#D 2D)*D 39 (30)4 (3*P)*P 2(70)4 oF (#F)3F 3(5)4 
(3° F)?F —3(10)4 2D  (2D)'D — (105)# 
(3F)*F 2D (#D)!3D —3(5)! 2G (#G)'G 106 (3)4 
2G (2'G)?G —5(2)8 (3*D)'D 6(10)4 @G (G)'G 5g (33)4 
#D (#D)3D —4(5)* 3(5)4 


and using (33)II we obtain 
2(2L+1)(d"vSL| P2(cos 
i<j 
= (-1)r" 


-(d*vSL||X 
SL" ||; sSL) 
+ 


From (27), (44)II, (80)II, and (73) we obtain 
finally 


(d*vSL| 

=(n/14)C 

(d*vSL|| SL") 

SL'"[d"-1(v' S’L)d SL) 
*(2L"+1)4/(2L +1) + 2 

(2L-+1)*]R*(dd, ds). (79) 


By means of this formula the interaction 
between the configurations d* and d*~'s was cal- 
culated for n=3, 4, 5; the results are different 


a He as dts Aa 
8(S)4 vF (3F)?F —2(15)4 
(°P)?P (210) —39(15)4 
(8P)*P 0 8F (8F)F 154(3)4 
(!D)?D — (35)# (3F)*F — (30)4 
#D (:!D)?D 3(S)4 — 194 (3)! 
#D (¢D)?D 3(10)4 56 (33)! 
(#D)?D 2(15)4 vG (15)§ 

2D (¢D)?D — (5)! 8G (¢G)°G — %(5 
(8D)?D —3(30)! 8G (8GPG —5e(1 
sD (8G)'G 5(2) 
sD —59(14)! 2H (¢H)*H 15)! 
(FPP el (S) 


from zero only if v’=v+1 and are given in 
Tables XX—X where the quantity 


H,= R*(dd, ds)/35 (80) 


was assumed as parameter. For n=3 our results 
agree with those given by Marvin." _ 

The interaction between the configurations 
d"~'s and d"~*s? may be calculated in the same 
way ; the result is 


e?/r;;|d"-*s*vSL) 
2S’+1 ] 


=(-1 


2S+1 


(81) 
"H. H. Marvin, Phys. Rev. 47, 521 (1935). 
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On the Scattering of Neutrons by Protons 


LAMEK HULTHEN 


Institute for Mechanics and Mathematical Physics, 
niversity of Lund, Lund, Sweden 


February 11, 1943 


ECENT experiments by Amaldi, Bocciarelli, Ferretti, 

and Trabacchi' give valuable information on the angu- 
lar distribution of fast neutrons scattered by protons. 
Measuring the scattered intensity J(@) in two directions, 
6=x/2 and x(@=scattering angle of the neutron in the 
center of gravity system), they found values of I(x) /I(#/2) 
between 0.71 and 0.52 for average energies of the incident 
neutrons varying from 12.5 to 14 Mev. 

We know, that “symmetrical’’ and ‘‘neutral’”’ theories 
must give rather different predictions on this angular dis- 
tribution.? From the symmetrical theory expounded by 
Mgller and Rosenfeld* we obtain for 14 Mev neutrons, 
considering the effect of S, P, and D states, a total scatter- 
ing cross section =0.66 X 10-* and 


I(@)~1—0.18 cos 6+0.45 cos? (1) 


thus I(x) /I(x/2) =1.63. This result can hardly be recon- 
ciled with the experience of Amaldi and his collaborators. 

Because the convergence of the power series in cos @ is 
somewhat slow, it seems necessary to estimate the influence 
of the neutron-proton states with greater angular momenta. 
Using a Born approximation, we get, assuming again the 
static interaction V, of Mgller and Rosenfeld:* 


where k= (ME;/2h*)*, and E; is the energy of the incident 
neutron. For E;=14 Mev and a meson mass= 200 m,, we 
have k=0.8x. Equation (2) then gives ¢=0.53 X 10-* cm? 
and I(x)/I(x/2)=4.8. The discrepancy mentioned above 
thus becomes still more striking. Besides, the comparison 
between (1) and (2) indicates that the usually accepted 
approximation, if we take only S and P states into con- 
sideration,* may be too rough for these high energies (k 
nearly equal to «). 

The effect of the non-static potential W,” has also been 
taken into account. It does change the angular distribution 
in the “‘right’’ direction, but the correction is too small to 
be significant. _ 

Comparing with theory, Amaldi, Bocciarelli, Ferretti, 
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and Trabacchi find that their results agree with the angular 
distribution deduced by Ferretti* from the neutral theory 
of Bethe,* but are inconsistent with the outcome of the sym- 
metrical theory according to Rarita, Schwinger, and Nye.* 
Still it should be noted that the distribution calculated by 
the latter authors from another special form of the neutral 
theory* gives practically the same value of J(x)/I(x/2) 
as that of Ferretti (because 1+) cos @=1+(b+<c) cos @ 
+c cos? 6 for @=x/2 and x, irrespective of c). 

It should also be observed that noncentral forces are not 
necessary to explain the effect considered here. In fact, the 
simple potential 

Vir) = —[Gi+G2(oo™) Je" /r (3) 


ensuing from a ‘“‘mixed”’ neutral theory,'® gives, to the 
same approximation as (2): 


(6) = (M/h*)*(3G2+G:*) sin? (0/2) P. (4) 


For k=0.8« this gives I(x) /I(x/2) =0.41. 

Thus the experiment of Amaldi and his collaborators 
seems to speak in favor of a neutral meson theory of nuclear 
forces and against the symmetrical possibility. This ap- 
pears to be the first really serious difficulty which the the- 
ory of Mgller and Rosenfeld has met. 

It is, of course, not possible to draw more extensive con- 
clusions until more detailed evidence of the angular dis- 
tribution is obtained. In this connection the importance of 
analogous experiments on the photo-disintegration of the 
deuteron? should again be stressed. 


1 E. Amaldi, D. Bocciarelli, B. Ferretti, and G. C. Trabacchi, Natur- 
wig 90,380 (1942) d H. A. Nye, Phys. Rev. 59, 209 
an ye, Phys. 
eld kindly called my attention to an earlier 
discussion of the sub: = E. Boeke f. 107, 683 (1937). 
Méller and L. densk. Selsk. math. -fys. 
Medd. 17, U 
4Dr. C. Mol ndly informed me that a similar calculation has been 
made by B. Ferret, Particulars are not known here y 


*C. Kittel and %: “Breit Phys. Rev. 56, 744 (1939); Rarita, Schwinger, 
and Nye, reference 2 

7 Reference 3, p. 42. 

§ B. Ferretti, Ricerca Scient. 12, 843, 993 (1941). 

*H. A. Bethe, Phys. Rev. 57, 261, 390 (1940). 

1° Reference 3, p. 34. 


Energies for Self-Diffusion in Zinc 


H. B. HUNTINGTON 
Washington University, Saint Louis, Missouri 
April 17, 1943 


XPERIMENTS by Miller and Banks! on self-diffusion 
in zinc have shown that the diffusion constants for 
flow in the basal plane and along the hexagonal axis have 
widely different temperature dependence. If one expresses 
the diffusion constant in the form Ae~@/®7, then in the 
basal plane A = 8X 10° cm*/day and Q = 3.1 X 10* cal./mole, 
and along the hexagonal axis A =4X10* cm*/day and 
Q=2.04X10* cal./mole. Accordingly diffusion is always 
more rapid along the axis because of the smaller activation 
energy involved. The anisotropy of the temperature- 
independent factor has relatively less effect. Some geo- 
metrical considerations have already been put forward to 
explain the difference between the A values on the basis of 
diffusion by vacancies.? 


The fact that Qpasai is larger than Qaxiai might be ex- 
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pected since the Debye temperature for lattice vibrations 
in the basal plane is higher than for the vibrations trans- 
verse to these planes. If the hypothesis of vacancy diffusion 
is correct, then it would appear that the diffusing atom 
halfway between two vacancies in the same basal plane is 
in a position of higher energy than a similar atom halfway 
between vacancies in adjoining planes. A comparison of the 
electrostatic energy of the two situations has been made on 
the assumption that the two valence electrons per atom 
are uniformly distributed in space. For a close-packed hex- 
agonal lattice the migrating atom would have in each case 
four nearest neighbors at a distance +/} times the lattice 
constant and the difference in energy would be small. The 
zinc lattice, however, departs considerably from close- 
packing (c/a is 1.86). The calculations for this case were 
carried out by the Ewald method* and showed that the 
configuration corresponding to basal diffusion lies higher 
by 0.055(2e)?/a or 2.7X10* cal./mole. This is about two 
and a half times the difference between the Q values. It 
might be expected that such simplified considerations would 
give too large a result, first because the two conduction 
electrons in metallic zinc probably behave quite differently 
from free electrons, and secondly because polarization 
effects will tend to reduce all electrostatic terms and differ- 
ences between such terms. 

1P. H. Miller, Jr., and Floyd R. Banks, Phys. Rev. 61, 648 (1942). 


by: s. Rev. 61, 325 (1942). 
Physik 64, 253 (1921). 


Erratum: Simple Method for the Investigation of 
Secondary Electrons Excited by y-Rays and 
the Interference of These Electrons with 
Measurements of Primary 6-Ray Spectra 


MEITNER 
Forskningsi: t for Fysik, Stockholm, Sweden 
” (Phys. Rev. 63, 73 (1943)) 


N page 76, first column, the lines which read “‘for the 
primary §-particles a range of 0.36-mm Al corre- 
sponding to 0.098 g/cm?,”’ should be changed to read ‘‘for 
the primary 8-particles a range of 0.36 mm corresponding 
to 0.13 g/cm?.” 


Mesotron Mass and Heavy Tracks on Mt. Evans 


Car. E. NIELSEN AND WILSON M. PoweELL* 
Department of Physics, University of California, Berkeley, California 
April 30, 1943 


BOUT ten fairly reliable observations of the mesotron 
mass have been published.' Most of these are con- 
sistent with a mass approximately 180 times the electron 
mass, but they do not constitute conclusive evidence for 
uniqueness of the mass. It is thus clearly desirable that 
more observations of the mass be obtained, both to deter- 
mine whether it is unique and to fix its value more pre- 
cisely. 
An expedition to Mt. Evans was made for the primary 
purpose of obtaining such observations. The ionization and 
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momentum of cosmic rays were observed simultaneously 
in a Wilson cloud chamber by photographing tracks curved 
by a magnetic field of 2500 oersted, and diffused sufficiently 
to permit photographic resolution of single droplets. A de- 
lay of ~0.15 second between reduction of the clearing 
field from ~20 v/cm to ~2 v/cm and expansion of the 
chamber was provided. Most of the photographs were of 
random expansions, and this time of reduced clearing field 
corresponds to the interval before expansion during which 
a track remains to be observed. It has been shown®® that 
condensation upon even a small fraction of the negative 
ions indicates a supersaturation sufficient to give condensa- 
tion upon substantially all positive ions; and this result is 
valid for any water-ethyl alcohol mixture in a broad 
range around the 70 percent alcohol mixture requiring 
minimum expansion. The number of droplets formed in a 
positive track with an associated negative is therefore 
equal to the number of ion pairs. An enlarged section of a 
typical mesotron track pair is reproduced in Fig. 1. 

Mass determination from ionization and curvature de- 
pends upon a knowledge of the relation between ionization 
and speed. The generally accepted form of this relation is 
expressed in the theoretical equation 

A 

This is the function used by Corson and Brode‘ to obtain 
mass from ionization and curvature. The experiments of 
C. T. R. Wilson and E. J. Williams® have yielded an experi- 
mental relation based upon counts of primary ionization of 
slow electrons. Their work is summarized by the equation 
I =1,8"'-*, which has been used by Williams and Pickup to 
derive mass values. For comparison it is interesting to 
note that the theoretical equation is well represented in the 
interval 0.20<8<0.90 by J =J98"'-*. Mass values derived 
from the experimental curve are larger than values from 
the theoretical curve. This uncertainty in interpretation 
would be removed and at the same time the possibility of 
error from inaccurate ionization measurement would be 
eliminated if one were to observe under identical conditions 
the curvature and ionization of slow electrons of the same 
speed as the mesotrons. This procedure would permit an 
unambiguous mass assignment: Two particles with equal 
charge ionize equally when moving at the same speed, 
and their masses are then proportional to their momenta 
or Hp values. Experiments to determine the ionization of 
slow electrons are in progress, but until their completion 
the mass values are subject to a possible error. For the 
present, computations from the Corson and Brode nomo- 
graph and the Williams curve are both given. 

Approximately 7000 photographs were taken on Mt. 
Evans, 14,100 ft. and 1800 were taken at Summit Lake, 
12,700 ft. The photographs show numerous heavily ion- 
izing tracks with negligible curvature produced by particles 
—presumably protons—much heavier than mesotrons. 
These heavy tracks were present in about eight percent of 
the pictures at 14,100 ft. and in about four percent of the 
pictures at 12,700 ft. Each picture records events in an 
interval ~0.3 sec., and in a volume ~1500 cm®. Only six 
heavily ionizing mesotrons have yet been identified. Since 
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Fic. 1. Section of mesotron track 6949, enlarged to about three times 
its size in the cloud chamber. Ion pairs are separated into positive and 
negative columns, with the positive on the left side. At the top a sharp 
electron track crosses the picture. 


they are quite conspicuous it is probable that not many 
more will be found by additional study of the pictures. 
Most of the heavy tracks thus do not represent slow 
mesotrons ;* in fact, the fraction of slow mesotrons may 
not greatly exceed the figure of 0.001 found by E. J. 
Williams at sea level. 

Four of these slow mesotrons yield mass values as 
follows: 


Track Ionization Hp Mass: C. & B. Mass: William’ 
8-21, 2 2.5 XIo 2.26 230+20 
1621 6.0 1.08 80+20 225420 
6941 6.2 1.11 240215 
6949 2.4 1.6 145+30 155+30 


The probable errors are estimated to include chamber dis- 
tortion of curvature and statistical uncertainty deriving 
from the finite number of droplets counted. Curvature 
was measured by plotting micrometer readings of coordi- 
nates along the track and fitting a curve to the points by 
the method of least squares. The curvature of the first 
track, for example, was determined from 13 points, with a 
probable error in fitting of 3.5 percent. The ionization of 
this track was determined by a count of 335 droplets, 
hence the statistical error is ~5 percent (or more, since 
each ion is not the result of a single independent event). 
We are grateful to all who have assisted with this experi- 


ment, and to the Carnegie Institution of Washington, the 
Fund for Astrophysical Research, and the Rumford Fund 
of the American Academy of Arts and Sciences for generous 
financial support. 
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3C. E. Nielsen, Phys. Rev. 61, "302A ( 
«D. R. Corson and R. B. Brode, Phys. Rev. -. 773 (1938). 
5 E. J. Williams, Proc. Roy. Soe. AIZZ, 194 (19 39). 
* Compare results of W. M. Powell, Phys. Rev. 61, 670 (1942). 


X-Ray Line Broadening by Cold-Working 
Alpha-Brass 
MAX PETERSEN AND C. W. TUCKER 


Lehigh University, Bethlehem, Pennsylvania 
May 3, 1943 


IEMANN and Stephenson have reported' that the 

x-ray diffraction line broadening produced by small 
amounts of cold-working of alpha-brass is enhanced by 
speed of working, while the internal friction thus intro- 
duced is not. The first of these facts, which might be of 
great importance to practical metal forming, is not con- 
firmed by our observations. 

Three like-tensile specimens of cartridge brass having 
sections 6.3 by 0.53 mm, sealed in nitrogen, were annealed 
for ninety minutes at 400°C, and showed fine-grained 
strain-free diffraction by back-reflection. Each was then 
cold-stretched about 7.5 percent, one by a falling weight 
at a speed exceeding 30 ft./sec., one by a hammer-blow of 
somewhat less speed, and one slowly in a small tensile 
stressing frame. The back-reflection patterns from these 
three strained specimens show no significant difference of 
line broadening. In each the Ni Ka doublet has become 
nearly irresolvable visually. If any difference can be found 
it is that both fast-worked specimens gave traces of coarse- 
grain irregularity in the photographs, as the other did not. 

These observations seem to remove the disparity be- 
tween line broadening and internal friction effects of cold 
work. The earlier observations may have been effected by 
initial differences between the specimens. 


1F, Niemann and S. T. Stephenson, Phys. Rev. 62, 330 (1942). 


The Forces Between Hydrogen Molecules 


HENRY MARGENAU 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 
May 1, 1943 


HE difficult part in an a priori calculation of inter- 
molecular forces is the determination of the repulsive 
exchange forces. For this reason, complete interaction 
curves have hitherto been obtained for a few simple 
spherical atoms only. Yet to understand even the most 
elementary properties of molecules, such as their size and 
shape, and also to locate the position of the minimum of 
their van der Waals forces, a knowledge of the exchange 
forces is required. Attracted by the fundamental nature of 
the problem, I have undertaken to calculate them for the 
simplest molecule and herewith report the results. 
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24 286 32 36 40 
Fic. 1. Interaction energies of two hydrogen molecules, in millivolts, 
plotted against distance between their centers (in ®: Relative orienta- 
tion of molecules is as follows: (a) ——, (b) — |, (c) | |. 


A function of the Heitler-London-Wang type was chosen 
to represent the electronic charge distribution of a single 
molecule. As is well known, this function yields the best 
molecular energy when the effective nuclear charge, Ze, is 
taken to be 1.166 electronic units.' In the present calcula- 
tion, Z was at first left undetermined. 

The interaction was investigated for three relative posi- 
tions of the molecules: (a) both molecular axes are parallel 
to the line connecting the centers; (b) one molecule is 
parallel, one perpendicular; (c) both molecules are per- 
pendicular to the center line. The customary approximation 
of neglecting double exchange integrals leads to results 
which are entirely erroneous in this problem, and it is 
necessary to retain all terms. Unwieldy four-particle 
integrals, however, can be reduced with good validity to 
known three- and even two-particle integrals. 

One feature emerging from the calculation is the fol- 
lowing. When Z is taken to be 1.166 (Wang’s value), the 
exchange forces are attractive. This emphasizes the fact that 
the effective nuclear charge does not have a unique physical 
significance and that its constancy may not be assumed in 
molecular problems where the charge symmetry is far from 

spherical. On further consideration it turned out that in 
the other two-electron problem, the interaction of He 
atoms, there is a limiting value of Z, namely (11/8)e, 
below which the exchange forces are also attractive in the 
interesting range of atomic separations. But in this case 
the normal value of Z is great enough to cause no difficulty. 
The physical reason for the inadequacy of the customary 
value of Z in the intermolecular problem is to be found in 
the circumstance that the electronic charge is strongly 
concentrated between the two nuclei of a molecule, leaving 
the nuclei bare to an abnormal degree. 

The calculations were therefore carried out for several 
different values of Z. Various considerations lead to Z ~ 1.4 
as the most reasonable choice. Hence the results will here 
be given for that case only. 

An H:-molecule represents a static quadrupole. The part 
of the interaction arising from this fact is almost negligibly 
small throughout. It can be expressed in the simple form 
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where A is the usual overlap integral between the atoms of 
a single molecule and d the distance between these atoms, 
The function f expresses the characteristic angular varia- 
tion for forces between linear quadrupoles. 

Finally the attractive van der Waals energy must be 
included. This may be taken from the work of Massey and 
Buckingham,? whose result, however, is larger than can be 
obtained from dispersion data on Hz vapor. When the 
dipole-quadrupole forces are added, the result is in fair 
agreement with that derived by the latter authors. 

Composition of all these effects yields the interaction 
curves (drawn for cases (a), (b), and (c) as defined above) 
which are given in Fig. 1. As would be expected, collinear 
molecules repel at larger distances and produce a shallower 
minimum than parallel ones. When the distances at which 
the minima occur, and also their depths, are compared 
with the empirical curves derived by Lennard-Jones and 
others’ from gas kinetic effects, satisfactory agreement is 
found. 


1S, C. Wang, Phys. Rev. 31, 579 (1928). 
31 1538). Massey and R. A. Buckingham, Proc. Roy. Ir. Acad. 45, 


Fowler and E. A. Guggenheim, Statistical Thermodynamics 
(Canibridge University Press). 


The Disintegration Scheme of Na™ 


L. G. Ettiott, M. Deutscu, AND A. ROBERTS* 
Massachusetts Institute of Technology, Cambridge, Massachusetts 
April 30, 1943 


HE gamma-rays of Na* have been studied by several 
authors. Various energies between 0.8 and 3 Mev 

were announced at different times. In 1941 Itoh! reported 
finding only two gamma-rays of 1.38- and 2.8-Mev energy 
respectively and equal intensity. Our experiments, com- 
pleted before we received Itoh’s paper, confirmed his 
results and therefore were reported only briefly at a section 
meeting of the Physical Society.? Since then these results 
were again questioned by C. E. Mandeville.* We have 
since repeated the measurements with improved accuracy. 


He 
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a (dotted groups the 1 38 and 2.76 
Mev gamma-rays of Na™ 
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Figure 1 shows the momentum distribution of secondary 
electrons ejected from a thin lead radiator (and its brass 
backing) observed in the magnetic lens spectrometer. The 
observed counting rates were not divided by Hp as is 
customary. This procedure was adopted to keep the ordi- 
nates of the curve more nearly equal throughout the entire 
range. Besides, we wish to emphasize the two photoelectron 
lines (marked by arrows) whose heights should not be 
divided by the transmitted momentum interval since they 
represent nearly monokinetic electron groups. 

The half-width of the transmitted momentum band is 
about 3.5 percent, but the observed width of the photo- 
electron lines is slightly greater because of the thickness 
of the radiator and the presence of the unresolved L photo- 
electron groups. The statistical counting errors are negli- 
gible for most of the points. However the probable error, 
as indicated by the scatter of the points is about two per- 
cent, due partly to small errors in correcting for decay of 
the sources and in normalizing the several runs of which 
the curve is a composite. 

All regions where photoelectron lines were suspected 
were studied in separate runs, and it seems quite certain 
that the intensity of any other gamma-ray above 0.34 Mev 
(if present at all) does not exceed ten percent of the in- 
tensity of the two observed lines. 

The distribution of Compton recoil electrons from a 
copper radiator was also studied. Its absolute intensity and 
shape deviate by several percent from the Compton groups 
from lead, probably because of the difference in scattering 
of the secondary electrons in the two radiators. Therefore 
only the regions near the end points are indicated in Fig. 1 
(crosses and dotted line). 

The energies of the gamma-rays, obtained both from the 
photoelectron lines and the Compton end points are 1.38 
+0.03 and 2.76+0.06 Mev. The ratio of the intensities 
may be obtained by correcting the heights of the lines for 
the variation in photoelectric cross sections and retardation 
in the radiator. Details of this correction will be published 
soon in a paper dealing with spectrometer technique. The 
ratio thus obtained is J,., : J2..=0.9 : 1, uncertain to about 
fifteen percent. 

While some of the divergent results previously obtained 
by other authors were probably due to misinterpretation of 
data, others may have been due to impurities. The remote 
possibility cannot be excluded entirely that there may be 
several gamma-rays of very low intensity in the energy 
range between the two intense lines. 

Coincidence measurements of beta- and gamma-rays and 
of gamma-rays with each other were also performed, and 
the results agreed with those of other workers," * showing 
that each beta-ray is followed by a 1.38- and a 2.76-Mev 
gamma-ray in cascade. Mr. W. C. Peacock of this labora- 
tory has recently repeated these experiments with carefully 
calibrated gamma-ray counters, and it appears quite cer- 
tain that the great majority (at least 90 percent) of the 
disintegrations proceeds by this scheme. As in the case 
of Y** the fact that the gamma-ray energies are in the 
ratio 1 : 2 has therefore no evident significance in terms 
of the disintegration scheme. 


Our thanks are due to Professors Evans and Livingston 
and the cyclotron crew. 


* Now at the M. I. T. Radiation Laboratory. 

as Itoh, Proc. Phys. Math. Soc. Japan 23, 605 (1941). 

2L. G. Elliott et al., Phys. Rev. 61, 99A (1942). 

3C. E. Mandeville, Phys. Rev. 62, 309 . 

«Feather and Dunworth, Proc. Camb. Phil. . 34, 442 (1938). 


Gamma-Rays from and 


C. E. MANDEVILLE 
The Rice Institute, Houston, Texas 
April 30, 1943 


EASUREMENTS on the energies and relative inten- 

sities of the y-rays from Na* have been repeated 

with a magnetic spectrograph which has been previously 
described. Recent measurements! seemed to indicate the 


2 4 6 


Fic. 1. Momentum distribution of the Compton recoils 
of the y-rays from Na®™. 


presence of quanta at 0.84, 1.31, 1.66, and 2.90 Mev. 
Subsequent improvements and modifications in the opera- 
tion of the spectrograph and experiments** on the y-rays 
from Sc and As’* have disclosed that the spectrograph 
was not satisfactorily designed for the measurement of 
y-rays of energy less than 1.0 Mev and with an intensity 
less than fifty percent of that of any high energy y-rays 
present when the first experiments! on the gamma-radiation 
from Na* were performed. In addition to measurements of 
the energies of the -rays from Sc* and As”*, observations 
as yet unpublished, extending over a wide range of energy, 
have been made upon a considerable number of the radio- 
active isotopes of other elements. On plotting the distribu- 
tions in momentum of the Compton recoils of the y-rays 
from these various radioactive isotopes, it became apparent 
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Fic. 2. Momentum distribution of the Compton recoils of the y-rays 


from La™“*, The symmetry of the distribution about the broken line ig, 


to be noted. 


that for a single y-ray, and for the geometry of this 
spectrograph, the momentum distribution of the Compton 
recoils is symmetrical about a line through the midpoint 
of the peak of the distribution and parallel to the N/Hp- 
axis. With this fact in mind, it soon became clear that the 
momentum distribution of the Compton recoils of the 
2.9-Mev y-ray from Na*‘ lacked the above-mentioned sym- 
metry; that is, the Hp coordinate of the peak of the dis- 
tribution appeared to be too large. 

NaF has been irradiated by slow neutrons from the 
cyclotron at Washington University, St. Louis, Missouri. 
The momentum distribution of the Compton recoils of the 
y-rays from Na* is given in Fig. 1. It is to be noted that 
the spectrum is composed of two electron groups corre- 
sponding to quantum energies of 1.38+0.03 Mev and 
2.94+0.06 Mev, respectively. The intensity ratio is 0.84 
and therefore not inconsistent with the view that the 
quanta are in cascade. Measurements with this spectro- 
graph are now in agreement with those of J. Itoh‘ and 
Elliott, Deutsch, and Roberts.* The distortion of the peak 
of the 2.9-Mev y-ray and the presence of what appeared 
to be the peak of a y-ray at 1.66 Mev reported in the first 
experiments on Na*‘ may be attributed to the presence of 
an impurity unknown in character, perhaps related to 
deuteron bombardments. A similar explanation might 
apply to the strong 2.04-Mev line reported by Richardson.* 

The symmetry of the Compton distribution of a single 
y-ray is illustrated in Fig. 2 where is presented the distri- 
bution in Hp of the Compton recoils of the y-rays from 
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La’, The y-rays appear to be monochromatic, and their 
energy computed from the end point of the momentum 
distribution is 2.04+0.04 Mev, in surprisingly good agree- 
ment with the value 2.00+0.05 Mev obtained by Weimer, 
Pool, and Kurbatov,’ using lead absorption. 

It was not at first realized that the peak of the electrons 
from the 2.9-Mev y-ray of Na*‘ had been misplaced, since 
its position seemed to agree very well with the position 
suggested by what appeared to be the peak of the 2.62-Meyv 
y-ray from Th(C+C”’). However, several y-rays of lower 
energy accompany the 2.62-Mev y-ray so that the exact 
position of its peak may be obscured somewhat by their 
presence. An examination of the curves for Th(C+C”) 
and Na* given by Curran, Dee, and Strothers® will reveal 
that the 2.00-Mev radiation from Na* reported by them 
may actually be part of the peak of a symmetrical mo- 
mentum distribution of Compton recoils from their 3.03- 
Mev radiation. The spectrograph of the writer is very 
similar in geometry and design to the one employed by 
them. 

One conclusion of this research is therefore that the end 
point of the momentum distribution of the recoils of the 
2.62 Mev y-ray from Th(C+C’’) may be used for cali- 
brating spectrographs of this type, but no particular sig- 
nificance should be attached to the shape of the curve. 

Shortly after completion of this repetition of the meas- 
urements on Na*‘, the writer was informed of the results of 
a further study of the y-rays from Na*‘ which has been 
conducted by Elliott, Deutsch, and Roberts.*® Their newly 
performed spectrometer and coincidence experiments indi- 
cate that the two quanta are in cascade and that the 
8-ray spectrum of Na*‘ is simple. 

The theoretical status of Na** appears to be very un- 
satisfactory. Objections to the cascade arrangement of the 
two y-rays have been discussed by Guthrie and Sachs," 
and possible complexity of the Na** 8-ray spectrum has 
been suggested by the recent results of Konopinski and 
Uhlenbeck." 

The writer wishes to acknowledge the energetic coopera- 
tion of the members of the cyclotron group at Washington 
University, St. Louis, Missouri. They prepared and sent 
to him the radioactive materials used in the course of these 
experiments. 
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3C. E. Mandeville, Phys. Rev. 63, 91 (1943). 

Phys. Math. Soc. Japan 23, 605 (1941). 
iott, M. Deutsch, and A. Roberts, or Rev. 61, 99 (1942). 
‘I; Reginald Richardson, Phys. Rev. 53, 124 (19. 
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